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ABSTRACT 


A large gravity gradiometer flown in a low earth orbit 
would significantly improve man's knowledge of the earth's gravitational 
field. Such knowledge would be of value to the geodesy, geophysics and 
geology fields, and would have application to orbital mechanics, naviga- 
tion, guidance, earth dynamics, and mineral prospecting. This report 
describes some preliminary mission studies and the design, fabrica- 
tion, and test of a breadboard model of an earth orbital, rotating 
gravity gradiometer with a design goal of 1 0 ^ sec”^ (0. 01 EU) in a 
35-sec integration time. 

The proposed mission uses a Scout vehicle to launch one (or two 
orthogonally oriented) spin- stabilized satellites into a 330-km circular 
polar orbit some 20 days before an equinox. During the short orbital 
lifetime, the experiment would obtain two complete maps of the gravity 
gradient field with a resolution approaching 270 km (degree 75). Indi- 
vidual point anomalies, although smaller than 270 km, would be mea- 
sured to 0.01 EU. The higher order harmonics would be sampled inde- 
pendently over 100 times during the mission to give an amplitude accur- 
acy after data reduction below 0.001 EU . 

The breadboard model of the gradiometer demonstrated a com- 
bined thermal and electronic noise threshold of 0.015 EU per data 
channel. Hughes Research Laboratories (HRL) identified the design 
changes needed to reduce the noise to less than 0.01 EU . Variations 
of the sensor output signal with temperature were experimentally 
determined and a suitable method of temperature compensation was 
developed and tested. Other possible error sources, such as sensor 
interaction with satellite dynamics and magnetic fields, were studied 
analytically and shown to be small. 



SEC TION I 


INTRODUCTION 


This final report co\tains the results of a 15-month, $110,000 
contract for the design, fabrication, and test of a rotating gravity 
gradiometer for earth orbit applications. The work was carried out 
under the ai spices of the NASA Advanced Applications Flight Experi- 
ments Program Office. Langley Research Center, Hampton, Virginia. 
The objective of the AAFE program ic to carry forward instrument 
development independent of flight opportunities. 

On this contract HRL developed a preliminary design of an 
experiment for measurement of the gravity gradients of the earth's 
field from orbit, and determined a set of mission, spacecraft, and 
sensor parameters that would achieve the desired scientific goals. 

HRL then executed a detailed design of a sensor structure, fabricated 
a breadboard model of the sensor structure and the important parts of 
the data processing electronics, and conducted tests of the sensor 
performance. 

Because of the relatively low level of effort, the testing program 
was limited to simulated signals. However, past testing with real 
gravity gradient signals (see Appendix B and Attachment A) had shown 
the validity of these simulations. The stated design goal of the con- 
tract was a sensor with a threshold sensitivity of 0. 1 EU (1 EU 

-9 -9 -2 -12 

= 10 gals/cm = 10 sec * 10 g/cm) with a 35-sec integra- 

tion time. As a result of preliminary design studies, HRL determined 
that a design goal of 0 01 EU would m-oduce a sensor that obtains 
useful scientific return from the experiment. A sensor was fabricated, 
having parameters that would result in a noise level of 0. 015 EU, and 
modifications that would result in a sensor with a noise level approach- 
ing 0. 007 EU were identified. 
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Prior testing on feasibility models of the rotating gravity 
gradiometer (1969 through 1970) have indicated that near thermal noise 
limited threshold signal levels could be attained with the sensor operat- 
ing in a sufficiently quiet (simulated free-fall) environment. Because 
of the low resonant frequency and large size of the earth orbital model, 
the cost of attaining such an environment was considered beyond the 
scope of the contract. Such a test could be run, and provided the rec- 
ommended improvements made in Section V of this report were com- 
pleted, such a test should result in demonstrating a threshold noise 
level near the sensor thermal noise limit (<0. 02 EU}* The estimated 
excess noise over the combined thermal and electronic noise of the sen- 
sor (<0.01 EU) would be the result of the residual sensor interaction 
with the laboratory vibrational environment. The foregoing estimate 
assumes a test involving a carefully isolated, magnetically shielded, 
and thermally controlled sensor, not rotating, with nearby rotating 
masses providing a dynamically varying gravitational gradient field. 

Because of the high vibration levels, concomitant with rotating 
this large size gradiometer in an earth-bound environment, an attempt 
to conduct a threshold sensitivity test in which the gradiometer rotates 
in a static gradient field is considered unfavorable. Such tests could 
quite readily be conducted, however, in a free-fall environment such as 
Skylab. 
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SECTION II 


DESIGN STUDIES 


The contract began with a series of design studies to determine 
the relevant parameters of the gravity field to be measured, the various 
possible missions that could be considered, and the effects of these 
missions on the satellite and sensor design parameters. The design 
of the gravity gradiometer is strongly dependent upon the particular 
mission and using vehicle, a condition that is especially true for this 
particular application to earth geodesy. The size and operational 
parameters of the sensor are determined by the orbital altitude and 
inclination, mission lifetime, and measurement requirements of 
geodesy. The sensor, in turn, has an effect on the spacecraft, 
especially the requirements for spin speed, temperature, and attitude 
control. 

The presently envisioned experiment is the result of these 
design studies. It uses a Scout vehicle with a 42-in. diameter payload 
shroud to launch one (or two orthogonally oriented) spin-stabilized 
satellites into a 330-km circular polar orbit some 15 to 20 days before 
the vernal or autumnal equinox. The satellite would carry a 76 -cm 
diameter gravity gradiometer with a sensitivity of 0. 01 EU at 35-sec 
integration time. The orbital lifetime would be short, but during that 
time the experiment would obtain at least two complete maps of the 
gravity gradient field with a resolution approaching 270 km (540 -km 
wavelength or degree 7 5). Individual point anomalies, although smaller 
than 270 km in extent, would be measured to 0.0 l EU while the higher 
order harmonic components would be sampled independently up to 
170 times during each orbit, thereby giving an amplitude accuracy 
after data reduction down to 0. 0008 EU for the higher order harmonics. 

HRL conducted brief design studies early in 1971, that resulted 
in the proposed general experiment design, and were verified 
in detail by a separate Phase A study* carried out by Jet Propulsion 
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Laboratory (JP L) during late 1971 and ear'; f97Z. The reader is 
referred to this more recent and more cn;n p reshens ive re r < 
detailed discussion of the mission parameters and two dil,- rent 
satellite designs, including power, weight, subsystem, and data reduc- 
tion estimate s . 

A. GRAVITY FIELD PARAMETERS 

To determine the gravity gradient sensitivity requirements lor 
the sensor, URL performed certain calculations of the gravity gradient 
field strengths to be expected in earth orbit. At the time of its design 
studies, little published data existed. Prior to the contract, HRL had 
undertook various studies to estimate the gravity gradient fields at 
altitude. For its estimates HRL used the Kaula "rule of thumb”^ 
that the strength of the various harmonic orders goes as 

- ^ - 1 - s(c + S 2 ) « - 1 - 0 - 10 - n > 2 . 

2n + 1 \ nm nm / 4 

m n 


A summation of this HRL effort was presented at the Washington AGU 
meeting in April 1971, and the reprint is included as Attachment A. 

HRL re-examined this work in its design effort (see Appendix B) and 
found it adequate for an initial examination of the problem. A calcula- 
tion was made of the gravity gradients to be expected for typical single 
anomalies (Fig. 1) and for various harmonic orders (Fig. 2), and the 
conclusion was drawn that in order to make a significant contribution 
to geodesy and geophysics, the gradiometer sensitivity should be 
better than the 0.1 EU which was the stated design goal in the con- 
tract statement of work. As a result of these studies, HRL adopted 
0.01 EU * ^ -ne -’esign goal for the sensor. 

.hese preliminary studies early in 1971 were followed by a 
number of studies by others that generally confirmed HRL’ s preliminary 
studies. A pacer by Sandson and Strange of Computer Sciences 
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Fig. 1. Vertical Gravity Gradient of Single 300-km 
Diameter Disk at 250-km Altitude. 



Fig. 2. Vertical Gravity Gradient of Various Harmonic 
Orders at 250-km Altitude. 
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Corporation calculated the gradient at altitude from 1° x 1° surface 
gravity data, and a 12th degree and order satellite gravity field at 
300 km altitude. The results (Fig. 3) indicated that accuracies better 
than 0. 1 EU were required to provide us. ful improvements to present 
geophysical knowledge. 

4 

A paper by Chovitz, I as, and Morrison of NOAA calculated 
the gravity gradients along hypothetical 300 -km altitude orbits selected 
in regions of dense 1° coverage. These results verified the validity 
of the Kaula "rule of thumb" to 30%. Typical results of their computer 
simulation were a magnitude of 0. 008 EU caused by the harmonics 
of the 6 1st degree only, and of 0.02 EU caused by the combined 
harmonics of the 61st through the 70th degree. 

A paper by Glaser of JPL^, using ext >oiations of Kaula' s 
rule of thumb, obtained similar results (Fig. 4) and also compa ed the 
relative accuracy of doppler, altimeter, and gradiometer techniques 
of obtaining gravity field data from orbital satellites (Fig. 5). As 
expected, each technique has its region of applicability and the three 
techniques should be considered as complementary rather than as 
competitive techniques. 

B. MISSION PARAMETERS 

To map the higher order harmonics of the earth's gravity field, 
it would be desirable to have the measurements take place at as low an 
orbit as possible. Because of the mathematical characteristics of the 
potential field, the resolution of any gravity measurement at altitude is 
roughly equivalent to the altitude. A low orbit, however, has a very 
short lifetime because of atmospheric drag, and a short lifetime makes 
it difficult to obtain the complete coverage of the earth that is also 
desired. 

Therefore, a low orbit with orbital rarameters is needed such 
that the orbital tracks interleave so complete coverage is obtained in 
a period shorter than the orbital lifetime, and where the track spacing 
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Fig. 3 . Unmodelled Anomalous Gravity and Gravity 

Gradient at Satellite Altitude (300 cm) (from 
Ref. 3). 
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Fig. 5. Regions of Best Sensitivity for Altimeter, 
Doppler, and Gradiometer by Harmonic Degree 
and Altitude (from Ref. 5). 
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is matched to the swath width (equal to the altitude). A set of orbital 

parameters is existent that fits these requirements fairly well. At an 

6 7 

orbital altitude of 270 km, an ’’integer orbit” * exists. The orbital 
track repeats upon itself after exactly 16 orbits. This orbit can be polar, 
with 16 orbits per sidereal day, or a sun synchronous orbit (at a slightly 
different altitude and inclination) with 16 orbits per solar day. If the 
altitude is slightly higher or lower, then the orbital track drifts so that 
the 16th orbit is displaced to either one side or the other of the first 
track. These offset orbits finally begin to repeat after a number of 
days when the drift has caused the satellite track to overlap the second 
ground track. Two of these orbits are of interest. They repeat after 
about 5 days, and their track spacing is approximately equal to the 
altitude. One is a polar orbit at about 320 km that repeats after 
79 orbits, and the other is a polar orMt with altitude of 220 km that 
repeats after 81 orbits. The track spacing between the half arcs for 
both orbits is approximately 250 km, so that there is a good match 
between the track spacing and the swath width. 

In reality, the orbital altitudes decay as a result of drag, so 
that these simple orbital path models are not followed exactly, HRL 
presently envisions launching into a 3 30 -km polar orbit and allowing 
the altitude to decay through these two altitudes where overlapping 
coverage is obtained. Hughes has chosen a polar orbit rather than a 
sun synchronous orbit in order to obtain full coverage of the earth and 
provide for calibration points twice per orbit at the two poles. The 
orbital lifetime estimated for the mission is approximately 30 to 50 days. 
The time spent near 320 km would be long enough to obtain good cover- 
age of the earth at that resolution (640-km wavelength or degree 62). 

As the altitude decreases, resolution will improve. There should be a 
substantial amount of coverage at around 220-km altitude with excel- 
lent resolution (440-km wavelength or degree 90), but some coverage 
will be lost because of the rapidly decreasing altitude and the fact that 
the track spacing at the equator of 250 km slightly larger than the 
sensor resolution. 
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Non-Eclipse Orbits (See Appendix E) 


It would be desirable to launch the gradiometer satellites into a 
polar orbit of the earth that does not cause the satellites to be eclipsed 
by the earth throughout the mission. The advantages of the non-eclipse 
orbit are the weight reduction and reliability increase available by 
elimination of batteries for electrical power during the eclipse portion 
of the flight. Also, the thermal control system required for the sensor 
would only have to contend with one state of thermal equilibrium rather 
than cycling between two. 

HRL has investigated possible non-eclipse orbits and has found 
that even despite the relatively low orbits under considerations, it is 
possible to achieve non-eclipse periods several times longer than the 
estimated lifetimes for these orbits. To attain these orbits only requires 
that a launch window constraint be placed on the mission. The satellite 
is launched 1 5 to 20 days before either the vernal or autumnal equinox 
(21 March or 21 September) into a polar orbit chosen such that on the 
day of the equinox, the orbital plane coincides with the terminator plane. 
At this point in time, the ecliptic and celestial poles of the earth are all 
in the terminator plane. The slow rotation of the terminator plane 
about the ecliptic poles causes a drift between the terminator plane and 
the orbital plane (~l°/day). However, simple calculations show that 
with this choice of launch time and orientation, it is possible to have 
non-eclipse periods in excess of 30 days even for orbital altitudes 
below 2 50 1cm. 

Non-eclipse orbits could also be chosen using sun-synchronous 
orbits lying near the terminator plane, the non-eclipse period is then 
theoretically infinite. The foregoing shows, however, that with this 
minor constraint on launch time, the advantages of a non- eclipse orbit 
may be achieved while the self-calibration and full coverage aspects 
of the polar orbit are retained. 
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c. 


SPACECRAFT PARAMETERS 


Because of the relatively low field strengths estimated for the 
higher order harmonics of the earth’s field, the sensor must oe made 
as large as possible. The arm length of the sensor is primarily 
determined by the maximum radius obtainable in the spacecraft, which, 
in turn, is determined by the payload envelope of the launch vehicle. 

If a Scout launch vehicle is used to keep costs down and reliability up, 
a number of launch shroud configurations that have been developed for 
this vehicle may be used. One of the largest in diameter is the 42 -in. 
diameter shroud mentioned in the Scout users 1 handbook. The allowable 
payload diameter for this shroud is 96. 5 cm (36 in. ). The cylindrical 
portion of the payload envelope with this diameter is 84 cm (33 in. ) 
long, which allows space for two cylindrical spacecraft 96 cm in. 
diameter by 42 cm thick. A very preliminary spacecraft design of this 
size is shown in Figs. 6 and 7. The sensor arm length with this space- 
craft configuration is about 40 cm. 

The front part of the payload envelope can be used for a spin-up 
and attitude control system that inserts the two spacecraft into orbit 
with the proper attitude and spin speed. After the payload attains 
orbit, the spin control mechanism increases the satellite spin speed to 
the desired rate (about 240 rpm) and orients the spin along the orbital 
track. After release of one spacecraft, the jets are used to torque 
the other spacecraft so that its spin vector is perpendicular to the 
first spacecraft. With the two spacecraft in this relative orientation, 
one craft measures the vertical gravity gradient and the cross -track 
horizontal gradient, while the other measures the along -track horizontal 
gradient and a redundant measurement of the cross-track gradient. 

After 1/4 of an orbit, the orientation of the two satellite spin axes 
relative to the orbital track changes, and the data output from the two 
sensors is interchanged. Although this is a relatively complex mode 
of data collection, it does allow for the measurement of more components 
of the gravity gradient at the same time. Most important, this mode 
of operation allows the cross -track gradient information to be obtained, 
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which can be used to tie the data together across the orbital tracks. The- 
cross-track gradient along one track can be used to predict the gravity 
field at the next track. This closure property of the data sets can be 
used to eliminate drift errors. 

A simpler version of the experiment is to launch a single satel- 
lite and torque the spacecraft spin axis so that it lies in the plane of the 
orbit. The advantages of this mode of operation are that the spin axis 
of the spacecraft does not change orientation with respect to the orbit, 
and the drag torques remain constant. In this orientation, the 
gradiometer measures the difference between the vertical gradient 
and the along -track horizontal gradient, and their orientation with 
respect to the local vertical. 

A major interaction of the sensor and the satellite is the 
dynamical interaction of the two mechanical structures. In its work 
prior to the contract, HRL performed a dynamic analysis of the inter- 
action of the sensor and the spacecraft (see Attachment B). For the 
design phase of the contract, HRL improved this analysis for the 
specific mission and sensor design. A report of this work constitutes 
Appendix C. 

As a result of this dynamic analysis, HRL determined a number 
of spacecraft parameters that should be controlled. The spacecraft 
transverse moments of inertia should be the same to 1%, and the align- 
ment of the spacecraft spin axis and the direction of the sensor torsional 

_3 

axis should be aligned to within 10 radians (the position of the center 

of mass of the sensor is not critical), and the spacecraft coning angle 

_ 5 

should be kept below 10 radians with nutation dampers. 

The results of the analysis are summarized in Table I, which 
shows the significant gradiometer errors caused by sensor-satellite 
dynamics. The table shows the sources of the error signals, their 
frequency, and their equivalent amplitude in EU. 

The dc torques between the sensor arms do not generate any 
steady state signals because of the finite resistance across the 
transducer. Those terms at or near the spin frequency will be cut by 
a factor of greater than 100 by the integrating circuits of the electronics. 
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TABLE V 


Significant Gradiometer Errors Due to Sensor Satellite 

Dynamics 


Error Ouputs 

Frequency 

Equivalent 

Amplitude 

(EU) 

Rotational Field 

dc 

40 

Du. ing Function 



(W .0) . ) 

i D 

2/k.k 0 ui 
12 s 

10 2 

Sum Mode Mismatch 

dc 

0.4 


/k.k n a) 
12 s 

10 _1 

Mass Unbalance 

dc 

5 x 10 5 

(r ) 

Umax 

a) 

5 

5 


dc 

20 


T368 


where 


= spin frequency 

/ kjk^ 1 = 0. 99 

Thi3 leaves only one term, the second term of the rotational 
field driving function, which is on the threshold noise level and close 
to the sensor detection frequency. If this term is traced back to its 
origins, it is found to be strongly dependent on the spacecraft coning 
angle. Any increase in this angle above the assumed 10”^ radians 
will cause signals above 10 2 CU. This signal varies as the square 
of this coning angle. 

Since the HRL design study effort, a complete and detailed 
Earth Physics Satellite design has been developed by JPL* which includes 
an extensive discussion of the desired satellite parameters and the 
methods of achieving them. The reader is referred to this study report 
for this more up-to-date design. 
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desired satellite parameters and the methods of achieving them. The 
reader is referred to this study report for this more up-to-date design. 

D. SENSOR PARAMETERS 

Most sensor parameters are determined by mission and space- 
craft constraints. The desirability of obtaining 0.01 EU sensitivity 
indicates the need for a sensor arm length that is as long as possible. 

A sensor arm length 76 cm from center to center of the end masses 
(86 cm overall) was selected as the largest arm diameter possible 
for the 96-cm spacecraft diameter, which, in turn, is dictated by the 
Scout payload envelope of 106. 5 -cm diameter. The chosen arm end 
masses were 2 kg each, this weight considered as being reasonable fer 
the size of the sensor. 

The 3 5-sec sensor time constant was derived by using the time 
required for the spacecraft to pass through one resolution element at 
the nominal altitude of 270 km at the orbital velocity of 7. 75 km/sec. 
This figure was considered a reasonable optimum between the 41 sec 
for 320-km altitude and the 29 sec for 220-km altitude. With this size, 
weight, and time constant for the sensor, the thermal noise caused by 
the Brownian motion of the sensor structure has an equivalent noise 
level of 0.007 EU . This ultimate lower limit was distressingly close 
to the design goal of 0.01 EU, but could not be lowered within the 
constraints imposed by the Scout payload envelope. This sensor system 
time constant is the smoothing time to be used in the sensor data 
preprocessing. The sensor output should be sampled approximately 
once every 1 to 5 sec to overcome digitalization noise, prevent aliasing, 
and pick up strong, short period signals resulting from dense localized 
anomalies . 

The sensor frequency of operation is not critical and is set by 
conflicting requirements. This frequency should be as low as possible 
to ease the spin speed stress requirements on the satellite structure, 
and should be high as possible to avoid the low-frequency noise in the 
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electronics and for ease in laboratory testing, where it is difficult to 
obtain adequate vibrational and acoustic isolation for mechanical 
structures below 10 Hz. The selected design frequency is 8 Hz, which 
implies a spin speed of 240 rpm (4 rps) for the satellite; although fast, 
this speed is not unreasonable. As is mertioned in Section III, the 
measured sensor frequencies are about 5 Hz. 

With the 8-Hz sensor frequency and 3 5-sec sensor time constant, 
r the desired sensor quality factor, Q, is 

Q = tt f t = 800 

Previous work indicated there should be little problem in obtaining a 
Q of this value with a structure such as this, since Q's from 
300 to 2000 had been formerly obtained. Section III shows that HRL. 
did not achieve this Q level in the sensor^ although a redesign should 
provide an increased Q. 
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SEC TION III 


GRAD IOMETER DESIGN, FABRICATION, AND TEST 


A. DESIGN CONCEPT 


During the contract study phase, several problem areas were 

anticipated which had to be accounted for in the sensor design. Primary 

among these was the fact that to obtain usable data from a 270-km 

altitude satellite orbit, the gradiometer must have a basic resolution 

3 GM 

level of 0.01 EU ( y ^3 )* This implies resolution of two parts per 
million in a background gradient field of 4500 EU . Therefore, several 
design areas of no great concern in previous sensor designs were 
considered significant for this application. These included: 

1 . Variations of resonant frequency and Q with 
temperature 

2. Reduced sensitivity to magnetic gradients 

3. Improved signal level 

4. Reduced electronic roise 

5. Reduced thermal noise 

6. Reduced inertial loading sensitivity 

The sensor design was therefore created with the foregoing items as 
primary bounding conditions. 

It is recognized, of course, that long-term variations in the 

signal can be eliminated because the sensor is recalibrated each time 

it passes over the north or south pole. Only those signal errors in the 

same frequency spectrum as the data to be obtained are of concern, 
th th 

e. g. , from the 10 to 80 spherical harmonic of the earth's field. 

The design drawings are attached to this final report as an 
appendix (Appendix A). The breadboard sensor, as it was assembled, 
is shown in Fig. 8. Although there are no magnetic shielding covers 
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shown in the photograph, covers were designed that provide adequate 
magnetic shielding for the sensor (see Appendix F). This cover 
design is shown in drawings xl 146, xl 147, and xl 149- The parts were 
not manufactured because of the high cost and because magnetic testing 
of the instrument was not planned under this contract. However, the 
design should provide the required magnetic shielding when the covers 
are lined with high permeability foil. 

The transducer was designed to provide a very high voltage level 
by using large-size, piezoelectric, bender transducers. Most of the 
sensor damping was found to be in the piezoelectric of the transducer; 
this had an adverse effect on sensor Q by keeping it below 125. To obtain 
the desired 3 5-sec total time constant it is desirable to have a sensor Q 
as high as 62 5. Such a Q lowers the thermal noise from 0.015 EU 
to 0.007 EU. 

Jt is suggested that for future models of this sensor, the amount 
of piezoelectric material be reduced. Such a reduction raises the Q 
and lowers the sensor scale factor, a trade off that leads to a more 
optimum unit and is one of the recommendations made in Section V. 

The signal preamplifier was found to have an equivalent noise 
level of 0.002 EU equivalent, which was well inside the requirement 
(see Section III-A). 

The size of the sensor, and to a large extent its weight, was 
determined by the thermal noise requirement. The breadboard was 
made deliberately heavy ( = 70 lb) to allow experimental testing in a 
1 g environment. For actual satellite operation, the weight could be 
probably cut to 1/3 of the present weight. 

The arm design and its support structure were determined by 
requirements to minimize the inertial loading effects on the transducer 
(see Appendix D). 

Materials were chosen to reduce sensor magnetic sensitivity 
and to minimize the termperature sensitivity effects (see Appendix G). 
Despite these efforts, HRL expected a considerable shift of xesonant 
frequency with temperature and decided that additional control could 


21 



be provided by using one of the transducers as an adjustable spring 
damper system by padding it with resistance or capacitance. Results 
of these tests are contained in Sections III-G and III-H. 

B. ASSEMBLY AND BALANCING 

Once the manufactured parts were in house, assembly proceeded 
with very little problem. Some difficulty was experienced in assembly 
of the transducer into the support brackets and better fixturing would 
have simplified that effort. 

Balancing of the arms of the first assembled sensor proved to 
be a difficult task. Data were inconsistent and varied with loading 
direction. The balance did not remain constant from day to day and 
balancing of the arms to within 0. 003 in. of the center of support was 
unsuccessful. 

Calculation of the stiffness of the transducer supporting 
structure (Fig. 9) disclosed that this structure was only a factor of 
2 stiffer than the transducer itself. Deflections in the structure were 
suspected as being the cause of spurious signals from the transducer. 
Following redesign of the structure (Fig. 10), little difficulty in balanc- 
ing the final sensor assembly to 0100014 in. (1 gram on the end of the 
arm). 

C. ELECTRONICS 

1. Amplifier Design 

The signal amplifier is a low-noise instrumentation amplifier 
consisting of an FET input section, a main amplifier section, and a 
signal adder section. The electrical schematic is shown in Fig. 11. 

The parts list follows. 

The FET input section consists of a matched pair of low-noise 
FET' 8 with shunting diodes on the gates for overload protection against 
large signals. 
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The main amplifier section consists of two low-drift Op-amps 
with negative feedback coupled to the source leads of the input FET’s 
for gain stabilization. 

The last section consists of a low-drift Op-amp operated as a 
subtractor for differential signal addition and common mode rejection, 

2. Amplifier Tests 

Gain, noise, temperature, and voltage tests were run on the 
electronics package. The nominal amplifier gain is 1815 over the 
bandwidth from 1. 5 to 500 Hz. The common mode rejection is better 
than 90 dB from 4 to 500 Hz with a value of 94 dB at 8 Hz. The shorted 
electronic noise level was measured and shows the expected variation 
with frequency (see Fig* 12). The noise is flat at 14 nV/Hz^^ from 
20 to 500 Hz and shows the expected 1/f behavior below 20 Hz. The 
noise level is 3 5 nV at the operating frequency of 5, 78 Hz with a 1-Hz 
bandwidth or an input noise characteristic of 3 5 nVA^Hz. If a total 
time constant of = 3 5 sec is assumed, the effective bandwidth 
would be 


B = ^ = 1. 4 x 10‘ 2 Hz 

Zr 

and the input noise with a 35-sec integration time would be 4. 3 nV. 
Using the scale factor of S = 2 p.V/EU calculated in Section I1I-D, the 
electronic noi9e equivalent gradient is 

(r..\ = - r ~- ~ 0. 002 L 

V ij/amp \JB S 

An additional long integration time test was run on the amplifier 
with shorted input. In this case, the output signal was run directly 
into the data processing electronics, and a noise trace was run, using 
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Fig. 12. Noise Output Versus Frequency. 
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an integration time of 3 1 sec. This noise curve is shown in Fig. 13. 
Here the output of the amplifier was = 0. 075 mV (rms). Using a 
total gain of G - 18, 150 and a transducer scale factor of 2 pv/EU 
(see Section III-D), the equivalent noise input level can be calculated 
at 


V 

= = 0. 002 EU. 

amp Ga 

D. SENSOR SCALE FACTOR 

1 . Expected Sensor Scale Factor 

The expected sensor scale factor can be calculated by evaluating 
the expected moment on the transducer for a 1 EU gradient input and 
its resultant voltage output by formulas from the piezoelectric 
catalogues. 

The torque level for a 1-EU input gradient is 

_3 

T = r)l r.. _ TT . = 6. 08 x 10 dyn-cm 
zz (1 EU) 



where 


n is the arm efficiency factor = 97.4% 

I is the arm polar moment of inertia = 6.24 x 10^ gm-cm^ 

- 9 -2 

r is the input gradient = 10 sec 

This torque appears as a force at the end of the* transducer where 
1/2 FL = T. The transducer scale factor is given in the piezoelectric 
handbook as 


V 3_ 1 

FL " 4 *WT 
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where 


-3 

j = 11.4 x 10 for the material (channelite 5500) 

W = 1 in. = 0. 0254 m 

T = 0. 232 in, = 0. 00589 m 
L = 3 in. = 0. 0762 m. 

Therefore, 


X. = 57. 15 V/n-m. 
FL 


Now the force at the end of the transducer is 


F = jtIj- -- 4 x 10’ 4 dyn 


and FL = 3. 04 x 10 ^ dyn cm = 3. 04 x 10 ^ n-m; consequently, the 
voltage output for a 1-EU input V = v/FL (FL) = 1. 73 x 10 8 V/EU. 
This assumes nonresonant operation. With a Q of 12 5, V/T would be 
(1.73 x 10' 8 ) (125) = 2. 17 x 10 -6 V/EU. 

2. Sensor Scale Factor Test 


a. Test Setup 

The sensor was calibrated for scale factor by adding a 
known mass to the end of one arm thereby creating a deliberate 
center-of-mass unbalance in that arm. The sensor was then subjected 
to linear vibration perpendicular to its spin axis and also perpendicular 
to the mass unbalance vector. The input acceleration level was 
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measured by calibrated geophones and the sensor output was noted 
through the data readout electronics. 


Geophone Output 
Geophone Scale Factor 

Unbalance Weight U 

Sensor Resonant Frequency f 
Sensor Output V q 


12. 5 mV (O-p) 

56 V/ in. / sec 
(<S>5. 78 Hz) 

81 g at 1 5 in. = 3086. 1 gm-cm 
5. 78 Hz 

9. 5 mV (gauge 1) (O-p) 

11. 25 mV (gauge 2) (O-p) 

10. 0 mV (O-p) avg. 


b. Input Acceleration Calculation 

The input peak acceleration car be calculated from the 
geophone output and scale factor 


peak 


V 

= _L 
Sg 


(2 jrf) = 8. 1 x 10 


■3. . 2 

in. /sec 


= 2. 06 x 10 


-2 . 2 
cm/sec . 


c. Torque Scale Factor Calculations 

The torque on arm 1 is given by 


T 


1 peak 


a peak 


63. 5 dyn-cm 


while the torque on arm 2 is zero 
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Th(* rct,o of sensor output voltage to peak input difference torque is 
thereto re 


V 

o 


T 


1 


1. 57 5 x 10 


V/dyn-cm. 


d. Gradient Scale Factor Calculations 

For a point mass gravitational anomaly (mass M, at 
distance R), the difference in torque on the two arms of a rotating 
gravitational gradiometer is given by 


Tj - = 2r) C( ) sin 2 u>t 

where C is a single arm polar inertia (6. 24 x 10^ gm-cm^ for this 
sensor) 

q - arm efficiency factor (0, 97 for this sensor) 

3 GM -9 -2 

T.. = y — r an equivalent gradient. = 10 sec for 1 EU 

lJ C R 

(This equation has been developed in Hughes literature many times 
before. ) 

The peak torque -gradient conversion for a 1- EU signal (I\j) 
can now be determined: 

(T 1 • T 2 ) 1 EU (peak) = L21 x 10 “ 2 <*yn-cm. 
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The scale factor can now be calculated, combining the torque 
and gradient scale factors 


S 

s 




peak 1 EU 



1. 91 x 10 


' 6 V/ EU 


or approximately 2 pV/EU. 

This is very close to the theoretical value calculated in 
Section III-D- 1. 

E. SENSOR NOISE LEVEL 

1. Thermal Noise 


The sensor thermal noise is calculated, using the equation 


2 

n 


2k T 

c Vt 


where t is the sensor time constant (Q/?rf) and T f is the total 
(filtered) time constant (3 5-sec max). 

The final sensor Q was 126 and f was 5.78 Hz; therefore 


and 


Tj = 7 sec, 





- 14 

2 8x10 ** 

0.97 6.24 x 10b(7) (35) 


{ 


0.015 EU 


I 
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It should be pointed out that <T g is defined as the gradient 
equivalent thermal noise standard deviation, based on the following 
torque equation for a gradiometer rotating in the x-y plane: 


T, 


U (r . 


yy 


F ) + a 1 sin 2u>t + (Z T t <x ) cos 2u;t 
xx s J xy c 


where r\ and C are as previously defined. 

r , r f and r are components of the general gravitational 

gradient tensor, and IF is equal to 7F (but statistically independent). 

( The total sensor output thermal noise IF + 7F - Jlv 0. 0212 EU ) 

The preceding definition implies that if operation on the output 

signal is conducted in such a way as to obtain T- . (as previously defined, 

3 l J 

3GM/2R ), the thermal noise level on this signal will be <r /2, or 

s 

0.007 EU. The following illustrative example is given in proof of this 
statement. 


Example: 

A spherical anomaly of 10^ kg whose center is at 300 km 
distance from Ihe gradiometer along the y axis has the gradient 
components 


r = r 

XX zz 


GM 

R* 


= - 0. 02 5 EU 


and 


yy 


= + 


2GM 


= G. 05 EU 


(r *y ■ ">• 
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Incorporating these values into the torque equation, a peak amplitude 
of 


T. - T = qC (^~ - ±yfz <T ) = 4. 54 X 10‘ 4 ±1. 28 X 10' 4 dyn-cm 
! 2 R * s 


is obtained. If this torque difference is multiplied by the voltage-to- 

-4 

torque ratio developed earlier (1. 575 x 10 V/dyn-cm), an output 
voltage of 


V = 71. 5 nV ±20. 2 nV. 
o 

is obtained. If this signal is phase split by a phase sensitive demodula- 
tion scheme, using the x and y axes as references, the sine and 
cosine components of this output voltage are 

V = 71. 5 nV ±14. 3 nV 
os 


and 


V = ±14. 3 nV. 
oc 

If V alone is now evaluated for the equivalent gradient and equivalent 
noise (I\.) by dividing mean voltage, and standard deviation by the sen- 
sor scale factor (S = 1.91 pV/EU), 

s 

r.. = 0 . 0375 EU ±0. 0075 EU 
ij 


3 

is derived where the mean gradient is equal to 3GM/2R , and the noise 

gradient is equal to <r / 2 as expected. 

8 
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z. 


Isolated Noise 


To check its noise level, the sensor was mounted on a 
pneumatically supported vibration isolation table. Data were taken for 
a 15-min run, using a 31 -sec integration time. The noise data are 
shown in Fig. 14. The RMS of these data is 8 mV, which is equivalent 
to a transducer output of 0. 008/ 1815 = 4. 4x10 ^ V or 2,2 EU 
(equivalent). 

It should be noted that this sensor was only moderately isolated 
from floor noise (resonant frequency of the suspension being in the 1- to 
2 -Hz region), and there was no acoustic isolation provided. The sensor 
output would increase by several orders of magni'ade with the sounds of 
door slams or heavy footfalls. It would be anticipated that at least a 
2 order of magnitude vibration level reduction in a spacecraft environ- 
ment would be expected. (Vibration level on the table was approximately 
10' 5 g RMS. ) 

F. RESONANT FREQUENCY AND Q DETERMINATION 

The first test performed on the sensor after assembly was that 
of running a frequency sweep by using ihs second transducer r.s ? driver 
to excite the sensor arms in the gradient detection mode and sweep 
the excitation frequency to obtain a resonant frequency curve. This 
test was performed on both sensor assemblies with the two different 
transducer support structures. 

The resultant Q curves are shown in Figs. 15 and 16. The 
first curve is with the old, weaker, transducer support structure and 
has a natural frequency of 4. 59 Hz with a Q of 138. The second curve 
(final version of the sensor) has a resonant frequency of 5. 65 Hz and 
a Q of 126. The resonant frequency difference is directly attributable 
to the modification of the transducer support structure. Although the 
Q in both these tests was greater than 12 5, a repetition of this value 
of Q was never achieved in later tests. Throughout the temperature 
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Fig. 14. Amplifier Noise (81-gram Unbalance). 
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Fig. 15. Sensor Q Curve (Old Transducer). 
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Fi g. 16 . 


Sensor Q Curve (New Transducer). 
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tests, the Q was between 70 and 85 at room temperature. The theory 
is that after the sensor was physically handled and moved, the bonding 
between the transducer halves was fractured in places, and certain 
cement joints and bond lines possibly weakened. The sensor stops 
were set to ±3°, but the sensor was roughly handled during moving, 
and this 3° limit may have been exceeded by actual bending of the stop 
support structure. 

G. TEMPERATURE TESTS 

Once the resonant frequency of the sensor was determined, 
detailed temperature tests were conducted on the gradiometer. The 
capacitive drive was used to excite the resonant mode, and resonant 
frequency curves were run with the sensor temperature stabilized at 
23. 5°C, 30°C, 3 5°C, 40°C, and 50°C. These curves were run on 
both models of the sensor (old and new transducer support structures). 
The data are shown in Figs. 17 and 18. Results of both tests are 
nearly the same. The frequency shift rate is C. 043 Hz/°C for the 
older sensor and 0. 038 Hz/°C for the newer one. This shift rate 
difference amounts to a percentage of »0. 66 %/°C. Q shift for the older 
sensor was 3. 5/°C at the highest part of the temperature range and for 
the newer, 1. 24/°C. These shifts are easily compensable by the trans- 
ducer padding methods discussed in Section III-H, 

H. FREQUENCY TUNING 

The sensor internal electromechanical parameters were calcu- 
lated using a capacitive load across one of the transducers, then the 
response of the sensor to various resistive and capacitive loads was 
determined; finally, the calculations were checked by using a 1 megohm 
resistor across the transducer and comparing the actual Q and f with 
those calculated values. 
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Fig. 17. Sensor Resonant Frequency 
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With no loading the old sensor parameters were 


4. 816 Hz, Q -- 60 . 


Adding a 0. 05 mf capacitor across one transducer the parameters 
changed to 


f n = 4.783 Hz, Q 
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The equivalent electromechanical circuit can be modeled as shown in 
Fig. 19. 



Fig. 19. Equivalent Electromechanical Circuit 
Model . 


which when simplified into a series circuit becomes 



tors -is 


'a 


Fig. 20. 

Model Simplified into 
Series Circuit. 
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where C varies, depending on the value C^. (When is zero, 

C = 6.461 x 10“ 9 F; when C T = 0.0 5 mF, C = 1 1. 016 x 10‘ 9 F. ) 

The general equation 


2 

co 



in the two cases (with and without C^) 


can now be used where C is the total capacitance of the circuit to solve 
for C c and 


C s = 2. 207 X 10' 1 ° F 


= 5.133 


10 6 H 


Using tne equation 

R = (R ex = 653, 290 £2 for R L = 10 M) 

Solving the equation gives 


R s = 1. 937 x 10 6 a 

Curves can now be constructed, showing the effect of any tuning 
capacitor on the natural frequency, and any resistor on Q. For these 
curves, the new sensor natural frequency (untuned) of 5. 75 Hz was 
assumed. (The new C g was therefore calculated as 150 x 10*^ F. ) 
These curves are shown in Figs. 21 and 22. 

I. COMPENSATION REQUIREMENTS 

If the slopes of the two tuning curves are now calculated at what 

- 8 

appears to be convenient generating points (C^. = 10” F, R^, = 100 K ft), 
it is found that 
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Fig. 21. Resistance Versus Sensor Q. 
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and 



-0.85 x 


10' 6 Hz/pF 


dQ 

dR T 


-2. i x io _5 /n 


The requirement of maintaining an accuracy of 0. 01 EU in the 
presence of a 4500-EU signal (3GM/2R^) has already been noted as 
implying an overall amplitude measurement accuracy of 2 ppm. This 
requirement can be used to establish th resolution required on the 
tuning components of the tuning circuit. Amplitude is a direct function 
of Q; therefore, a 2-ppm variation in Q represents a 2-ppm variation 
in amplitude. This 2-ppm requirement on Q, therefore, necessitates 
a resolution in the tuning resistor of 


AR 


2 x t0~ 6 x Q ^ 
2. 1 x 10" 5 


AR 


5. 75 U out of 1 0 il 


which should be attainable. 

The amplitude variation with natural frequency is somewhat 
more complex because this variation is also a function of Q. 

The steady state amplitude near resonance is given by 


X = 




7 
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at resonance 


X 


R 


P 



Therefore, 



l + 


*4 



2 


or in terms of deviations from the resonant condition, 


1 




J 
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if AX/Xr must be held to 2 ppm then 

Aw „-i /AX . . ia -3»-1 

“T Q \/x- = '- 4 * 10 Q 

n > rv. 

or in tabular form: 


Q 

Ago 

CO 

n 


. 5 

50 

2. 8 x 10 


- 5 

100 

1.4 x 10 


-6 

200 

7.1 x 10 

400 

3. 5 x 10' 6 

600 

2.4 x 10" 6 


With a high Q system (Q = 600) Aco/go^ must be held to ^2 ppm which 
implies a frequency resolution of 

Af = 2 x 10‘ 6 w = 1 .6 x 10‘ 5 Hz . 
n 

The slop:? of the frequency versus capacitive curve was 0.85 x 10 ^ Hz/ 
pF which implies a capacitance resolution of AC ~ 19 pF . Again, one 
which appears feasible. 
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SECTION IV 

CONCLUSIONS 


In this report, a set of mission and spacecraft parameters has 
been established that appear to be within the state of the art. Addition- 
lly, supporting studies by JPL indicate the attractiveness of such a 
mission as an alternate and/or verification test for the radar altimeter 
earth physics satellite. It has also been shown that such a gradiometer 
mission is feasible. The gradiometer design parameters are such that 
with some modification a resolution of 0.01 EU at 35-sec integration 
time could be reliably obtains 1 from such an instrument operating in a 
spinning spacecraft. 

Certain aspects of the sensor design need further study. Pri- 
mary among these is the desirability of reducing sensor damping in 
order to improve the sensor Q and thereby obtain an improvement in 
sensor thermal noise threshold. The adjustment of Q and resonant 
frequency by padding a piezoelectric spring was demonstrated, but an 
effective system to provide the adjustment of the tuning components has 
yet to be designed. 
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SECTION V 


RECOMMENDATIONS 


As a result of its experience with the breadboard model of the 
rotating gravity gradiometer, HRL has seven prim pal recommendations 
for future work on the sensor. 

1. The amount of piezoelectric in the transducer structure 
should be reduced by half to two-thirds and replaced by aluminum to 
raise the natural irequency to 8 Hz or more, and probably increase the 
Q significantly. This ratio should be adjusted to obtain best «^gnal to 
noise, as when the piezoelectric material is reduced, the voltage 
output per EU drops. 

2. The bellows method of decoupling the transducer struc- 
ture from the arms (except for the desired torsional coupling) failed. 
Because bellows have a number of undesirable cros s -coupling modes, 
a better design solution is needed. 

3. The electronics for frequency and Q control of the 
censor (See Section III-H) should be constructed and tested. 

4. The vertical separation of the center of masses of the 
two arms caused some cross-coupling problems. A design concept 
tha^ brings the centers of mass of the arms closer together, but 
retains the biaxial torsional suspension for each arm, was recently 
developed by HRL under Air Force Contract F19628-72-C-0222, and 
should be considered for the next NASA design. 

5. The arm length and end masses of the sensor should be 
increased as much as possible to lower the sensor thermal noise level 
well below the desired 0. 01 EU threshold sensitivity level. This 
would imply the use of a different launch vehicle than the Scout. 
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6 . 


The proposed magnetic shielding should be fabricated 
and its effectiveness demonstrated in the laboratory. 

7. A very low frequency (0. 1 Hz) vacuum enclosed 

isolation suspension should be designed and fabricated so that the 
thermal and electronic noise level of the sensor can be demonstrated 
on the ground. 

This test could be combined with a gradient detection test 
using rotating masses in the vicinity of the sensor to calibrate the 
sensor by gravitational torque excitation rather than the mass unbalance 
torque excitation method used in the completed program. 

There are a number of recommended areas of study for the 
spacecraft, such as the development of adequate methods for measure- 
ment and control of temperature, spin speed, attitude, spacecraft 
dynamics, and data digitalization and reduction. These are adequately 
discussed in the JPL study. 
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APPENDIX A 

ASSEMBLY INSTRUCTIONS 

L bUB ASSEMBLIES 

1 . Transducer Assembly 

Preparation of Transducers (L-B-936308) 

Using an cxactc knife or equivalent, remove 3/8 x i. 0 in. 
of plating at each end of the transducer (mounting areas). Sand clean 
with 320 grade silicon carbide paper. Silver is to be removed from 
both sides at both ends. No conductivity is to remain on transducers 
at mounting areas. Cement 0. 015 x 0. 2 50 in. phenolic strips to each 
side at ends of transducer. Cement 0. 015 x 0. 060 in. phenolic strips 
across edge of transducer to prevent shorting to mount. Accurately 
measure slot of mount and cement brass shimstock (selected X. 12 5) to 
sides of previously mounted phenolic to provide not more than 0.001 in. 
or less than 0. 0005 in. side clearance. This ensures perpendicularity 
and parallelism of transducer to support arms. Note 60°C max for 
cure. 

b. Preparation of Plate Support (L-C-936315) 

Remove balance weights to allow accessibility while 
cementing transducers in position. Cement 0. 020 enameled wire in 
0. 062 holes, leaving approximately 1 in. exposed each end. Remove 
enamel from each end, but no closer than 0. 050 of part body, to ensure 
infinite resistance. Clip outboard ends to 1/4 in. after tinning. 

Inboard ends are to be clipped to convenience when attaching to gauge 
after assembly of gauges to support arms. Balance weights are to 
remain off until completion of transducer assembly. 
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c. 


Assembly of Transducers to Plate (L-C-936314) 


Clean parts in acetone. Assemble support plate, using 
bolts and spacers. Check parallelism to 0. 001 in. at edge of plate; 
dress spacers until this dimension is achieved. Check fit of transducers 
with shims in their respective slots. They must fit smoothly and with- 
out distortion. Orient transducers to slots so that center terminal 
strip is located in short post slot adjacent to cemented wire terminal. 
Using wooden shims and wedges, block transducers to position. Do not 
wedge tightly . Transducers should be sufficiently free to allow removal 
with shims in place. For convenience, one gauge may be cemented in 
place per operation. After fit is established and transducer end is 
coated with epoxy, it may be inserted in place and cured for two hours 
at 50° to 60°C. When both gauges are initially set, a tinkers dam 
(coated with release agent) is inserted in slct to hold additional epoxy. 
Slot is to be filled flush with o. d. of support. Cure temperature is feO°C . 
Gravitational forces require a separate curing operation at each end. 
Central terminal strip of gauges is soldered to 0. 020 in. wire and 
flushed before final epoxy operation is performed. Connect both sides 
of transducers in parallel to adjacent terminal with 0.005-in. diameter 
enameled wire. Install counter weights, using Locktite-G on threads. 
Tighten securely and leave in secure storage until ready to complete 
assembly. 

Preparation of Adapter Arm (L-D-936323) 

Install posts with large dia diametrically to the C/L of 
two smaller dia posts. Use Locktite-G on threads. On large dia post, 
slots must match holes in plate. Cement 2 pieces of 0.020-in. dia 
wire in place. Allow l/4-in. max overhang at the plate end. Allow 
min of l/2 in. at post end for attachment of flexural leads. 

e * Transducer — Final Assembly of Supporting Elements 

Add supporting plate (L-D-936323) to transducer plate 
(L-C-9363 1 5), using spacer. Orientation of support plate s with C/L 
larger dia post aligned with C/L of shorter transducer support. 

Check concentricity on level block of each of the four plates at four 
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points on each circumference to max of 0. 002 TIR misalignment. 
Further errors of parallelism will be corrected during attachment of 
sensor arms. 

2. Sensor Arm Assembly 

Insert arm weights with arbor press. If hole diameter is on 
low end of tolerance and diameter of weight is on high end, the arm can 
be heated to max of 78°C for relief. Caution : Apply heat by oven only. 
Uneven or pinpoint heating will distort arm. Place arm in secure 
storage until needed for final assembly. 

3 . Installation of Transducer Assembly in Housing 

Insert assembly in housing. Orientation is established by holes 
in each end of housing. Use small parallel clamp to secure transducer 
in housing and prevent damage by impact. Place in secure storage 
until needed for final assembly. 

Note : Use acetone for final deg* casing and cleaning of all 
parts and assemblies. Connect transducer leads to plate terminals. 

II. FINAL SENSOR ASSEMBLY 

Mount support plate on leg supports and orient with base 
terminals that are outside evacuation area seals, facing up. Mount 
transducer housing assembly in such a manner that ends are equidistant 
from C/L of plate. Install torsion bar to housing. Note : Make certain 
that dowel is slide fit; otherwise, disassembly can. become difficult. 

Using 4. 5-in. blocks, place on support plate in such a manner 
that the ends of arms are supported. The sensor arm may be lowered 
over torsion bar. Note : Do not mate arm hole to torsion bar diameter 
or allow it to bend. When arm is in place, torsion bar support plate 
may be added to assembly. An error of 0. 005 in. to 0. 010 in. will be 
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absorbed by glue line. Use sufficient epoxy to form a barrier at 
inboard end of enjoinment. Note : Refrain from using excessive cement; 

otherwise, it may drip down on sensor. Once lower end is sealed off, 
remaining void may be filled ana surplus removed with razor blade. 

Cure at room ambient for min of 16 hr. Very carefully lower arm 
onto register and dowel, and secure. 

Note : Locking screw and supports must be in place before 

at taching arms or rotating sensor assembly . Rotate plate 180° and 
perform like operation to other arm. When both arms are secured, 
orthogonally measure clearance of adapter arm post to arm. Note : 

Make certain weight of transducer assembly is resting on lower arm. 
Measure each post clearance and divide shim into two equal parts; a 
±0. 0005-in. difference is allowed. Distortion of transducer assembly 
must be held to absolute minimum. Total value of each pair of shims 
should not exceed ±0. 000 5 in. of measured gap. 

III. FINAL ASSEMBLY 

Install capacitor drive. Set gap at 0. 020 in. ±0. 002 in. with 

sensor arms locked Release arm lock and check clearances. If 

error appears, reset stops until released position is same as locked 
position. 

IV. DISASSEMBLY - REMOVAL OF TRANSDUCER ASSEMBLY 

Orient sensor with housing cover up and release transducer 
from both arms. Remove bolts from torsion bar and arm. Block arm 
above register with 4. 5 -in. blocks after removal of arm locking 
brackets. Exercise caution in removal of screws in torsion bar case 
flange. Release torque from two screws in torsion bar support and 
leave sufficient torque to hold position. Remove remaining screws. 
Remove original two screws without exerting force on torsion bar. 
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Insert two razor blades, diametrically opposite each other, 
between torsion bar and case. Tap lightly if necessary. Lift sufficiently 
to insert two screwdrivers or case jacks and remove. Caution: Do not 
apply any force on support. Secure end of torsion bar to support with 
strip of tape and place in a secure storage. Disconnect transition wires. 
Remove housing cover. Transducer now may be safely removed from 
housing. Normally, it is not necessary to disturb opposite arm. If 
required, the same procedures apply. 

V. REASSEMBLY - REPLACE TRANSDUCER ASSEMBLY 

Clean cavity of housing. Insert transducer assembly, replace 
cover, and remove tape from torsion bar. With blocks in place, lower 
arm and torsion bar into place. Use light hand pressure to push 
torsion register into place. Replace screws in torsion bar support and 
lightly torque into position, avoiding as much as possible any preload. 
Secure torsion bar flange to housing. Lightly tap to provide nominal 
preload condition. Torque all screws tightly. Replace arm and secure 
to flange. Secure transducer to arms. Use same procedure as in final 
assembly. Reconnect t ransition wires and re-establish capacitor 
drive gap. 
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ABSTRACT 


Improvements in the knowledge of the earth's gravitational field would 
be of great value to the scientific fields of geology, geodesy and geophysics, 
and would have significant applicatic to orbital mechanics, navigation, guid- 
ance, earth dynamics, and mineral prospecting. The usual techniques of sur- 
face surveys with gravimeters have been augmented in recent years with 
orbital surveys using data derived from doppler tracking of satellites. These 
techniques have improved our knowledge significantly, but both are entering 
the regime of diminishing returns in terms of cost, time and effort required 
to obtain new data. Gravity gradient instrumentation for directly measuring 
the variations in the earth's gravity field from orbit has been under develop- 
ment for a number of years, but has not yet emerged into flight hardware 
status* Computer studies and experimental demonstrations with real gravita- 
tional fields in the laboratory have shown that a rotating gravity gradiometer 
flown in low earth orbit would be able to significantly improve our knowledge 
of the earth's gravitational field. Such instrumentation, flight proven for 
near-earth orbit mission, would also have wide application for measurement 
of the gravity fields of the moon, the terrestrial planets, the finer details of 
the outer planets and their satellites, and would significantly improve the 
scientific return from Comet/ Asteroid Rendezvous and Docking (CARD) mis- 
sions. The objective of this project in the AAFE program is to study the 
instrumentation and mission parameters required for improved measurement 
ef the earth's gravity field from near earth orbit, design and fabricate a sensor 
optimized for this mission regime, and show by operational demonstrations of 
laboratory prototype hardware that the instrumentation requirements are 
attainable. The work on the project commenced on 22 June 197 1. The mission 
studies are nearing completion, and the design parameters for the sensor are 
emerging. This paper is a preliminary report of the studies to date, combined 
with a summary report of relevant experimental and analytical studies carried 
out prior to the c ommencement of the project. 



ROTATING GRAVITY GRADIOMET ER 


Method of Operation 


The rotating gravity gradient sensor that has been developed at the 
Hughes Research Laboratories is a device for measurement of the second 
order gradient of the total gravity potential field. ^ The sensor configuration 
consists of a resonant cri aform mass-spring system with a torsional vibra- 
tional mode (see Fig. 1). In operation, the sensor is rotated about its tor- 
sionaliy resonant axis at an angular rate w which is exactly one-half the tor- 
sional resonant frequency. When a gravitational field is present, the 
differential forces on tne sensor resulting from the gradients of the gravita- 
tional field excite the sensor structure at twice the rotation frequency.^ The 
differential torque AT between the sensor arms at the doubled frequency is 
coupled into the central torsional flexure. The strains in this flexure are 
sensed with piezoelectric strain transducers which provide an electrical output. 

Since the rotating gravity rradiometer moves through the gravity 
gradient field and obtains a continuous sample of the field components in its 
plane of rotation, the output of the gradiometer contains two independent 
measurements of certain components of the gravity gradient field tensor. The* 
two measurements appear as two sinusoidal signals in quadrature 

,2 

*T=-^ l(r xx - r yy ) cos 2«t + 2r xy Sin 2ot] . (I) 

One output is a measurement of the difference between two of the diagonal 
components and the other measures the cross product component of the gravity 
gradient tensor in the coordinate frame of the sensor. Alternatively, the output 
can be represented as an amplitude which measures the difference between th** 
principal components of the gradient tensor a phase angle that represents 
the orientation of the reference frame of \\ > ipal gradients to the sensor 

reference frame. 


ATs£ t «tt -r„ „»««« 


For the simple ». >.amplf shown in Fig. 
by the mass M at the distance R is 


?.a) . (2) 

r. differential torque induced 


AT , i GM_ffi£_ 2ut 
4 R } 


(3) 


The angular resonant deflection between the two quadrupoles of the sensor 
rotating at one-half its torsional resonant frequency with an associated quality 
factor Q is therefore 


p 


2 ATQ - GM Q . , . 

• j - 3 — r s, -= sin cut , 

I (2u») IT (2«r 


(4) 
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Fig. 1. A Opc-atiOi >i Torsional Gravity Grrdiometer. 
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where I = mi / 2 is the quadrupole inertia. The angle 0 is extremely small. 
Surface gradients produced by the earth will produce angular deflections of 
B ' 3 x 10"6 rad in typical orbital torsional sensor designs, while useful 
threshold signals of 10"11 gal/cm (0.01 Eotvos unit (E.U.)) produce angular 
responses ot 10'^ rad. 

Although the. ~ reflections are small, they are easily measured by ^ 
utilising piezoelectric strain transducers attached to the torsional flexure. 
The threshold deflections of 10“^ rad produce voltage outputs of 10“® V from 
typical transducers. These voltage levels are easily measured by modern 
amplifiers. 


Present Development Status 


The ultimate objective of ou r work v. rotating gravitational gradient 
sensors is the development of a class of rugged sensors cf high sensitivity 
and precision which may be used to measure accurately and rapidly the details 
of the gravity field during airborne or orbital surveys and as a component in 
an inertial guidance system fo remove the effects of gravitational anomalies 
on the ultimate system performance. 

The objectives of the research program initiated in 1962 was to investi- 
gate the engineering feasibility of the basic concept, to develop sensor struc- 
tures which would operate a high sensitivity level both in free fall and in 
1 G environment, tc measure the censor's sensitivity to gravii aiional fields, and 
to investigate the sources of noise produced by the relati r 1 op-tne sensor. \ 
torsionally flexible structure utilizing piezoelectric rcadwafwas found to be a 
suitable design and offers a significant improvement over other possible gradi- 
ometer designs. It has demonstrated the capability of being operated in an 
earthbound laboratory environment (see Fig. 2) while still maintaining a high 
sensitivity and low signal-to-noise ratio. The present noise level of this sen- 
sur is tl E.U, (U at an integration time of 10 sec) and it limited by back- 
ground noise in the laboratory. Using this sensor, we recently arried out an 
experimental simulation ,here we measured in real time gravity gradient 
fields that had exactly the same magnitude and time variation as the gravity 
gradient signals that would be expected from an orbiting vehicle around the 
moon. * 


Figt e 3 shows the output of the sensor curing the passage at 60 cm 
distance of two masses spaced 62 cm apart. The first mass was 14.4 kg and 
the second was 15. 5 kg. The two gravity signals are resolved and if is possible 
to see the magnitude diffe. ence in the two masses. The three curves are the 
total signal amplitude, the sine output, and the cosine output. The dotted lines 
are the predicted outputs from tne computer simulation. 

At the present time, the development effort on the sensor is heading in 
two different directions. One program, sponsored by the Air Force, is for the 
development of an airborne gravity gradient measurement system. The major 
thrust of the development effort is to design a suitable hajrd mounted bearing 
and drive system that will spin a compact sensor (10 cm diameter) at the 
desired speed without * troducin* excessive amounts of noise and to design a 
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vibration and rotation isolation system that will isolate the sensor system 
from the aircraft noise and motion. The noise level of our present system 
(see Fig. 4) is about 6 E.U. with a 30 sec integration time. ' The goal of 
the program is to develop a gradiometer system capable of measuring gravi- 
tational gradients at the 1 E.U. level with a 10 sec integration time on a 
moving base, such as an aircraft or submarine. A similar design is being 
considered for deep space planetary and CARD missions. ^ 

The other program, under NASA sponsorship through the AAFE pro- 
gram, is for the design of an earth orbiting gravity gradient measurement 
system.-’ For the orbital case, the optimum method is to fabricate a sensor 
with a relatively low resonant frequency (2 to 8 Hz) and very long arms, attach 
it directly to the spacecraft and spin the spacecraft itself at the desired spin 
speed (1 to 4 rev/ sec). This mode of operation has two significant advantages. 
There are no bearing noise problems that are the primary source of difficulty 
in earthbound operation, and most important, since the spacecraft is rotating 
along with the sensor, the gravity gradient field of the spacecraft is stationary 
in the frame of reference of the sensor and the sensor does not sense the 
gravity field of the spacecraft, only the gravity gradient field of the earth. We 
have fabricated some prototype designs (Fig. 5) optimized for lunar orbital 
use and calibrated the sensor with ac gravitational fields from a dynamic 
gravity field generator. The measured noise lever* of 5 E.U. at 30 sec was 
limited not by the sensor but by the difficulty of isolating room vibrations at 
this low frequency (2 Hz), The objective of the AAFE program is to develop 
a sensor system optimized for earth orbit applications and capable of measur- 
ing gravitational gradients at the 0.01 E.U. level with a 35 sec integration 
time. 


Gradiometer Noise Limit 


The fundamental sensitivity of any sensor is determined by the thermal 
noise limitation. For the past ten years, we have been developing gravitational 
sensors working near their thermal noise limit. 8 Because this basic limit 
is dependent upon energy considerations, its calculation depends only upon very 
general parameters of the sensor, such as its temperature, mass, effective 
length, and time of integration. The results can then be applied to all gravity 
gradient sensors, regardless of their detailed design. 

In our torsional sensor the thermal signal-to-noise energy ratio is 
obtained by comparing the gravitational gradient signal energy in each of the 
two orthogonal outputs of the sensor to the 1/2 kT of thermal noise in each of 
the two degrees of freedom of the sensor. 

The minimum gradient that can be measured in each output of a ther- 
mally limited sensor with an effective arm radius r and total effective mass m, 

is; 


r - £M- _i_ /*L i/2 

1 ~ r 3 " rirr 'm' 


(5) 
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where r = 2Q/u is the 1/e time constant of the sensor and r\ is the ratio of the 
quadrupole inertia to the total inertia of each arm. 


i 


For what might be the desired sensor for earth geodesy, 3 one with a 
total arm mass of 8 kg (2 kg each), an r) =0. 95, and an arm radius of 40 cm, 
the thermal noise equation gives us a kT limit of 0.007 E.U. for a 35 sec inte- 
gration time. We should be able to obtain a measured noise very close to this 


thermal noise limit with a properly designed structure and electronic matching 
circuit. 5, ' * 8 
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EARTH GRAVITY MEASUREMENTS WITH ORBITING GRADIOMET ERS 


Harmonic Representation of the Earth's Geoid 


An objective of geodesy is to determine the variations of the earth's 
gravitational potential, which can be expressed in terms of spherical harmonics: 


r M oo . .n oo 

V = — — 1 + £ (“) \ P (sin f ) { C cos m\ + S sin m\} 

r L n=2 ' r/ m to nm nm nm J 

where a is the mean radius of the earth, P n m i 8 the normalized Legendre 
polynomial, C nm and S nm are the coefficients of the harmonic terms, and 
(r, p, X.) are the coordinate positions of the instrument. 

In the present satellite geodesy programs, orbital perturbation methods 
of obtaining the gravitational potential harmonics have led to the determination 
of the harmonics through the fourteenth degree and order. 9 l n theory, this 
technique can be extended to obtain all higher orders of the gravitational 
potential; however, it is anticipated that it will be difficult to obtain the higher 
order components. 

The advantage of gradiometer techniques in obtaining the higher order 
harmonics of the earth’s gravitational field is straightforward. Terms with 
increasing n correspond to small scale features on or near the surface. 
Although the contribution of these harmonic components to the gravitational 
potential is quite small, their contribution to the gravitational force gradient 
at a point above them is a substantial fraction of the gravitational gradient of 
the entire earth. 


A typical term in the gravitational potential 

v --5M te\ n+1 

V nm a \r/ 

gives rise to a radial gravity of 


P (sin f) C cos m\ 
nm nm 


^ ^ C n,n 


n+2 


and a radial gravity gradient of 

r rr . 0 = <„♦!> (n+2) £* (i) " t3 P nm C Bm co. mV . 


(6) 


( 7 ) 


( 8 ) 
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For the radial doppler velocity we take the time integral of the radial 
acceleration 

Av = f g dt = — — (-) P C sin m\ (9) 

r J ®r n av \rl nm nm 

where we have used the fact that the maximum spatial periodic variation (m=n) 
has a time variation due to the orbital velocity v given by 

cos m\ = cos n\ = cos ( n — | = cosi— t) . (10) 

max V a / \a / 

If we assume that^the strength of the harmonic components follows the 
statistical law 1 * S nm ~ C nm ~ l0~->/n^, and that (2n+l) terms contribute to 
each order, we can calculate the doppler velocity, gravity and gravity gradient 
as a function of the harmonic order. These are plotted in Fig. 6 for 250 km 
altitude. 

The doppler velocity data in Fig. 6 are correct, although they differ by 
two orders of magnitude from what would be calculated from Fig. 5-7, page 
5-28 of the Williamstown report, as presently published. In recent corres- 
pondence, William M. Kaula has brought attention to the fact that the right hand 
ordinate of Fig. 5-7 in the Williamstown report should read 10" 2 mm/sec 
rather than mm/sec. 

Figure 6 gives the rms strength of the signals at each harmonic order 
averaged over the entire orbit. If we assume a maximum sensor integration 
time of about 30 sec, which is determined by the time it takes the measurement 
system to pass through one resolution element (250 km), then there will be 
approximately 170 independent samples of each of the higher harmonic orders 
per orbit. The accuracy of determination of the rms harmonic component 
would therefore be the measurement system sensitivity for the 30 sec interval 
divided by 13, the square root cf the number of samples per orbit. Thus, a 
0.01 E,U. gradiometer system would be capable of better than an 8* measure- 
ment of the higher orders out to degree 75. 

Figure 6 indicates that if satellite-to-satellite doppler tracking tech- 
niques attain their anticipated sensitivity level of 0.05 mm/sec at 30 sec, 
doppler tracking will be able to extract gravity data up to degree 50, and if a 
gravity gradiometer with an 0.01 E.U. sensitivity at 30 sec can be flown it 
will contribute significant information out to degree 75. The comparative 
signal-to-soise of the two techniques crosses over at degree 35. We thus see 
that the two techniques are complementary rather than competitive since below 
degree 35 doppler tracking has a better signal level while above degree 35 the 
gradiometer gives better data. 

The average strength of the higher order gravity variations predicted 
in Fig. 6 and Fig. 5-7 of the Williamstown report use a statistical model 
based on the autocovariance analysis of a large variety of samples of gravim- 
etry, 12 and are estimated to be correct within ±30%. A statistical model 
assumes that the phases of the various harmonics are not correlated, whereas 
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we might expect some correlation in phases to occur at the position of signifi- 
cant geophysical anomalies, such as mountain ranges. To obtain some feeling 
for this possibility, we have also looked at the gravity fields to be expected at 
altitude for reasonable mass anomalies on the surface. 


Localized Anomalies 


Although a statistical analysis of the gravity field contributions of the 
various harmonic orders can give us a general idea of the rms signal strengths 
averaged over an orbit, it does not give us a good picture of the time history 
of the signals expected over specific anomalies where the phases of the har- 
monic orders combine to produce an impulse type response. In an attempt to 
get a better feeling of the signals to be expected from localized anomalies, a 
massive disc mast model was used to generate gravity data and the signals 
expected for both a single disc and a periodic array of discs were calculated. 

For the single disc model we chose a disc radius of 150 km or disc 
diameter of 300 km. The disc mass was chosen so that the gravity at the 
surface was 10 mgal. A plot of the results is shown in Fig. 7, which indicates 
that a disc with diameter 300 km, thickness 15 'an, mass of 1.7 x 10^ kg 
and density difference of 0.016 gm/cc will create at an altitude of 250 km the 
following signals: 

• Vertical gravity of 1.5 mgal peak 

• Vertical gravity gradient variation of 0. 11 E.U. 

• Vertical doppler velocity shift of 1.0 mm/ sec. 

The analysis of a single disc is, however, not a close analogy to the 
periodic variation in the gravity field that is implied by the usual harmonic 
representation of the field. The disc model analysis was therefore expanded 
to a calculation of the signals expected over a periodic array of positive and 
negative disc anomalies. The mass (positive or negative) was assumed the 
same as in the single disc analysis. The curves in Fig. 8 are extrapolations of 
the data from the center portion of the disc array to eliminate end effects. The 
periodic signals were 

• Vertical gravity ±0.65 mgal 

• Vertical gravity gradient ±0.1 E.U. 

• Vertical doppler velocity shift ±0.08 mm/ sec. 

We notice that the gravity gradient magnitude is almost the same as for 
the single disc. This is because the gravity gradient signal, being the spatial 
derivative of the acceleration, has a sharp cutoff, and the signal from an adja- 
cent disc of opposite mass acutally contributes slightly to the total signal. The 
magnitude of the periodic vertical gravity signal has decreased slightly from 
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Fig. 7. Gravity Signals of Single 300 km Diameter Disc at 250 km 
Altitude. 
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Fig. 8. Gravity Signals at 250 km Altitude from Periodic Disc Array. 
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the single disc signal. The broad signature of the vertical gravity signal 
causes signals from adjacent discs of opposite sign to partially cance'.. 

Finally, notice the very large decrease, over an order of magnitude, 
in the vertical doppler velocity signal from the single disc case to the periodic 
disc case. This is because the doppler velocity signal is the integral of the 
acceleration signal and the integration tends to smooth out the periodic varia- 
tion that we are looking for. 
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MISSION AND DESIGN STUDIES 


Introduction 


We have found in our research that the design of a gravity gradiometer 
is often strongly dependent upon the particular mission and using vehicle. This 
is especially true ir this application to earth geodesy. The size and opera- 
tional parameters of the sensor are strongly determined by the orbital altitude 
and inclination, mission lifetime, .and the requirements of geodesy. The sen- 
sor, in turn, has an effect on the spacecraft, especially the requirements for 
spin speed and attitude control. We have carried out some preliminary mission 
studies, and are beginning to define a set of mission, sensor and spacer raft 
parameters that are compatible with the overall mission objective of t lie accurate 
measurement of the earth's gravity field. 

The presently envisioned mission would use a Scout with a *12 inch pay- 
load shroud to launch two orthog'" illy oriented, spin stabilized satellites into 
a 330 km circular polar orbit some 1 5 -20 days before the vernal or autumnal 
equinox. Each satellite would carry an 8 kg, 80 cm. diameter giadiometer with 
a sensitivity of 0.01 E.U. at 35 sec integration time. The orbital lifetime 
would be short, but during that time we would obtain at least two complete maps 
of the vertical gradient and the horizontal gradients, both along and across the 
orbital track, with a resolution of about 270 km (540 km wavelength or degree 7^). 


Orbital Parameters 


To map the higher order harmonics of the earth's gravity lield, it would 
be desirable to have the measurements take place at as low an orbit as possible. 
Because of the mathematical characteristics of the potentU 1 field, the resolu- 
tion of any gravity measurement at altitude is roughly equivalent to the altitude. 
This is because the field strength of the gravity contributions due to the higher 
orders begins to fall off exponentially with altitude when the half wavelength 
becomes less than the altitude. However, a low orbit has a very short lifetime 
due to atmospheric drag, and a short lifetime makes it difficult to obtain the 
complete coverage of the earth that is also desired. 

What is needed is a low orbit with orbital parameters such that the 
orbital tracks interleave so complete coverage is obtained in a period shorter 
than the orbital lifetime, and where the track spacing is matched to the „wath 
width (equal to the altitude). It turns out that there does exist a set of orbital 
parameters that fits these requirements fairly well. At an orbital altitude of 
2 70 km, there exists what is called an "integer orbit.' 13. M The orbital track 
repeats upon itself after exaitlv 16 orbits. This can be a polar orbit, with 
16 orbits per siderea? d^y « r a sun synchronous orbit (at a slightly different alti- 
tude and inclination) with 16 orbits per solar d%y. If the altitude is slightly 
higher or lower, then the orbital track drifts so that the 16th orbit is displaced 
to one side or the other of the first track. These offset orbits finally begin to 
repeat after a number of days when the drift nas caused the satellite track to 
overlap the second ground track. There are two of these orbits that are of 
interest to us. They repeat after about 5 days, and their track spacirg is 
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approximately equal to the altitude. One is a polar orbit at about 320 km which 
repeats after 79 orbits and the other is a polar orbit with altitude of 220 km 
which repeats after 81 orbits. The track spacing between the half arcs for both 
orbits is approximately 250 km, so that there is a good match between the track 
spacing and the swath width. 

In reality, the orbital altitudes will decay due to drag, so that these 
simple orbital path models will not be followed exactly. We presently envisage 
launching into a 330 km polar orbit and al. owing the altitude to decay through 
these two altitudes where we get overlapping coverage. We have chosen a 
polar orbit rather than a sun synchronous orbit in order to obtain full coverage 
of the earth and to provide for calibration points twice per orbit at the two poles. 
The orbital lifetime estimated for the mission io somewhere around 30-50 days. 
The time spent near 320 km would be long enough to obtain good coverage of 
the earth at that resolution (640 km wavelength or degree 62). As the altitude 
decreases, we will get better and better resolution. We should get a substan- 
tial amount of coverage at around 220 km altitude with excellent resolution 
(440 krn wavelength or degree 90) but we will lose some coverage due to the 
rapidly decreasing altitude and the fact that the track spacing at the equator of 
250 km is slightly larger tha r the sensor resolution. 


Non-Eclipse Orbits 


It would be desirable to launch the gradiometer satellites into a polar 
orbit of the earth which will not cause the satellites to be eclipsed by the earth 
throughout the mission. The advantages of the non-eclipse orbit are the weight 
reduction and reliability increase available by elimination of batteries for elec- 
trical power during the eclipse portion of the flight. Also, the thermal control 
system required for the sensor would only have to contend with one state of 
thermal equilibrium rather than cyling between two. We have investigated 
possible non-eclipse orbits and find that even despite the relatively low orbits 
under considerations it is possible to achieve non-eclipse periods several times 
longer than the estimated lifetimes for these orbits. To attain these orbits 
only requires that a launch window constraint be placed on the mission. The 
satellite is launched 15-20 days before either the vernal or autumnal equinox 
(21 May or 21 Sept) into a polar orbit chosen such that on the day of the equinox, 
the orbital plane coincides with the terminator p’ane. At this point in time, the 
ecliptic poles and the celestial poles of the ear are all in the terminator plane. 
The slow rotation of the terminator plane about the ecliptic poles will cause a 
drift between the terminator plane and the orbital plane (M°/day). However, 
simple calculations show that with this choice of launch time and orientation, it 
is possible to have non-eclipse periods in excess of 30 days even for orbital 
altitudes below 250 kmlS, 

Non-eclipse orbits could also be chosen using sun-synchronous orbits 
lying near the terminator plane, and then the non -eclipse period is theoretically 
infinite. However, the above shows that with this minor constraint on launch 
time we can achieve the advantages of a non-eclipse orbit while keeping the 
self-calibration and full coverage aspects of the polar orbit. 
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Launch Vehicle and Spacecraft Configuration 


Because of the relatively low field strengths estimated for the higher 
order harmonics of the earth's field, it is necessary to make the sensor as 
large as possible. The arm length of the sensor is primarily determined by 
the maximum radius obtainable in the spacecraft, which in turn is determined 
by the payload envelope of the Launch vehicle. If we limit ouselves to a Scout 
launch vehicle in o. der to keep costs down and reliability up, we can use a num- 
ber of launch shroud configurations that have been developed for this vehicle. 
One of the largest in diameter is the 42 inch diameter shroud mentioned in the 
Scout users handbook. The allowable payload diameter for this shroud is 
96. 5 cm (36 in). The cylindrical portion of the payload envelope with this 
diameter is 84 cm (33 in) long, which would allow space for two cylindrical 
spacecraft 96 cm diameter by 42 cm thick. A very preliminary spacecraft 
design of this size is shown in Figs. 9 and 10. The sensor arm length with this 
spacecraft cohfiguration would be about 40 cm. 

The front part of the payload envelope could be used for a spin -up and 
attitude cobtrol system which would be used to insert the two spacecraft into 
orbit with the proper attitude and spin speed. After the payload had attained 
orbit, the spin control mechanism would increase the satellite spin speed to the 
desired rate (about 240rpm)and orient the spin along the orbital track. After 
release of one spacecraft, the jets wouLd be used to torque the other spacecraft 
so that its spin vector was perpendicular to the first spacecraft. With the two 
spacecraft in this relative orientation, one would be measuring the vertical 
gravity gradient and the cross-track horizontal gradient, while the other would 
be measuring the along -track horizontal gradient and a redundant measurement 
of the cross-track gradient. After 1/4 of an orbit, the orientation of the two 
satellite spin axes relative to the orbital track would change, and the data output 
from the two sensors would be interchanged. Although this is a relatively 
complex mode of data collection, it does allow for the measurement of more 
components of the gravity gradient at the same time. Most importantly, this 
mode of operation allows us to obtain the cross -track gradient information 
which can be used to tie the data together across the orbital tracks. The cross - 
track gradient along one track can be used to predict the gravity field at the 
next track. This closure property of the data sets can be used to eliminate 
drift errors. 

A simpler version of the experiment would be to launch a single satel- 
lite and to torque the spacecraft spin axis so that it lies in the plane of the 
orbit. The advantages of this mode of operation would be that the spin axis 
of the spacecraft would not change orientation with respect to the orbit and 
the drag torques would remain constant. In this orientation, the gradiometer 
would measure the difference between the vertical gradient and the along-track 
horizontal gradient and their orientation with respect to the local vertical. 


Sensor Parameters 


Most of the sensor parameters have been determined by the mission and 
spacecraft studies to date, although a few remain to be determined. The sensor 
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arm length has been chosen at 40 cm radius (80 cm diameter) as the largest 
arm radius possible for the 96 cm spacecraft diameter, which in turn is dic- 
tated by the Scout payload envelope. The arm end masses have been chosen at 
2 kg each as a reasonable weight for the size of the sensor and the payload 
capability of the Scout booster. The sensor time constant has been chosen at 
35 sec by using the time required for the spacecraft to pass through one resolu- 
tion element at the nominal altitude of 270 km and orbital velocity of 7.75 km/ 
sec. This figure was chosen as a reasonable optimum between the 41 sec for 
320 km altitude and the 29 sec for 220 km altitude. This sensor system time 
constant is the smoothing time that will be used in the sensor data preprocessing. 
The sensor output will be sampled more often than this (every 1 -5 sec) to over- 
come digitalization noise, to prevent aliasing, and to pick up strong, short 
period signals due to dense localized anomalies. 

The major sensor parameter that is yet to be determined is the sensor 
resonant frequency. We are presently .ompleting a study of the interaction of 
the spacecraft dynamics with the sensor dynamics and the choice of the resonant 
frequency will probably be determined as a result of optimization studies carried 
out using this dynamics interaction study. It is suspected, however, that the 
frequency of the sensor will be about 8 Hz, which implies a spacecraft spin speed 
of 240 rpm (4 rps). 
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CONCLUSION 


l 


We have carried out preliminary studies of the application of the 
rotating gravity gradiometer to the measurement of the gravitational field 
of the earth from orbit. Although the studies are still in their preliminary 
phases, a sensor design and mission profile have emerged that show promise 
of substantially improving our knowledge of the earth's field. The studies will 
continue with the objective of determining the optimum sensor, spacecraft and 
mission parameters. These will then be used in the design, fabrication and 
test of a laboratory prototype to show that the instrumentation requirements 
are attainable. 
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ABSTRACT 

Prior analysis has established that error signals gener 
ated in a Hughes Rotating Gravity Gradiometer (RGG) mounted 
in a spin stabilized satellite can be classified as three 
types: angular rate errors, angular acceleration errors, 
and mass unbalance errors. Equations of motion are developed 
for a dynamic model of the sensor satellite system, and errors 
are evaluated using parameters developed for an Earth orbiting 
experiment. The errors are shown to produce inaccuracies less 
than 0.01 E.U. for such a system. 



T 

l . 


INTRODUCTION 


The purpose of this report is to present a more general, 
improved version of the dynamic analysis of the Hughes rotat- 
ing gravity gradiometer mounted in a spinning satellite in 
planetary orbit. The original Hughes analysis was published 
in Research Report 441, Dynamic Analysis of the S e cond Order 
Grad iometer , by R.W. Peterson, July 1971. 

This analysis is preceded by a discussion and justifi- 
cation of the changes and generalizations stated in the 
Peterson report. It then completes a dynamic analysis with 
the new assumptions, makes some numerical assumptions about 
the sensor-satellite configuration, and develops numerical 
values for the errors generated by each of the dynamic 
sources . 

Evaluations of aerodynamic torque and aerodynamic drag 
are also made and included in the analysis. 

In conclusion, it is shown that dynamic interaction 
problems may be overcome with judicious control of sensor - 
satellite parameters. 
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TECHNICAL DISCUSSION 


II . 

A. Point of Departure from Previous Analyses 

Peterson's approach to the evaluation of the effects 
of satellite dynamics on the function of the gradiometei: 
appears to be a very useful analytical road. However, some 
of the assumptions made are open to question. 

1 have accepted Peterson's basic model of the gradiom- 
eter, the assumption of infinite rigidity about the k axis, 
and the ignoring of cross products of inertia. But I must 
disagree with the form of the inertia tensors for each gradiom- 
eter arm which Peterson proposes. 

Peterson has: 


I 1 

0 

0 


0 

0 

0 

0 

0 

0 

and $2 = 

0 

*2 

0 

0 

to 

0 

*1. 


0 

to 

0 

V 


Equations (3) and (4), RR 441 


This form of the principal inertia tensor assumes that each 
arm is dumbbell-like as in the following diagram: 


1405- I 


kk 
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However, even the dumbbell model fails to be represented 
by Peterson's suggested tensors, because the principal moment 
about 57 for Arm 1 is neglected, and the principal moment 
about ii for Arm 2 is ignored as well. 

I suggest that the principal inertia tensors for the 
arms be kept as general as possible, in line with the diagram. 

1403- 2 

kk 



conformity with Peterson's notation but use superscripts tc 
indicate principal axis: 



I* 0 0 


r 

o 

o 

.H <N 
M 

1 

Arm 1 = 

0 1^0 
0 0 

Arm 2 = 

0 l\ 0 

0 0 1^ 


Peterson's* coefficients become: 


a 


1 




a 


2 




12 


D 1 + D 2 

7 * 77 * 

1 1 l 2 




+ K 


+ I 


2 

3 

2 


*RR 441, p. 6. 
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b> 


"0 


D, 


u>‘ = K, 


IV + 1 3 1 
• L 1 1 2 

1 

D l + ° 2 ' 

I 3 I 3 

T 

1 ? + 1 1 

L 1 

2 J 


L 1 2J 

r -» 



- 

1 J + 1° 
• L 1 2 

4 . 

K 1 + K 2 

3 3 

T J 1° 

L 2 J 

7 

I 3 + I 3 

ll 2J 


Peterson's equivalent gradient signal (eq. (20)) becomes 
(with the new tensor forms) : 



The most significant changes come in through the torques 
M 2 which must now be rewritten. 

Pecerson's eqs. (15) and (21) become: 


M 1 

» 

\M 

II 

V 

+ E * 

M ul - 

[I ?“k + *1 “i»j 



M 2 

= k • 

M g2 

+ 1c • 

"u2 - 

ll 2“k + I 2 “i"j 

- *2 



M i 

3 


M 2 

7? 

~2 


r ij + r wi ' W k + 


r ij + P y2 ' “k + 



W i W j 
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where 


Let 


r*. 

13 


k • M 


3k 




1c • II 


3l 


yi 


k • M 


A 


y2 


k • H 


V2 



and 



Equation (2) is the resulting expression for r^. This 

equation is analogous to Peterson's eq. (23) . As in Peterson's 

analysis we take the mass unbalance terms to be of opposite 

sign for the maximum error (i.e., r,«-r„-r ). How- 

1 yl 2 y2 ymax 

ever, I see no reason to make r^j * - Ty as in Peterson's 
analysis. 


r 




— — s — a v 

<■ ♦ « l2 s ♦ eJ 2 j [s 2 ♦ £ s ♦ «• 


Gravity Gradiant Signal 


X Y 2 ~ c Vl ) * * B 1 Y 2 


Rotational Piald Error 


Yj) + S(« 2 t x - 


tton Unbalanca Error 

♦ |* a u> ♦ »u 2 ♦ « 2 > ♦ <», ♦ )J r | i _, 

Angular Accal oration Error 

♦ [ («! - « 2 > » ♦ <•? * »*»] -J 


8 :v]'tj 

•l v 2*] “i“j 
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( 2 ) 



The three gradient error sources (as identified by Peterson) 
are: 

• Rotational field w.w. 

i D 

• Arm mass unbalance r 

ymax 
2 2 

• Sum mode mismatch (3 1 - (J_)<I>. acting on angular 
acceleration input <i£. 

Only w. affects F exactly as in Peterson's analysis. Note 
k e 2 2 

also that a dc error of the form (f^y^ ” ^l Y 2^ now a PP ears 

as a result of having y^^ f 1 and y 2 ? 1. There is still a 

potential large error due to excitation of frequency & 2 * 

1 have followed Peterson's basic approach of estimating 
the satellite motions in the principal inertial axes of the 
satellite and then transferring the effects of these motions 
to the sensor in terms of the orientation of the sensor with 
respect to the principal inertial axes of the satellite and 
with respect to the center of mass of the satellite. But 
I have radically altered the expression for w.w. and J> v 
found by Peterson by incorporating essential changes in 
the dynamical analysis. By far, the most fundamental change 
is to include the effect of the satellite's precession. 

I believe Peterson has inadvertently eliminated this effect 
by calling the nodal angle zero (an angle which he calls 
p. 10 of RR 441) . In order for the satellite to precess, 
the nodal angle, defined by the intersection of the satellite's 
equatorial plane with the (approximately) invariable momentum 
plane and by some arbitrary direction in the momentum plane, 
must be cyclic at the precession frequency. 

Peterson also makes the assumption that a special 
relation holds among the satellite's principal moments, 

A, B, and C. He assumes k^ » k 2 or (C-B) /A - (OA)/B. The 
latter implies C « A + B, which seems quite arbitrary. Peterson 
later estimates k^ * (B-A)/C « 0.1, which, with his former 
assumption, implies A ■ 0.45 and B » 0.55 if C is 1.0. I 
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r 


have chosen instead to assume general principal moments, 

A, B, and C even though I feel that choosing A = B would 
be practical enough. 1 have chosen a precession frequency 
expression based on the coning angle 0, spin frequency w , 
and principal inertial moments A, B, and C. Strictly speak- 
ing, a torque-free body exhibits uniform precession only if 
A = B (A and B are transverse moments of inertia) . However, 
since the body is in this case not torque-free (due to 
aerodynamic torque) and A will probably nearly equal B, I 
have confidence in estimating a uniform precession. Essen- 
tially, I have considered a uniformly precessing satellite 
and superimposed aerodynamic torque on its motion. 1 have 
also considered the effect of the linear acceleration caused 
by aero-drag on the mass unbalance error terms. In this 
analysis, I have chosen to ignore the solar radiation torque, 
the magnetic torque, or the micrometeorite impact torque. 

With the exception of the magnetic torque, calculations would 
show that the latter torques would be less than the gravity 
gradient and the aerodynamic torques in near earth orbit. The 
following points summarize my changes to Peterson's dynamic 
analysis of gravity gradient effects: 








Do not assume (C-B)/A = (C-A)/B. 

Assume general principal satellite 
moments C > B > A. 

Include effect of aerodynamic torque 
on spin: M = aerodynamic torque. 

Include satellite precession. 

Change to conventional Euler angle 
notation. 


Include linear acceleration from 
aerodynamic drag. 
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B. 


Dynamic Analysis 
1 . General Equations 


The local orbital coordinates are X, Y, and Z; 
Y aligned with local vertical, Z normal to orbit plane, and 
X tangent orbit, as illustrated in Fig. 1. 



Fig. 1. Local Orbital Path Coordinates. 

The principal axes of the satellite are x, y, and z 
which are related to the approximately inertial space frame 
axes X, Y, and Z by the Euler angles ip, 9, and (f> . The 
geometrical relationship of these angles and the axes is shown 
in Fig . 2 . 

The spin of the satellite, which is approximately 
constant at <I> , is about the z principal axis and describes 
the angular motion <fr. The Z axis is not only normal to the 
orbit plane but describes the constant direction of the 
satellite angular momentum vector for torque-free motion. 

The angle 9 is the coning angle of the satellite precession. 
For torque-free motion, 8*0. The angle ip describes the 
satellite *8 precession; 4* is the precession frequency. 

(Note: Peterson calls ip, $ and then considers $ ■ 0, implying 
no precessional motion.) 
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Z SATELLITE 
a EQUATOR 


ORBIT PLANE 
(INVARIABLE MOMENTUM 
/ PLANE OF SATELLITE 
( SPIN ALSO) 


Fig. 2. Euler Angle Relationrhip of Satel- 
lite to Inertial Space Coordinates. 


The standard transformation between the space axes 
and the principal inertia axes is: 


(coM coM * tin* cosd sind) (-sin* coM * sind cos8 coM) (sind coM) x 

(coM sind ♦ siM cosd coM) (-Bind sind ♦ cosd cosd cosd) (-sind cosd) y 


(sind sind) 


(sind cosd) 


The Euler dynamic equations describing the satellite 
motion under gravity gradient torque alone are contained in 
Peterson's eq. (26): 


• w + to x • to] “ 3Hq [Y x * Y) ] 


inertial angular velocity of the 
satellite in principal axes coordinates 
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( 


A 0 0 


0 B 0 


0 0 C 


= principal inertia tensor of the 
satellite 


R = satellite distance from center of Earth 

0 

140 = mass of Earth 

SIq = angular velocity of satellite in orbit 

a 0 « % 


2 

For surface orbit Rg = 6,371 km, ft* = 1.54 x 10 sec 

The dynamical regime of the satellite includes only 

small coning precession angles, 0 << 1. Hence, cos 0 ~ 1. 

We therefore can approximate Y: (Note that <J> = w t and 

s 

later that = (it.) 


Y = 


(cos u) t sinip + sinw t cosiji) x 
s s 


sin (w t + \j0 x 
s 


(-sin w t sin\|> + cost*) t cosiJO y = cos (w t + y 


-0 COS<|l 


-9 cosi|» 


Y x ($ g • Y) is found by the determinant: 


sin (u t + 4>) 

S 

A sin (w t + (J>) 


COS (u) g t + lj>) 

B COS(b) t 4* \|l) 


-0CO8(|» 
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i 


» 


1 


7 x (4> • Y) 

s 


(B - C) 0COS4) cos(co t + 4>) X 

s 

(C - A) 9cos^ sin(u t + iJj) y 

b 

(B - A) cos(co t + 40 sin(aj t + 40 z 
s s 



I 

[ 

i 

L 

[ 

[ 

l 

t 

t 

i 

i 

I 


u>) 


co + co x ( 4> 


The resulting dynamical equations are: 



(Ail 

X 

+ 

C03 03 

y z 

- Bu) u) ) 

y z 

X 

35 

(Bu 

y 

+ 

Aoj o) 

X z 

- C03 'U ) 

X z 

y 


(Cw 

+ 

3 

3 

m 

- AoJ 03 ) 

z 


z 


x y 

x y 

m 


co + k, co co 

x 1 z y 

co - k_oo co 

y 2 z x 

co + k_co to 

z 3 x y 


* - 3ftgk^0 cos4> cos(co g t + 40 
2 

= 3^2^20 cosijt sin(co s t + 40 
« (3/2) B^sin 2(co s t + 40 * 


where 


k 

k 

k 


1 

2 

3 



(3) 


To solve these highly nonlinear equations would oe 
impossible without some simplifying assumptions. I have 
solved the first two equations for co and u by first 

^ y 

assuming w approximately constant, co ■ co . Then the product 

Z Z 8 

oo co may be substituted in the co equation and co will be 
** j z z 

known. Because the solution for u is done after u and u> 

z x y 

are found, there is an opportunity to include the predominant 
effect of aerodynamic torque, which is a decrease in the spin 
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angular velocity. The u equation for gravity gradient and 

z 

aerodynamic torque is: 


(I) + k-u uj 

z 3 x y 


| «Qk 3 sin2(w g t + *) 


M 


where 


M = aerodynamic torque about z axis. 

After the approximation = us has been made, the x 

t g 

and y equations are in matrix form: 


s 

k, b) 

1 s 


U) 

X 


2 

-SQ^k^y cossj; cos(w g t + i|») 

; k 2“ s 

s 


“y. 


2 

3JJ 0 k 2 8 cosij; sin(w g t + s|>) 


where s = d/dt. 

Solving for u and u indicates that the equations 
x y 

to be solved are a second order linear differential equation 

in u and a second order linear differential equation in u . 

x y 

The homogeneous part of the solution is found from the simple 
solutions to the following equations: 

(S 2 + *i k 2 w s S)w x * 0 
(S 2 + *i* 2 w s s)w y - 0 

The homogeneous solutions are of the form 

u» x ■ Cj sin /lcpqu> a t + C 2 cos /kjkjW^t 
w y - C 3 sin /k^kjW^t ♦ cos ✓k^kjWgt 
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The above must also solve the original homogeneous equations: 


u) + k.u u = 0 
x 1 s y 


li “ k_U) GJ = 0 

y 2 s x 


To find the relationship between the coefficients of the 
homogeneous part of the solution, we combine the last two 
equations and obtain: 


C, /kTiTTtjj cos /k, k_ui t - C~ /k, k_a) sin /k. k~w t 
1 1 t s 12 s 2 12 s 12 s 

+ C-k..u> sin /k, k_u t+ C-k^ cos /k, k_ u t = 0 
31s 12s 41s 12s 

C 3 *^1^2 w s COS ~ ^4 ' / ^1^2 u s S ^ n ^2 w s t 

- C..k 0 u) sin /OTTw t - C .k-w cos /iTjcTw t = 0 
12 s 12 s i 2 s 12 s 


We find that 


C 1 "^1*2 + C 4 k l * 0 


^ 1^2 + C 4 *^ 1^2 = ® 


=> c 


- ’''*1*1 


and 


C 3^1 ~ ^2 = ® 


C 3 *^1*2 ~ C 2 k 2 


! > C. 


/snr 

* 1*2 


Thus , the resulting homogeneous solutions are: 


« x * C^sin /kpqu. s t + c 2 cos 

“y * — Up c l co * + — kH- c 2 ,tn /E I E 2“. t 
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and C 2 are arbitrary. The nonhomogeneous part of the 
solution for u> and oj may be determined by the method of 

^ y 

undetermined coefficients. The following assumed forms of 

w and u (nonhomogeneous) were uscjd successfully. (A, B, 
x y 

and C are not inertia tensor members.) 


a) = A sint sin(u t + t) + B cost cos(u t >- t) 
y s s 


w = C sint cos(u> t + t|>) + D cost sin(u> t + t) 
X s s 


(6) 


Differentiating the above terms of w and u (remember, 

• x y 

t * tt) : 


li = A- 1 ' cost sin(u t + t) ♦ A(u> + t) sint cos(w t + t) 

y s s s 

• . 

- Bt sint cos (u t + t) - B(u + t) cost sin(w t + t) 

5 S S 

• • 

u> * ct cost cos(u t + t) - C (w + t) sint sin(w t + t) 

s s S 

• • 

- Dt sint sin(w t + t) + D(u» + t) cost cos(w t + t) 

o 8 8 


Substituting the m , w . and w , w forms into the nonhonoge- 

x y x y 

neous differential equations for u> and u we get: 

x y 

• • e 

iS equation cos cos: Ct + D(w +t) + k. w B » - 30*k, 9 

8 X 8 U X 

• • 

sin sin: -C(w # + t) - Dt + k-jU A * 0 

<*> y equation cos sin: At - B(w + t) - k 2 w # D * 30 0 k 2 9 

sin cos: A(« s + t) - Bt - k 2 w # C - 0 . (7) 

The latter equations may be put in matrix form for easier 
visibility of the solutions: 
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0 

k l“s 

i 

( Ug + 

T* 1 


"- 3!, o k l 9 " 

Vs 

0 

-(u g + ip) 



B 


0 

* 

• 

-(U) g + !(/) 

0 

” k 2 U> 8 


C 

- 

3^ k .e 

“s + * 

4 

“Vs 

0 


D 


0 


2 

Adding the first two rows: (A + B) k.w + (D - C) u> = -3Q rt k.0. 

la soi 

Subtracting row 3 from row 4: (A + B) u + k 0 ai (D - C) = 

<5 S Z S 

-311^20 • -Solving for (A + B) and (D - C) : 


A + B 


-3nQ0< 2 (1 - k x ) 

-7 1 ” k l k 2^ 


D - C 


-3fl‘0 k. (1 - k.) 

— •7Tr--Tjkp- 


Adding rows 3 and 4: 


IA - B) (u> # ♦ 2 <ii) - k 2 w g (C + D) « 3njk 2 9 


(Sa) 


Subtracting row 1 from row 2: 

(A - B) kjW g - (u> # + 2i) <C 4 d) - 30^0 . 

Solving for (A -• B) and (C + 0) : 


A - B 


3ftjk 2 e [k^ - <u> # + 2*)] 
tk 2 k lW J - (w # ♦ 2i) 2 ] 


C D 


-30*^6 [k 2 u> g - (w § 2*) 1 

[k 2 k ^ - (u s 4- 2tJ) 5 1 


(9b) 
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2 . Estimate of Precession Frequency iji 

The nonhomogeneou^ parts of u and w involve 
• ^ y 

= 4>t as arguments of sines and cosines; hence, it is 
important to know . Even though our satellite is not 
torque-free, we can get an estimate of ^ by examining torque- 
free motion. 

With the assumption that 0=0, the angular velocities 
and accelerations in terms of Euler angles are: 


oj = sin9 sin<t> 

9 

Wy = sin0 cos<)> 

• • 

U) = <f» + \I) COS0 = U) 
z s 


(i = \(>4> sin0 cos<j) 

9 9 

<*>y = sin0 sin4> * 

u = 0 
z 


Taking one of the torque-free Euler dynamica 1 . equations; 


Ail) + u) u) (C - B) = 0 
x x z 


Substituting the values of <S , u , and u> above: 

XX z 

* • • • • A 

sin0 cos<t> + sin0 cos^ + y sin0 cosO cos<J)] (C - B) = 0 . 
The result is: 

♦ =* (^T^) 

Since 0 << 1 and 6 ■ u : 

8 

♦ * ( -B--" g fi ) (“I 
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I 


Note that for c > B > A, ij> <0 for > 0. This means that 
our satellite will undergo retrograde precession. 

3 . Evaluating the Product 

The product is present in the gradient 

error driving functions <I>, and to.w., so it is important to 

K x ] 

evaluate. 

First, we simplify the nonhomogeneous parts of and 
Uy with the basic trigonometric relationships: 


sin 


p + q 
2 


cos 


P ., -1 Si 

2 


j (sin p + sin q) 


sin 


E-l-S 

2 


cos 


E_L_a 

2 


= J (sin p - sin q) 


cos 


p_+. a 


cos 


E--_a 

o 


= J (cos p + cos q) 


sm 


E— t a 

2 


sin 


E -_SL 

2 


j (cos q - cos p) , 


if P " = ip and P = w t + ip then p = (u t + 2iJ>) 

Z 4 S S 


and q = -uj t 
^ s 


The simplified equations are: 

w = (cos uj t - cos (u t + 2ip) ] + y [cos ( u> t + 2^) 

Y 4 5 S 4 9 

+ cos w t] 
s 

w = £ (sinfu.t + 240 - sin u t] + E [sin (w t + 2i|>) 

X 4 S S 4 S 

+ sin w„t] . (11) 

s 

Thus (nonhomogeneous solutions) : 
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0) 


A + B . A - B , , , 0 , v 

r COSO) t - 7Z COS(U) t + 2ty) 

Z S • z s 


to = C — t - D sin(to t + 2\p) + - - sin oj t i 

X Z S 4 SI 

(12) 

These equations make use of the relationships between A, B, 
C, and D previously derived. 

Adding the homogeneous parts to the nonhomogeneous 
parts of and and solving for the product 0 ^ 01 ^: 


CJ U) 

x y 


_ (A + B) (C + D) 



4 


(A + B) 

(D 

- C) 


4 


(A + B) 
2 

C 1 

cos a) 

(A + B) 
2 

C 2 

cos u> 

(A - B) 

(C + 

D) ( 

4 



(A - B) 

(D - 

C) , 


cos a) t sin (u t + 2\b) 
s s 


cos w t sin a) t 
s s 


12s 


cos (w t + 2ip) sin (u t + 2 iJj) 
s s 


cos (w t + 2 iJj) sin u t 
s s 


- — - 0 - C. cos (w t + 2 i|i) sin /kTkTw t 

2 1 s 12s 

- — — C~ cos (u> t + 2^) cos /OTu) t 

2 2 S 12S 


/kTk 

~k 


K V 

^-2 (cos /k^kjWgt) D | sin (ai g t + 2i JO 


✓nr 

K 1 K 2 


C, ~ - -y - c .l cos /kXut sin to_t 


k x W 1 2 




s 


/kTk 


|-2 cJ cos /Kpkjw s t sin ✓k^kjw g t 
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4 


/k l k 2 


C.C. cos v'k. k. u t 
12 12s 


2k, 


C- (C + D) sin /k, k_oa t sin (go t + 2 ty) 
2 JL Z S S 


l/ hS 


2^ C 2 (D - C) sin /k^k 2 (jLi s t sin u g t 


+ — 


•'Vi 

V~ 


2 * / k 1 k_ _ 

C,C„ sin /k,k o 0) t + — =-=- C_ 

12 12s k, 2 


cos /k, k-uj t sin /k, k_u> t 
12s 12s 


( 13 ) 


By using the sine and cosine relationships already 

described, it is possible to reduct each term of the expanded 

w w expression to a recognizable frequency combination of 
x y 

sines and cosines. 

The full expansion of 2S: 

= (A _ + B) 8 , (C + ,_ P) . [sin2 (u g t + 41) + sin 2 ty) ] 


. (A + B) (D - C) ... _ (A - B) (C + D) 

+ sin 2u g t ~ § 


sin (2u t + 4i|>) - — P — q ~ ^ ~ — [sin2(w t + ip ) 

s o s 


- sin 2’|> ] - ~ 2 k ~ 2 - sin(2/F|Tc^ w g t) 


/kTk 


iq- 2 Ci c 2 oos < 2 /k I k I “ 3 « 


(A + B) 


C, [sin (1 + /k^kj) w g t - sin (1 - /k, k 0 ) w„t] 


1"2' s 


141 





l 


+ ih t c 2 [cos (1 + ^k^) o) s t + cos (1 - 


/i qk? o) s ti 


(D I C) C 1 [sin (1 ♦ A ? 9 V + sin (1 - 


✓k,k n (D - C) C. 


“s tl + 


[cos (1 - /k,k') w t 


cos 

(1 + ’ 

(A 

- B) 


4 

sin 

(24» + 

jA_ 

- B) 


4 

cos 

(24- + 

/k 

?2 (C 


12 s 

^ c 2 2 sin < 2 /^ 


[sin ( 24 - + (1 + /k^k 2 ) ^ s t) 


Jlrill [sin ( 24- + (1 + » 8 t) 

4k l 

+ sin (24- + (1 - /k^kj) u) s t ^ 

+ iSSlil * °— [cos ( 24 - + (1 - /k]k^) “ s fc > 

4K 1 

- cos ( 24 - + (1 + •'kjkj) Wgt) ] • (14) 

4. Angular Misalignment of Gradiometer; Calculation 


of u . w . and w. 

1 J 

Following Peterson's formulation, the gradiometer 
case axes T, J, and 1c are assumed to be canted with respect to 
the- satellite principal axes by the matrix relation: 
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a -e 


3 and a are small angles 

a << 1 
3 << 1 

For a better understanding of the geometry of this 
relation, I have separated the above matrix into a product 
of two matrices that perform the rotations of the x, y, and z 
axes into the i, j, and k axes shown in the following. 


a -3 


1 0 

0 1 

0 -3 


0 -al Wotes a 3 << 1 


*' i" 


ROTATION ABOUT Y AXIS ROTATION ABOUT X* AXIS 

The angular misalignment results in the following relations: 


u i 

* U) 

X 

- 

aw 

z 



** 

- 

aw 

z 


•n 
3 . 


+ 

3w z 

, and 


“y 

+ 

w. 

• 

u k 

N 

3 

n 


au x " 6a) y 


“k " 


+ 

adi - 3di 
x M y J 
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By multiplying 


we have the useful relation: 


2 

u).u. = u> at + u (@u - au ) - aB « 

lj xy z x y z 


where 


CO =0) 

z s 


( 15 ) 


From the previous expansion of to u> , we can list the ampli- 

x y 

tudes and frequencies of the rotational field driving func- 
tion j , as shown in Table 1 . 

Finding the sum mode mismatch driving function 
involves the two equations 


(16a) 

k.to co . (16b) 

3 x y 

The sum mode mismatch driving function therefore con- 
tains all of the frequency components that the rotational 
field function had, except their amplitudes are multiplied 
by the factor . However, it also contains the additional 
components listed in Table II. 


<I>, = <!) + ati - QiL 

k z x y 


14 


“z = I n 0 k 3 Sin 2{ *s t + + c 5 - ‘ 


C. Aerodynamic Torque About Spin Axis 

According to a technical analysis by K.R. Johnson on* 
the effect of aerodynamic torque on satellite rotation, the 
torque caused by the differences in velocities with respect 


•"The Effect of Dissipative Aerodynamic Torques on Satellite 
Rotation" by K.R. Johnson, Journal of Spacecraft and Rockets, 
Volume 5, Ho. 4, April 1968, page 408. 
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TABLE I 


Rotational Field Driving Function, 


oj . ui . 

i J 
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Additional Components of Sum Mode Driving Function 











to the atmosphere of different portions of a rotating satel- 
lite causes secular variations in the precession rate, spin 
rate, and coning angle. I have chosen to integrate into my 
analysis only the secular variation in spin rate caused by 
aerodynamic torque. The choice was only based on the desire 
for simplicity and the fact that the coning angle and 
precession rate change introduce ’’perturbations on 
perturbations. •' 

Johnson's eq. (44) indicates the time variation of 
spin rate. 

* - t/T (fc 

$ ~ 4* o® r 


where is the initial spin, and xs a parameter to 
be defined below. 

Hence , 


$ 



-t/T 

e 




M (aero) 
z 

C 


0 



Johnson indicates that his p formulation for a 
cylindrical satellite holds only when a parameter p, which 
he defines, is either greater than 2 or less than 1. 


P 



3R 2 A 
— 

+ r 0 )C 


K « radius of cylinder 

r n ■ half of satellite 

length 

A ■ transverse moment of 

inertia 

C * z axis moment of 

inertia . 
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For our proposed satellite designed to fit into the 
Scout upper stage, R = 38 cm, r^ = 20 cm. Hence, p = 2.36 
(A/C) . Depending on our A/C ratio, the analysis may or may 
not be valid. Assuming 

* > 

C 

an unlikely configuration due 
or 


we proceed to determine $. 

According to Johnson, 


3C 

R 2 ) <s in 6> av 

density 

velocity 

one precession period of the 
satellite spin axis and 
air flow 

determine 


(0 240 km altitude) 

r Q ■ 20 cm 

R ■ 38 cm 

C » 7 . 1 kg m 2 




4p V 0 r 0 R ( 


P - 


<sin 6> 


av 


upper atmospheric 

satellite orbital 

the average over 
angle between the 
the direction of 


I used the following values to 


♦ i 


<sin 6> 


av 

P 


- 4 Hz 

* 1 (6 * tt/2) 

» 2 x 10" 13 gm/cm 3 

■ 7.4 x 10 3 cm/ sec 


2 

2.36 ' 

to stability considerations, 
1 

2.36 
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i 


-4(2 x 10" 13 ) (7.4) x 10^) (7.6 x 10 2 ) (4.33 x 10 3 ) 

3(7.1) 


= 9 x 10 ^ rad/sec 3 = ~ 


Arm Mass Unbalance Driving Function 


T is produced by the linear acceleration of the 

gradiometer center of mass acting on the individual offsets 
of center of mass for each arm from the pivot center. The 
basic correction to Peterson's derivation which I would make 
is the inclusion of the effect of aerodynamic drag. Peterson's 
model has the gradiometer center of mass offset from the 
satellite center of mass by the vector T. Peterson's expres- 
sion for the linear acceleration of the senscr center of 
mass is then: 


a = u x H + u x [u x £.] 


I = and u = 


uj . (coordinatized in sensor 

case axes) 


The effect of aerodynamic drag would come in as an additional 


term: 


a»a D +uxl+ux (wx!) 


Peterson's expression for I*. needs no change: 

UZuaX 


? (k • (i x a) | 


Umax 1 
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I 


3 3 

Where I * 1^ * n 8 mass of one gradiometer arm, 
and 7 is the offset of the sensor arm center of mass from 
the pivot center. 

6 ." 
l 

Y 

a D is initially coordinatized in satellite principal axes 
and it must be put in sensor case axes: 

1 0 -°] IV 

0 1 a a y 

° l \ [ a 2_ 

The total linear acceleration of the gradiometer center of 
mass is: 




(* x - «*,» ♦ “j‘k • Vj + - w j l i> - ( “k £ i - “AV 

t» y ♦ ♦ v i - y k ♦ <v*j l k - 

(«* x - e» y ♦ y ♦ ♦ i“iVi ' “iV ’ ( “j‘k ' “j“kV. 


Peterson's analysis ignores the aerodynamic drag contribution 
to Syntax which i,: 

^ymax^Aero Drag " ^i^ a y + ® a z^ ” ^j^ a x ~ oa *^ ± 

The aerodynamic drag will be along the tangent to the orbit 
(approximately) in the direction of the negative X axis. 
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Ignoring precession and coning angle. 



Hence, 


cos k> t sin u t 0 

r-a..l 

s s 


X 

-sin w t cos u t 0 


0 

s s 



0 0 1 


0 

_ 


_ _ 


(T ) , = [6 (a„ sin u t) + 

Umax aero drag i X s 


V a x cosu, s tn r 


( 18 ) 


E . Graphical Evaluation of Error Terms 

Figure 3 plots the magnitudes of seven complex func- 
tions, formed by replacing S by jw versus frequency in rad/3ec. 

2 2 2 2 2 
Uq/Q is indicated, as are (S^ 9 - u> + jwa^) and (“g “ « + 

ju (w 0 /Q) ) . A combined curve of the three factors is formed. 

The three faccors below are also plotted as: 

2 2 2 

Mass Unbalance: ( @ ^ + - 2u> + jio(a^ + a^) 

2 2 2 

Rotational Field: “ ®1 Y 2* " w ^ Y 1 “ Y 2^ + jutc^i " a l T 2^ 

2 2 

Angular Acceleration: (0^ - & 2 > + (o^ - a^) jo 

In Fig. 4 the error numerators are multiplied by the combined 
curve to give tne dimensionless frequency re ponse . * To find 
the absolute error magnitudes at specific frequencies, it is 
only necessary to multiply the three curves by the magnitude 
of <L k , and respectively. 


Reference the Appendix. 



MMNtTUDC OF MULTIPLIER (UNITS VARY WITH CURVE) 


I40S-T 



Fig. 3. Complex Functions In Dyna- 
mic Equations. 





In order to construct Figs. 3 and 4, it was necessary 
to make numerical assumptions about the gradiometer system. 
While an effort was made to make reasonable estimates, the 
ideal gradiometer design process might include a computer 
evaluation of the errors under varying assumptions of 
parameters . 

1 . Design Assumptions 

a . Sensor Assumptions 

ij ^ 1^ = 4.5 x lO^gram cm^ 
g 

K q = 4.6 x 10 dyne cm/rad 

^ Kj = 4.6 x 10 6 7 dyne cm/rad 
Q (untuned) = 70 

= 45.4 rad/sec = 7.23 Hz 

2 -2 
Uq = 2060 sec 

Wq/Q = 0.65 sec”^ 

“0 

^ = 29.4 sec ^ 

6 2 

Dq = 1.17 x 10 g cm /sec (From Above Assumptions Plus 

D x -v D 2 = (1/2) D q and w Q /Q 
Expression, Section II -A. ) 
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5 2 

5.8 x 10 g cm /sec 


D, D„ = t: D n = 

1 2 2 0 


a l % a 2 ^ <X 12 ~ sec 


-1 


0 1 ^ S 2 % 0 12 = 10-2 sec ” 2 


^1 ^ ^2 % ^12 = 2,2 sec 


-1 


Yl % y 2 = °‘ 7 


b. 


Difference Factor Assumptions 


2 x 10 -4 (1 

g at 30 cm) 

(10 -2 ) (0.13) 

= 1.3 x 10 -3 sec -1 

-2 

-2 

10 z (10.2) 

= 0.102 sec 


6 1 ' S 2 

(a 2 Y 1 - a^) = 10 -2 (0.13) (0.7) = 9 x 10 -4 sec -1 
(S 2 Y 1 " 8 1 Y 2 ) = 10 " 2 ( 10 - 2 )(°« 7 ) = 7 x 10 - ' sec -2 


c . 


Gradiometer — Satellite Assumptions 


Misalignment: a = 8 = 10 -3 rad 

Coning Angle: 0 = 10 - ^ rad 


n 2 = 1.5 x 10 -6 sec -2 


U) g = 22.7 rad/sec 


tji % -2o) e * - 45.4 rad/sec 

k. 


s 

C - B 


A 

B - A 
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Assume k^ = 0.99 and k 3 = 0.01; the result is that 
A = 0.497 C and B = 0.507 C; then k 2 = C - A/B - 0.99. 

2 . Angular Rate Error Evaluation 

Before evaluating the error magnitudes, estimate 
the constants of eqs. 9(a) and 9(b): 



1 - 0.99 
1 - ( 0 . 99) 2 


0.01 

0.02 


(A + B) = (D - C) 


-3O 2 0k 9 (l - k.) 

V 1 “ k l k 2 ) 


(Equality in the special 
case of k^ = k 2 ) 


0.5 


10 


-12 


sec 


-1 


(19) 


and 

(A - B) 


(C + D) 


(k^ - (w g + 2*) ] (3fi 2 k 2 0) 
(k 2 k l^ ‘ ~s + 2 ^ 2 ) 


10 12 sec -1 


( 20 ) 


Another set of constants which must be estimated is and 

C 2 that are connected with the homogeneous solution for 

w and w of eq. (13) . 
x y 

An estimate of the magnitude of both and C 2 comes 
from the standard expression for the magnitude of the 
transverse angular velocity of a torque-free body of 
revolution: * 


* 

Introduction to Space Dynamics by William Tyrell Thomson, 
John Wiley 6 Sons (N.Y.) 1963. 
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c 

w transverse ” A w s *" an 
= 2(25) (10 -5 ) 

-4 

- 5 x 10 rad/sec 

Ci ^ C 2 - 5 x 10 4 rad/sec . (21) 

From Table 1, Fig. 4, and the preceding assumptions 
and estimates, we can construct Table III: 

3 . Angular Acceleration Error Evaluation 

The angular acceleration error has all of the 
frequency components of Table III, but the amplitude of 
these error components is a factor of ^( = 0.01) less than 
the corresponding rotational field errors. In addition, the 
angular acceleration error includes the data of Table IV. 

4 . Mass Unbalance Error Evaluation 

To find the mass unbalance errors it is 
reasonable to use Peterson's simplification (eq. (43) in 
RR 441) , with a term added to take into consideration the 
effect of aerodynamic deceleration. 


r 5 

Umax 


m6£. 


[uk + ok + (ok + Wj)u> k ] 


mfiil . 


— t2w. + w. + u). + 2w? + 2w.w.J 

j K X 3 K X 3 


+ 


26 a x m 


@ Freq. 


u> 


s 


(Last term is eq. (18) of 
this report.) 
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TABLE III 

Rotational Field Errors 





Error Equivalent 

Error Frequency, rad/sec 

Order of Magnitude of 

Multiplier 

_2 


Amplitude 

(Fig. 4) 

aec 2 

E.U. 

dc » 0 

10" 6 (4 x 10 2 ) - 4 x 10' 4 aec* 2 

10“ 4 

4 x 10" 8 

40 

w # - 22.7 

louo" 3 ) no' 12 ) • 10 * 14 i«c' 2 

10- 7 

1°‘ 21 

10- 12 

2u» # - 45.4 

io" 25 aac -2 

5 x 10* 5 

5 x 10* 30 

5 x 10" 21 

2* = - 90.8 

10' 25 aec* 2 

IO* 6 

* ID’ 31 

lO' 22 

/k, k_ u» * 22.5 

10' 3 (2 x 10) (5 x IO -4 ) - IO -5 

10- 7 

ID' 12 

ID' 3 

12 ■ 

-2 





aec 




1 /C I K 2 "« * 45 ‘° 

(5 x IO -4 ) 2 - 2.5 x 10* 7 aac -2 

5 X 10" 5 

10 x IO -12 

10' 2 

(i ♦ /iqiq) w # * 45.2 

XO' 14 aac' 2 

5 x 10" 5 

5 X 10" 21 

5 x 10' 12 

a - y^r 2 ) w # * 0.227 

10*“ aac’ 2 

ID* 4 

10' 2U 

io' 11 

2(U # ♦ i) * - 45.4 1 

ID' 21 aac- 2 1 

5 x 10" 5 

5 X IQ" 28 

5 x : )“ 17 

(u + 2>) * - 68.1 

10(10 -3 ) (IO* 10 ) - IO* 12 aac” 2 

5 x 10' 6 

5 x IO* 18 

5 x 10" 9 

O 

2 (“_ ♦ 2*) * - 136.2 

10 -21 aac -2 

5 x 10' 9 

5 x 10" 29 

5 * 10-“ 

• 

2* + (1 ♦ /k^) * - 45. S 

2* ♦ (1 - /k^) ” 90>6 

10 -“ aac 2 
10- 12 aac- 2 

5 X 10' 5 
ID* 6 

5 X 10' 21 
10' 22 

5 X IO* 12 
10' 13 


4 T36b 


II 

B 

H 

il 












I 


I 


m = Mass of one sensor arm 

6 = Magnitude of sensor arm C.M. offset 

from pivot center 

= Magnitude of sensor C.M. offset from 
K spacecraft C.M. — k axis direction 

H. = Magnitude of sensor C.M. offset from 
** spacecraft C.M. — i-j plane 

a y = Aerodynamic deceleration tangent to 
x orbit 



The dimensionless parameters 


m6 £. 


m6 2. 


and 


should be approximately equal. Estimating their size: 


m 'v* 5 kg 

I 3 * 4.5 x 10 6 g cm 2 
6 ^ 10” 2 cm 

^ Jl. . ^ 1.0 cm 

% m ^ij (5 x 10 3 ) (10“ 2 ) (1) % 10 -5 

I 3 I l 4.5 x 10 6 

As indicated by the above I*. expression, there will be 
the full complement of rotational field and angular accelera- 
tion error frequencies but at amplitudes a factor of 
-5 

2 x 10 times the Fig. 4 mass unbalance multipliers. 
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The error caused by aerodynamic drag has magnitude: 


2 x (10~ 2 ) (5 x 10?) 
4.5 x 10 6 


(a x ) (5 x 10“ 2 ) ^ 5 x 10" 9 sec' 2 


= 5 E.U. 


@ U) 


where 


(5 x 10' 2 ) 


- Mass unbalance multiplier at w 


(2 x 10~ 13 )(7.4 x 10 5 ) 2 (3.04 x 10 3 ) 

6.85 x 10 4 

-3 2 

5 x 10 cm/sec ; 

2 

also, the 2w, term in the r expression is dc and of 

k ymax 

magnitude: 

10' 5 ( 2u 2 ) ( 5 x 10' 2 ) % 5 x 10' 4 sec' 2 
s 

= 5 x 10 5 E.U. 

The remaining mass unbalance errors are in the terms: 

2 2 

(u. + oj.), (w. + w.), and (w. + w.)w. . Rather than labori- 

X ] X ] x 3 K 

ously expanding the latter, it is easier to make certain 

observations. The magnitude of any error coming from 
2 2 

(uk + (i) j ) can be taken as roughly equal to the mass unbalance 
errors from 2u>^Uj which were explained previously. Also: 

“k^i + u j* * w s^ w x + w y + " »)« a l 

and 

(w i + d»^) ■ + (0 - a) d» z 


and 


1 (,v2ac d 


M 


s/c 
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The (w + u )w and (<I> + w ) parts each have amplitudes of 

x y 12 * y 

the order of 10 due to the A, B, C, and D coefficients. 

-17 

The error caused by these terms will thus be of order 10 

—2 —16 —2 
sec for (<L + (L ) and 10 sec for u> (w + w ), times 

x y s x y 

the Fig. 4 mass unbalance multipliers (all less than 1.0). 

2 

The remaining errors are (3 - a) w and (3 - a)<I> . 

-4 s -3 z 

(3 - a) could be as small as 10 when a ^ 3 = 10 . The 

2 -4 

error from (3 - a)u) g would be dc and of order (10 ) 

(4 x 10 2 ) (10“ 5 ) (5 x 10~ 2 ) sec -2 = 2 x 10 _8 sec -2 = 20 E.U. 

The errors from (3 - ot)u> will be (10 4 ) (10 8 ) = 10 8 times 

z 

the amplitudes within <L . The latter amplitudes are listed 

z 

in the second column of Table IV. 


Sunmary of Mass Unbalance Errors 

-2 

(values in sec ) 


Less than (2 x 10 -8 ) x (column two. Table 

111) (from u.ut, terms' 

13 

—9 

Less than 10 x (column two. Table IV) 

(From (3 - o)w terms) 
z 


Aerodynamic drag: 5 x 10 
Other dc errors: 


-9 


@ 0 ) 


2w 2 : 5 x 10" 4 


(3 - a) u> 2 : 2 x 10” 8 


-14 

• Errors from (u + w )u> : 10 

x y s 

• Errors from (<L + <S ) : 10~ 15 

x y 

2 2 

• Errors from («. + u>.) : less than 

(2 x 10" 5 ) x (colunm two Table 111) . 
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III. 


CONCLUSIONS 


We may wow look at this wealth of dynamic outputs 

and eliminate from consideration all these with magnitudes 

-11 -2 

less than 0.01 E.U. (10 sec ) because these fall below 
the anticipated threshold noise level of the instrument. 

The remaining outputs are listed in the following 
Table V. 


TABLE V 

Significant Gradiometer Errors Due to Sensor Satellite 

Dynamics 


Error Ouputs 

Frequency 

Amplitude, 
Equivalent E.U. 

Rotational Field 

dc 

40 

During Function 

2/lc.k- ui_ 
12 S 


3 

3 

10 2 

Sum Mode Mismatch 

dc 

0.4 

(w k ) 

/k,k~ u> 
12 s 

o 

1 

I—* 

Mass Unbalance 1 

dc 

5 x 10 5 

‘W 

s 





dc 

20 


TS68 


We observe that dc torques between the sensor arms 
do not generate any steady state signals because of the finite 
impedance across the transducer. In addition# those terms 
at or near the spin frequency will be cut by a factor of 
greater than 100 by the integrating circuits of the 
electronics. 
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I 


This leaves only one term; tc m 2 of the rotational 

field driving functior, which is on the threshold noise level 

and is close to the sensor detection frequency. If we trace 

this term back to its origins, we find that it is strongly 

dependent on 6 (the coning angle) . Any increase in this 

angle above the assumed 10 ^ radians will cause signals 
_2 

above 10 E.U., which fluctuate as the square of this 

coning angle. 

Therefore, it must be determined what control can be 
placed on 8 to maintain it below 10 ^ radians while data are 
being taken. 
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APPENDIX: DIMENSIONLESS FREQUENCY RESPONSE FUNCTIONS AS 

PLOTTED IN FIG. 4 


Mass Unbalance 


(S 2 (2) + S(a x + a 2 ) I 8 2 ♦ B 2 ) ] -- 


(S 2 + a 12 S + B 2 2 ) 


. 58 .. 


"0) 


[(3 2 + 3 2 ; - 2uj 2 + ] W (a 1 + * 2 ) 1 ^ 

2 2 
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APPENDIX D 
AAFE ARM DESIGN 


Several requirements are imposed on the design of the arms of 
the rotating gravity gradiometer (Fig. D-l). Of primary importance is 
the fact that the structure must be isoelastic, i. e. , the bending and 
longitudinal stiffnesses of the arm must be equal, and as large as 
practical. Also, the rotational stiffness (tangential motion of the end 
masses while the center is held fixed), as well as the flapping stiffness, 
i. e. , rotation about an axis normal to the spin-axis, should be as large 
as practical. 

For the AAFE sensor, it is desirable to have the ratio of arm 
width to length as small as possible to keep weight down and moment 
of inertia efficiency up, and to allow for electronic packaging space 
between the arms. 

Several arm design configurations have been considered; the 
one described and analyzed herein, however, appears to best meet all 
the foregoing requirements. 

A. LOAD ANALYSIS 

1. Bending Load 


Pg is the ma acceleration force acting on the arm end mass. 
The tensile force, R, in the upper arm strut causes a compressive 
force, F, in the upper support web. From the defined geometry, it 
can be seen that 


R 


B 



B 
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Fig. D-l(c). :isor Arm Schematic. 
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Longitudinal Load 


tors- 1* 



Fig. D-2. Longitudinal Loading Diagram. 


Tor a longitudinal acceleration load, the cross support will 
carry no load. Each support web will bend to react with P^/2 induced 
force. The tension in the arm strut is: 
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Fig. D-3. 

Arm Strut Tension 
Di agram. 


B. DEFLECTION ANALYSIS 

For a bending load, we may analyze only half the structure, 
since it is assumed symmetrical. 
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Fig. D-4 . 
Schematl c. 
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It will be assumed that the horizontal resistive force provided by 
bending of the support webs is negligible. Hence, the deflection picture 
may be composed of two portions: axial deflection of the cross supports 

plus axial deflection of the arm struts. The cross supports will con- 
tribute to the end vertical deflection as follows: 
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Fig. D-5. 
Bending Load 
Di agram. 


To compute the end vertical deflection caused by axial deflection 
of the arm struts, note that the locus of the end point deflection is the 
intersection of two arcs of a circle swung from the point of attachment 
of the arm struts to the cross supports: 
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Fig. D-6. 

Deflecti on Di agram. 



Simultaneous solution of the equations of these circles yields 


\ 



Note that the x component of this deflection is second order 
and therefore negligible. 

The total vertical deflection of the arm end is therefore 




' 6 1 + E 2 




Aj and A^ can be computed in terms of the vertical load, P, 
as follows. The axial deflection of a member is A = P£/AE where 


P = axial ^orce 
£ = length of member 
A = cross section area 
E = Young's modulus . 
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For a tapered member of the following dimension 
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Fig. D-7. Tapered Member Dimen- 
sion. 


the axial deflection is given by: 


PI , ^2 A PI 

« d 2 - d i'“ d i = A e££ E 



Hence, A, is given by 

T 

a & bs ■-> 

A 1 = A^E = 2h A" E P B 

and 

R l l 2 

A 2 ~ A^E = '2h AjE P B 
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The total vertical deflection 6g is then: 


6 b = 5 A i + I *2 


( 


u 2 .3 

— + — L | p 

2 . „ ,, 2 . „ I B 


2h AjE 2h A 2 E 


6 - P B 

( jel \ 

B ' 2h 2 E 

l A 1 A 2 / 


For longitudinal end loads, the deflections of the arm struts and 
the deflection of the support web comprise the total deflection. Because 
of symmetry, only half the structure need be analyzed. 
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Fig. D-8. 
Total Deflec- 
tion Diagram. 


The total end deflection, 6^, is composed of the support web 
bending deflection, Aj plus the arm struts axial deflection, 
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The bending deflection, 6^ is: 

I 

I 

. . 

! 6 1 " 6EI 


The deflection A^ may be derived as follows: 



By using the Pythagorean theorem 

(t + A) 2 = (b + 6 2 ) 2 + h 2 

and 




e 


2 


b 2 + h 2 




subtracting, and ignoring the second order terms, 


S 2 * I* 


I 


17 ? 
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A is now given by 


and 6^ is 


A = 


Ti 

A 2 E 


p L * 
2h A 2 E 


6 , = t- 


P L r 

2HAJE 


p l £ 


2bh A 2 E 


6 l + 6 2 


6EI 


p l £ 

2bh A £ E 


or 


P L I 3t 3 
l. “ I bhA, 



where I is the section moment of inertia of the support web. 
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CROSS-SECTION PROPERTIES 


2073 - 2 * 
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Fig. 0-10. Tapered Membe 



D. 


NUMERICAL EVALUATION 


Assume the following arm dimensions: 


a = 

5 in. 

f = 

to be determined 

b = 

10 in. 

h = 

3 

c = 

1. 25 in. 

e = 

>/h 2 + b 2 = 11. 18 

d i = 

0. 5 in. 

u = 

0. 5 in. 

= 

1.5 in. 




The width of the support web, f, will be determined by making 
6g = 6^ for Pg = P^, i. e. , an isoelastic configuration. The cross- 
section properties are: 


. 2 


Aj = w d^ = 0. 5 x 1. 5 = 0. 75 in. 


(1. 5 - 0. 5) x 0. 5 0. 5 

* * V “ v* 


'eff 


In 3 


1 . 1 


455 in. 


d 2 r - , . 

i = = j r f = °* 125 f in - in ln - > 


The end bending deflection ir-. 


B 

*B/E 


1 /ab 2 A f 3 \ 1 / 500 . 11. )8 3 \ 

17 \^T + X7/ s TI y oTTs ITffr J 


-1 


0. 0555 (657 + 3, 680) = 38 + 204 = 242 in. 

J •— relative deflection of 
/ arm strut or 84. 5% 

^ relative deflection of 

cross support or 15. 5% 
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The end longitudinal deflection is: 


P l /E " 6 \ bh A 2 


t)'S( 


3x11.18 
10 x 3 x 0. 455 


TTs' \ 

. 125 f / 


0. 1667 I 307 + 


15. 625 


= 51 + 


2. 605 


Equating the two normalized deflections 


242 = 51 + 


2. 605 


f3 = 2 ^ 5 ! = °- 01362 


f = 0. 239 in. 


°L 

p^7e 


= 51 + 191 = 242 




relative deflection of support 
web or 79% 


relative deflection of arm strut 
or 21%. 


I 



Absolute spring ratio for aluminum arm: 


K 



p l /e 


x E 


10 x 10 6 
242 


41, 300 lb/in. 


f 


n 


JL\/li - 1 a/4 1, 300 x 454 x 386 

2ir V m ” 2ir V 2000 


303 Hz . 


E. ROTATIONAL STIFFNESS 


tors-27 



Fig. D- 11. Moment Loading Bending Diagram. 


For the indicated loading, the cross supports carry no load. 
Bending deflection of the support webs and axial deflection of the arm 
struts account for the total arm end point deflection, 6^. 

1. Loads Analysis 

The loads induced in the structure are similar to those pre- 
viously derived. The tension (or compression) in the arm struts is 

R - p 

K R 2h *R * 
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The force acting to bend the support webs is 


F 


R 


b 

2h 


P 


R 


2. Deflection Analysis 


The total deflection of the arm end point 


» & R . is 


R 


= 6 1 + 6 2 


where 

= the vertical deflection of the arm end point caused by 
bending deflection of the support webs 

6^ = the vertical deflection of the arm end point caused by 
axial deflection of the arm struts. 


2073 * 2 % 
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Fi 9 . D- 12. 
Bendi ng 
Diagram. 


6 


1 


is computed as 



where 


or 


f r c ' 


b c- 


3 El 6hEI R 


6h 2 El R 


>2 is computed as before, i. e. 


6 = — A 

°2 h 2 


where 


Rf £ 

A^ = axial deflection of arm struts = ^ — g = ^hA E F R 


whence 


6 2 ■ 7 


2h E 


The total deflection is 


6 R = 6i + 6 2 


/b 2 c 3 + £ 3 \ 

^-5T- + 7K- Z ) 


2 / Eh 
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For the arm design parameters (see page 180), this rotational 
stiffness can be evaluated as follows: 


I = 0. 125 f = 0.125 x 0.01362 = 0. 001703 in 


. 4 


A = 0. 455 in. 2 
eff 


'R 


R/E 


1 / b 2 c 3 £ 3 \ 1 

hM^T- 2A 2 /=9 


1 / 100 x 1. 25 11. 18' 

9 l 6 x 0. 001703 2 x 0. 455 


= | 1 19, 130 + 1 535 j 


130 + 1535 = 2130 + 171 = 2301 


L 


relative deflection of 
arm strut or 7. 4% 

relative deflection of support 
web or 92, 6%. 


Absolute spring rate (for aluminum): 


K R = ^ * E s iS mr- -- 4350 lb/i “- 


Natural rotation mode frequency: 


i = —\f— = -i- \/ 

R 2ir V m 2, V 


4350 x 454 x 386 
2000 


= 98 Hz 


F. MASS AND MOMENT OF INERTIA ANALYSIS 


Weights and mass momenta oi inertia of all component parts of 

the struefure are computed, assuming the arm material is aluminum 

3 

with a weight density of 0. 1 lb/in. . 
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1. Weight Estimate 

a. Arm Str ut 

1 1 5 + C 5 

W = pV = 0.1x9x|x 1> ^ - 5 = 0.45 1b 

b. Cross Support 

W = pV = 0. 1 X 10.25 x^x 1.5 = 0.77 1b 

c. Support Web 

W = pV = 0. 1 x 2.25 x 0.430 x 1. 5 = 0. 145 1b 

d. Central Web 

1. Ring W = 0. 1 x^4. 25 2 - 3 2 J x 1.5 = 1.07 1b 

2. Webs W = 2 x 0. 1 2. 75 x 1. 5 = 1.24 1b 

Total center web 2. 31 lb 

e. End Masses 

W = * 2 (includes aluminum ring) = 8. 80 lb 



1 

1 


Sumn. i'v 

(lb) 

Arm Struts 

1. 80 

Cross Supports 

1. 54 

Support Webs 

0. 58 

Central Web 

2. 31 

End Masses 

8. 80 

Total Arm Weight 

15. 03 lb 


2. Moment of Inertia Estimates 


a. Arm Strut 

Assume it is a line mass with linearly varying weight 
distribution over its length: 



20 T 3 -t» 




Fig. D-13. 

Mass Distribu- 
tion Diagram. 
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A - \ 2 

v 4 = taper ratio = -j 


dm = P A o [l - £ x] 1 dx 


XX 


r +2h/3 

■if 

y -h/3 


^ * 2 p^ 1 ' h x ) 


dx = 


h 4 £ 

0. 096 pA -2-^ 
b 


4 4 

= 0. 096 X m^ = 0. 096 x|x 0. 45 x y = 0. 0525 lb-in. 2 
£ b d £ \o £ 


I = I + Wy 

XX XX 7 

o 


= 0. 0525 + 0. 45 x(|x 3 J = lJJ52JLb-in. 




I yy = j x 2 dm = ||\ 2 = 0. 0525 = °- 583 lb-in. 2 


I = I + Wx 2 = 0. 583 + 0. 45 x (5. 125 + -|) = 17. 46 lb-in. 2 

yy yy rt \ 3 / 


Because of the "line mass" assumption, I = I + I , 

r zz xx yy 


I = 1.852 + 17.46 = 19. 31 lb-in. 2 

zz 
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b. Cross Support 

Assume a "line mass" distribution: 


I = 0 
xx 

o 


I = mr I 2 = 0. 77 x 3 2 - 6. 93 lb-in. 2 
xx 


T m r 0. 77 . 2 

I = — t-t- = — v x 10. 25 - 6. 75 lb-in. 

yy 12 12 ===== 


I = I + I = 13. 68 lb-in. 

zz XX yy . . _ , — 


c. Support Web 


‘xx 12 


mlZ + mr 2 = 0. 145 x + 1. 625‘ 


12 


— > 


I s mr 2 = 0.145 x 5 2 = 3. 62 lb-in. 2 

yy ______ 


I s 0.44 + 3.62 = 4. 06 lb-in. 2 
zz 


d. Center Web 


1. Rins 


2 , 2 
r l + r 2 

I = I = m - 1 — -r — 
xx yy 4 


- ( 2 . 125 2 + 1. 5 2 ) = 


I =21 = 3. 62 lb-in. 

ZZ XX - 


. 44 lb- in. 2 


1.81 lb- in. 2 
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Webs (Both) 


£ 2 , ,4 2 2 

m Yz UZ4 x Tz 


0. 41 lb-in. 2 


St 2 

I = m -pr + mr = 1 

yy 12 


, 4 (H 


+ 3. 3 2 I = 14. 28 lb-in. 2 


I = 0.41 + 14. 28 = 14. 69 lb-in. 
zz — . 


End Mi. .sea 


1. Aluminum Arm End 


= •>5 K ! - d 1 2 )« = 


3 2 - 2 2 \ x 1. 0 x 0. 1 


= 0. 392 lb 


, 2 . 2 , 

m r l + r 2 , £ 2 

4 + 12 


0.3, 2 UW^ 


0. 351 lb-in. 2 


I = I + m£ 2 = 0.351 4 0. 392 x 15 2 = 88. 55 lb-in. 2 
yy xx ■ a - 


I = m tLllL t . 0 . 392 1.5 2 + l.q 2 

zz 2 2 


4 0. 392 x IT 2 = 88. 84 lb-in. 2 


W = 4. 4 - 0. 392 = 4. 008 lb 


1= -fit = — 4x ^0 08 2>l25 . n> 


ptrd 0. 6 x ir x 2 
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XX 

-G*£) 
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= 4. 008 + 

2. 125“ 
12 

I 

yy 

I + m2 2 = 

XX 

2. 51 + 4. 008 

X 15^ 
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zz 
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00 
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Moment of Inertia Summary 



= 2. 51 lb-in. 2 


904. 31 lb-in. 2 


903. 80 lb-in. 2 


TABLE D-I 

Moment of Inertia Summary 


Item 

Quantity 

I 

XX 

I 

yy 

I 

zz 

I 

XX 

I 

yy 

Arm Strut 

4 

1. 85 

17. 46 

19. 31 

7.40 

69. 80 

Cross Support 

2 

6. 93 

6. 75 

13. 68 

U> 

00 

O' 

13. 50 

Support Web 

4 

0.44 

3. 62 

4. 06 

1.76 

14. 50 

Center Web 
Ring 

1 

1 . b . 

1. 81 

3. 62 

1. 81 

1. 81 

Center Webs 

1 

0.41 

14. 28 

14. 69 

0.41 

14. 28 

Alum. Arm 
End 

2 

0. 35 

88. 55 

88. 84 

0.70 

177. 10 

Tungsten 

Masses 

2 

2. 51 

904. 31 

903. 80 

i 

5. 02 

1808. 62 


Totals for one arm 


30. 96 

2099. 61 


Inertia efficiency ratio: 


Srv ‘ X xx 2099. 61 - 30. 96 
^ = i — — 2T2ir?9 — 

zz 
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DETERMINATION OF A NON-ECLIPSE ORBIT 
FOR GRADIOMETEK SATELLITE 



PACT, BUW NOT SUM® 

APPENDIX E 

DETERMINATION OF A NON-ECLIPSE 
ORBIT FOR GRADIOMETER SATELLITE 

A. SUMMARY 

It appears possible to launch the gradiometer satellite into a 
polar orbit of the earth, which will not cause the satellite to be eclipsed 
by the earth even though only a relatively low orbit has been z chieved. 
Depending on the orbit altitude, it is possible to achieve non-eclipse 
periods several times or more longer than the approximately one-week 
period necessary for complete surface coverage. The advantages A the 
non-eclipse orbit are the weight reduction and reliability increase asso- 
ciated with no dependence on chemical batteries. Also, because the 
gravity gradient sensor output is thermally sensitive, the required 
thermal control system would, therefore, have only one state of thermal 
equilibrium to contend with rather than cycling between two. The trade- 
off does not seem too severe. A launch window is created which 
previously was not a mission constraint. 

Figure E-l shows the earth at the time of the vernal equinox 
(around March 21). Th* r i-:„» a^en looking down from the north 
ecliptic pole (NEP). The -* t i L .nd sovth celestial poles (NCP and SCP) 
are each ~23-l/2° away lr - -I.e a- ant ecliptic poles and at the 
times of vernal and autumnal a l ie on the terminator . 

The terminator rotate.- c ■». rclockwise about the polar axis to 
the ecliptic (passing through NJu.-'i a a the earth orbits the sun. The 
terminator rotates with the earth's orbital angular velocity, which is 
about l°/day. Considering the time around the vernal equinox: If the 
satellite is placed in a polar orbit plane that coincides with the termina- 
tor plane at the time of the vernal equinox , then no eclipse of the 
satellite will occur at the time of the vernal equinox. If the satellite 
is launched into the latter orbit some days before t* vernal equinox, 
the orb u plane and terminator plane intersect at an angle, 0 , whose 
vertex is at the NEP (Fig. E-2). 
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Fig. E-l. Earth at Vernal Equinox. 


tors - it 



Fig. E-2. Earth Position Several Days Be 
fore Vernal Equinox. 
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The main question becomes: How large can 0 be for a given 

altitude of the specified orbit with the satellite not eclipsed as it travels 
over the dark area of the earth? An answer to this question gives the 
non-eclipse period (maximum), which is centered at each equinox and 
is oi duration: 


Time of non-eclipse 


29° 

1°/DAY 


B. DETERMINING MAXIMUM 9 

One method for finding 0 is to p/cect the circular orbit 
onto the terminator plane and determine i the projection anywhere 
dips within the circular intersection the earth and terminator plan*;. 
This formulatior. tacitly arsumes that the shadow of the earth is cylin- 
drical when in reality the umbra 1 shadow is cornea' with a cone half- 
angU of about 0. 5J°. However, since the cone half-angle is so small, 
the simple approximation of a cylindrical hadow will be used: 

Let Xj, yj be the rectangular orbital coordinates of the satel- 
lite in the orbital plane (Fig. E-3). Then 


tors* st 



Fig. E-3. 
Orbit Plane. 
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Earth Radius = 6371 km 


Xj = (R q + h) cos <t> 


(R + h) sin 
o 


h = Orbit Altitude 


Let the terminator plane projections of these coordinates 
be x, y. Then 


x = x, cos 8 = (R + h) cos «t> cos 8 
1 o 


y = Yj = (R Q + h) sin <f> 


The condition for non-eclipse is: 


VJ77 ’> 


which results in: 


[K + h ) : 


[R^ + h) 2 cos 2 $ cos 2 8 + |r q + hj 2 sin 2 4»j > R c 


Squaring to get 


2 2 2 

cos 4* cos 6 + sin $ > 


K + h )' 


Then 


1 + coa 2 $ (cos 2 0 - 1) > 


1 


R o R 
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2 2 2 2 
Since h /R q << 2 h/R o> neglect h /R q and make approximation: 


1 + cos^ $ (cos^ 0 - 1) ^ 1 - 

o 


The res”lt is: 


cos 4> < 


2h 


R (1 - cos** 0 ) 

o 


(non-eclipse condition) 


Since 0 S cos <$> < 1, the non-eclipse condition must satisfy the most 
severe constraint: 


Hence, 


2h 


R o (1 - cos* 0) 


> 1 


2 a ^ , 2h 
cos 0 > 1 - — — 

o 


cos 6 





(non- eclipse condition) 


The duration vs. altitude of non-eclipse orbits is shown in 
Table E-l. 
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TABLE E-l 


Duration vs Altitude of Non-eclipse Orbits 


h (km) 

h/R 

o 

(1 - h/R o ) 

0° 

max 

Duration =20 /1°/Day (Days) 

250 

0. 0392 

0. 9608 

16° 6' 

32 

275 

0. 0431 

0. 9569 

16° 53’ 

34 

300 

0. 0471 

0. 9529 

17° 39’ 

35 Launch Window is 

Equal to Duration Minus 

325 

0. 0510 

0. 9490 

18° 22’ 

37 Minimum Required 

Mission Time 

350 

0. 0550 

0. 9450 

19 5 

38 

375 

0. 0589 

0.9411 

sO 

0 

<T_ 

40 Launch Windows Occur 

Twice per Year 

400 

0. 0629 

0. 9371 

20 26 

41 


C. SOLAR CELL POWER GENERATION IN NON-ECLIPSE ORBIT 

Since the satellite will have its spin axis perpendicular to the 
non-eclipse orbit plane, the bulk of the power will have to come from 
solar cells mounted on the sun-facing end of the cylinder (Fig. E-4). 

The tentative 76 cm configuration to fit the internal diameter of 

2 -3 

Scout has a 4550 cm area for the cylinder end. At 61. 5 x 10 watts 

2 

per cell, each of which has 4 cm area, the power available from the 
end solar cells is: 


P END = 
0 = 0 ° 


4550 cm 2 x 61. 5 x 10~ 3 watts 
4 cm 2 /cell ccl1 


70 watts 


For the 0 excursions to be expected over mission lifetime 
resulting from terminator rotation, the latter power will decrease by 
(1 - cos 0) or from 4 to 7% for 250- to '50-km orbits. 
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By contrast, were the current configuration to go into an 
eclipsing orbit with the sun vector perpendicular to the cylindrical 
wall, a substantial power available reduction would result. The rectan- 
gular cross section = 76(40) = 3040 cm^ leading to: 

3040 . _ 

P = Veen x /0 = 47 watts . 

4550 

The 47 watts would hardly be sufficient to operate the space- 
craft and recharge batteries. But 70 watts should be sufficient for 
operation in a non-eclipse orbit. 
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APPENDIX F 

RGG TORQUES RESULTING FROM PARAMAGNETIC ARMS 

A. INTRODUCTION 

This section is based largely on Ref. F-l. 

In the cgs electromagnetic system, the force on a magnetic 

pole is 


F = Hm 


(F-l) 


where 


F = force on pole in dynes 
H = field strength in oersteds 
m = pole strength in emu. 

The flux emanating from a magnetic pole is 

4> = 4irm maxwells (F-2) 

If a 8 oft ferromagnetic or paramagnetic bar is introduced into a uniform 
field, the bar will align itself with that field, if the bar is unrestrained. 
Magnetic poles will be induced in the bar. If it is assumed that all the 
flux from the induced poles emanates from the end of the bar (an over- 
simplification, but it does provide simple equations with practically 
acceptable accuracies for relatively long bars), the flux balance will be 

4>. - 4> = 4irm (F-3) 

D O 


205 



where 


4>, = flux in the bar 

b 


4> q = flux in the volume occupied by the bar with 
the bar removed 


47rm = flux caused by the induced poles in the bar* 


Dividing Eq. (F-3) by the cross-section area of the bar, A: 


^b ^o 4?rm 
A " A = A 


But, /A is defined 
b 


and 


^/A = = flux density in the bar 


$ q /A = B q = flux density in free space 


(F-4) 


I 

I 

I 

I 

3 

] 

3 


and 


B q s H in the emu system. 


Rewriting Eq. (F-4), 


n t t 4rm 

B b • H - nr 


11 

11 

ii 

ii 

« 

ij 
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and 


B 


b 


H + 


4nrn 

A 


B. PERMEABILITY AND SUSCEPTIBILITY 

Dividing Eq. (F-5) by H 

B . 4Trm 

H “ 1 HA 


(F-5) 




II _ . 4trm 
H " HA 


(F- 6 ) 


p is defined in the emu system by Eq. (F- 6 ). 777 - = K = the magnetic 

HA 

susceptibility of the substance of the bar. It is common in scientific 
work to define the susceptibility per unit mass X , 


X = Kp (F-7) 

where p is the density of the substance. X is also frequently based on 
either the atomic weight X & or on the molecular weight X . In this 
discussion, X is based on density exclusively. Thus, 

p = 1 + 4irK 

= 1 +■ 4*Xp . (F-8) 


The permeability of some materials is shown in Table (F-I) 
from Ref. F-l, F-3, and F-5. 
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TABLE F-I 


Magnetic Properties of Pure Materials 



X £> 
x 10 6 

p 

H- 

Aluminum 

+0. 63 

2. 68 

1. 0000212 

Beryllium 

-1. 0 

1. 85 

0. 9999767 

Copper 

-0. 086 

8. 96 

0. 9999900 

Tungsten 

+0. 28 

19. 3 

1. 0000675 

Uranium 

+2. 6 

19. 1 

1. 000624 

Elgiloy 



1. 00005 

Air 



1. 00000036 

Permalloy 



100, 000 

Iron, technically pure 



6, 500 

Cast iron 



600 


C. TORQUE ON A BAR MAGNET 

If a bar magnet of length, i , with pole strengths, tm and -m, 
is introduced into a uniform field of strength, H, at an angle, 6 , to 
the field, the following is derived from Eq. F - 1 and the geometry 


Torque = mlH sin 6 . 


(F-9) 




I 


D. TORQUE ON A PARAMAGNETIC BAR 
It is seen from Eq. (F-3) 

4> b - 4 > q = Aitm (F-3) 


that a pole of strenth, m, is induced in a bar when it is introduced into 
a field. From Eq. (F-6) for a field of H oersteds 


M- 


1 + 


4nm 

HA 


(F-6) 


or 




. 47rml ! _ . 4ffm i 

HAT ' H Vol 


( F — 1 0) 


or 


H (p - I) 


4wml 

Vol 


and 


mi 


(u - 1) H Vol 
4 it 


(F- 1 1) 


This "magnet" of pole strength m and length t can be substituted into 
Eq. (F-9) and 


Torque = ll» - ») «| Vol . qiu 6 , . (F-121 


209 



Two elements have not been taken into account the derivation of 
Eq. (F-12): the demagnetizing effect and th?. : eduv < 1 projected area 
of the ends of the bar when the latter is not parallel the field. The 
reduced area coefficient is simply cos 0. 

When the bar is introduced into the field, the induced poles tend 
to counteract the initial field. Thus, the effective field H is 

H - H’ - NKH = H' - NXpH (F-13) 


where 


H = effective field 

H' = initial field in a vacuum 

N = demagnetizing factor. 

N varies from nearly zero for long thin rods to 4# for flat plates. For 
paramagnetic materials with mass susceptibilities in the ranges shown 
in Table F-I, 


H » H' 


even if N is as large as 4», its maximum value. 

Thus, the final torque equation for a paramagnetic bar introduced 
into a magnetic field of strength H is 


or 


Torque = U- H ---Vpl * in JLSSl Q 


(F-14) 


Torque = XpH 2 Vol sin 0 cos 6 


(F-15) 
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where 


Torque = dyn cm 


X - emu units 


P = g/cc 

H = field strength in oersteds 


E. FORCES ON A PARAMAGNETIC MASS IN NON-UNIFORM FIELD 


The force on a homogeneous mass, M, in a non-uniform field 
is derived by Williams (Ref. F-l). 


F = 

F = 
X = 
M = 
H = 
dH/dx = 


dH 

X MH AX 


dynes 

emn 

mass 

oersteds 

oersteds/cm . 


F. 


MAGNETIC TORQUES ON THE AAFE RGG 


Torques can be introduced into the sensor either by the unequal 
attraction of the end masses or by the magnetic poles introduced in the 
paramagnetic arms. Because of the symmetry of the end masses, poles 
induced in them will not introduce torques. 

The orbit of the AAFE sensor is such that the spin axis of the 
sensor is at all times orthogonal to the earth's spin axis. That is, 
the plane of the sensor arms lies in a longitudinal plane through the 
earth. 
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HRL notes state that "the earth can be considered «.s a dipole 

25 

magnet with a. pole strength of aoout 8. 1 x 10 emu. (The field 

gradient in a longitudinal plane) varies from about 0. 005 7/ft altitude 
at the equator to 0. 01 7/ ft altitude at the poles. " 

Taking the worst case at the pole 


= 10 7 oersted/ft (1 Y = 10’ 5 oersted) 

= 3.3 x lO -9 oersted/cm . (F-17) 

Since the gradient cnanges only by a factor of 2 to 1 as the sensor 
revolves around the earth, the maximum torque resulting from the 
field gradient occurs when the arms are at an angle, 6, to the field 
and the field differs because of the separation of the masses. The 
torque on one arm is 

Torque = fjYj - {^2 = ^1 ’ t ?) Y co * 0 

= - XM,H 2 |g) Y coaO 


but 



H. + Y 


dH 

sin 6 


and 


Torque 


X.M, 


dH 
dX v 


cos 8 


( H i • H . 


y sin 8 


Torque 


XM, v Z 


cos 8 sin 8 



(F- 18) 
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Even though and V are large for this sensor, y^x) * s SO smai ^ 
the torque caused by this effect is insignificant even in this sensor. 
The induced pole alignment torque is from Eq. (F-15) 


Torque 


XpH 2 


Vol sin 6 cos 6 


(F- 1 5) 


Assume 

X = 0. 63 x 10 ^ for alumin’im 
P = 2. 68 g /cc for aluminum 
H = 0. 5 oersted 

Vol = 60 cm long x 7. 5 cm thick x 5 cm wide 
= 2250 cm 3 
6 = 45° 


Torque = 0. 63 x 10' 6 x 2. 68 x 2. 25 x 10 3 x 0. 5 2 x 0. 707 2 
Torque = 4, 8 x 10 ^ dyn cm (F-19) 


because of paramagnetic aluminum arms. 

The AAFE sensor has a polar moment of inertia for one arm 
2 

of I = 0. 622 N msec and an inertia efficiency of about 97%. Thus the 

2 

effective inertia, I , equals 0. 6 N msec . 


Torque resulting from gradient = 

Torque = 0.6x3/2x0.01x10"^ 
- 5 * 

Torque = 9 x 10 dyn cm 
because of 0. 01 EU gradient input. 



= 9 x 


I 3/2 r 


e eq 

- 12 

10 N m 


(F-20) 
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Thus, the peak torque resulting from paramagnetic aluminum 
arms is 5. 34 times the peak torque caused by the desired gravity 
gradient resolution. This torque will vary by a factor of about 2 to 1 
while rotating around the earth. 

Table F-I shows that it makes little difference whether aluminum 
or beryllium it used. The Xp product is essentially the same for both 
of them, but the polarity of the torque would be reversed. 

Since the field strength will vary by about 2 to 1 as the sensor 
orbits the earth, it seems obvious that the sensor should have a magnetic 
shield. This is discussed in Section F-H of this appendix. 

The altitude of the orbit, 250 to 300 km, puts the sensor 
generally in the Van Allen belt. Variations in this belt will also cause 
significant variations in H; therefore, magnetic shielding is again 
required. 

G. EFFECTS OF FERROMAGNETIC IMPURITIES IN SENSOR ARMS 

The effects of ferromagnetic impurities are exceedingly difficult 
to estimate. Iron seems to be an impurity in almost all commercial 
and even technically pure alloys. However, a number of compounds 
of iron and other ferromagnetic materials are only slightly paramagnetic. 
If the impurity is not in a combined form, the effect is particle size, 
heat treat, and cold work dependent. 

The problem may be bounded by assuming that the iron inclusion 
has the permeability of cast iron; it is evenly dispersed and is 1% by 
weight. 


p cast iron * 600 = 1 + 4irXp 
(Xp) iron = 599/4ir = 47.7 
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If the iron is 1% by weight, the volume will be only 


Vol iron = Vol a t x 0. 01 x 2. 68/7. 85 
= Vol a i x 0. 0034 

Assuming the susceptibilities add directly, the torque resulting from 
a 1% free iron impurity is 


Torque = ( X ai p a£ + 0.0034 X^p^) H Vol sin 0 cos 0. 


The torque is increased by a factor of approximately 


>-3 


3.4 x 10 x 47. 7 
0. 63 x 10' 


A Torque * ~ = 2. 6 x 10‘ 


The total peak torque will be 

Torque = 4. 8 x 10 2.6x10^ = 124 dyn cm (F-21) 

because of 1% free iron impurity in the AAFE aluminum arms and the 
earth's field. 

H. MAGNETIC SHIELDING OF RGG 

Since the AAFE sensor arms rotate in a longitudinal plane 
through the earth's field, the sensor can be shielded by a cylindrical 
shield with the axis of the cylinder along the spin axis of the sensor. 

The cylinder should be as long as convenient, but the ends need not be 
covered if the shield extends 5 or 6 in. beyond the arms. 
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Reference F-2, Section 5-33, states that the induction inside 
such a cylinder, B^, is uniform and has the following ratio to the 
external field B: 


_B 

B. 

1 


4ub 2 + (u - l) 2 (b 2 - a 2 ) 
4pb 2 


(F -22) 


where 

p = relative permeability of the shield. (This is the 
same as p in emu units given in Table F-I and 
the previous equations. ) 

b = external radius of the shield 

a = internal radius of the shield 

For thin shields of thickness, t, and if p 5>1, as any good shield will 
be, this can be approximated 

^ = 1 + °* ^ - . (F-23) 

i 

For the AAFE sensor, assume t = 0. 01 in. , b = 15 in. and 

4 

p = 10 . (Such high initial permeability will be quite difficult to obtain 
in practice. ) For this condition 


B _ , . 0. 5 x 10 4 x 10" 2 

B - i + f5 

l 


= 4.35 (F -24) 
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The flux density inside the shield has been reduced by a factor of 
4. 35 to 1. 

As was shown in Section F-G, since the induced RGG torque 
varies as the square of H (B and H are essentially equal in the arms 
and free space), this shield should reduce the induced paramagnetic 
torques to nearly acceptable values. However, since high initial 
permeability is difficult to obtain in situ, a shield thickness of 0. 03 in. 
is recommended. This would theoretically provide a shielding factor 
of 10 and a magnetic torque reduction of 100. 

One more calculation is required to ensure that the shield 
material is not saturated by the magnetic flux it is conducting around 
the shielded volume. It can be assumed that all the flux intercepted 
on one -half the cylinder face is concentrated and conducted through the 
shield. 


4> intercepted = H x 2b x l 
Cross section area of shield A. 


At the point of maximum flux density 


A = 2t l 


B shield 


± 


A 


H2bl 

Zth 


0. 5 x 15 

0.01 


750 G . 


(F-25) 


(F -26) 


(F -27) 


This liux density is easily acceptable. 

The shielding ratio for a spherical shell is also given in 
Ref. F-2 as 



atfcj Li (I » - 1> 2 n> 3 - 1 3 ) 

*b 3 


(P-28) 
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which can be approximated for these particular conditions as 

£ = 1 + 2/3 p £ . (F-29) 

i 

For a 0. 03 -in. -thick shell with p = 10, 000, this is 
B 

■g- = 14. 5 shielding ratio for a spherical shell. (F-30) 

i 

Thus a spherical shell provides a little more shielding. 

It probably will be wise to cover the drum heads of the shield 
as a precaution. 

It should be noted that multiple -layered shields provide better 
shielding than a single shield. This layering is significant only if the 
spacing between the layers is a significant fraction of the shield radius. 
In the present case, this is obviously impractical. 
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CONCLUSIONS 


The AAFE RGG sensor should be designed as follows: 


1. The aluminum arms should not have a relative 
permeability of greater than about 1. 0000 50, 
and no ferromagnetic inclusions can be allowed. 

2. The end masses must be symmetrical about an axis 
parallel to the sensor spin axis, i. e. , cylindrical, 
as they are. The permeability can be as high as 

1. 001, but the end masses must be homogeneous. 

No ferromagnetic inclusions or residual magnetism 
can be allowed. 

3. A magnetic shield approximately 0. 03 -in. -thick is 
required. This is visualized as a three -layer wrap 
of 0. 010-in. -thick "Co-Netic” magnetic shielding 
tape. This tape is considered to have the necessary 
high initial permeabilities at the very low flux densi- 
ties involved. This capability is being investigated. 
It should be noted that this tape in the 15-in. width 
weighs 0. 62 lb /linear ft, and will add about 20 lb 

to the sensor. One telephone quote estimates the 
cost at $18/ ft. 

The Rawson flux meter in the HRL laboratory has 
full-scale ranges from 0. 0005 to 5. 0 G. HRL, 
therefore, has no difficulty making tests to ensure 
that adequate shielding has been attained. 
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ABSTRACT 

In this report It Is shown that a simple 
compensation or calibration approach to the re- 
duction of the R6G sensor temperature sensitivity 
appears to be Impractical. Automatic sensor com- 
pensation and adaptive tracking of the sensor 
natural frequency to the satellite rotation speed 
Is proposed. This approacn appears to be feasible. 
Future work and necessary tests are discussed. 
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I. SENSOR DESCRIPTION 

A . General Overview 

The RGG sensor is composed of two isoelastic (in the 
plane perpendicular to the axis of rotation) arms that have 
concentrated masses at their extreme ends. These arms are 
attached to the case by end pivots, also called support 
pivots. The end pivots offer little torsional stiffness 
(they contribute approximately 1 Hz of the 8 Hz sensor 
resonant frequency) but they are quite rigid in the longi- 
tudinal and lateral directions. 

The system main torsional stiffness is provided by 
two piezoelectric bimorph transducers arranged to operate 
in a pure bending mode. The use of the piezoelectric mate- 
rial to provide the main torsional stiffness is desired in 
this sensor, because the signal energy input is so low 
tuat most of it must be transformed by the piezoelectric 
transducers in order to obtain a usable output signal 
level . 

H . Sensor Characteristics Summary (Preliminary) 

The following sensor characteristics are summarized 
largely from the Monthly Progress Reports^ - and from the 


notes 

of C.C. Bell. 


1 . 

RGG 


Arms : 




Material 

Aluminum 


Radius to CG of mass 

0.381 m 


22? 



Mass : 


Material 

Sintered tungsten 

Each end mass 

2.0 kg 

One Complete Arm with Masses: 


Polar moment of inertia 
about spin axis 

2 

0.624 Nm sec 

Mass of one complete arm 

6.8 kg 

Inertia efficiency = A-B/C ^ 

97% 

Sensor Resonant Frequency 

8 IIz 

a) = 50.24 rad, sec 
n 


Sensor Stiffness 


K t = <4(1/2) = (k o + = 788 

dm/rad (6800 in lb/rad) 

2 

K end pivots, each = (6.28) I 

= 24.6 Wm/rad 
(218 in lb/rad) 

K transducers = 775.7 ilm/rad 


K each transducer = 387.8 Nm/rad 

2. Environment 


Temperature 


In pivot region 

28°C ± 0 . 25°C 

In arm region, initial 

24°C + 1°C 

gradually shifting to 
at the end of 30 days 

15°C ± 1°C 

Magnetic Field (approximate) 


Horizontal at equator 

0.6 G 

Vertical at poles 

0.7 G 

Acceleration (any axis) 

< 10" 4 G 

Spin Speed 

4.0 rps 
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STATEMENT OF THE PROBLEM 


II . 


The phase and amplitude of the output of the RGG is 
known to be extremely sensitive to operating temperature. 
This sensitivity is due to the variation in material dimen- 
sions, variation of the modulus of elasticity, and the vari- 
ation of the transducer piezoelectric coefficients with 
temperature. It is shown in Ref. 2, p. 473, eq. (VI-h-15) , 
that the undamped natural frequency of the sensor will 
change : 


Af 

f 


- 5 - 2 


(3a 


pivot 


+ e 


pivot 


- 2a ) 
arms 


AT 


( 1 ) 


where 


a . = temperature coefficient of linear 

P expansion of the pivot material 

a = temperature coefficient of linear 

arm , , . , . . 

expansion of the arm material 

e 0 . . = temperature coefficient of the 

” modulus of elasticity of the pivot 

material 

6 = Af/f , the change in frequency 

divided by the natural frequency. 


The change in the sensor natural frequency due to 
(1) will cause a change in the amplitude of the output as 
shown by (2) , if the rotation speed is maintained constant. 



1/2 


1 + 


2 2 
4Q*V+ 


- 1 


26 


( 2 ) 
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I 


In addition to the amplitude reduction due to operate i off 
the natural frequency of the sensor, the characteristic 
amplitude of the sensor itself changes due to variation in 
moraent of inertia and spring rate. This change, aga n from 
Reference 1, p. 474, is: 



(2a 


arm 


"* a pivot 


e ) AT . 

pivot 


(3) 


The total change in output amplitude is given by the sum of 
(2) and (3) . 

The above equations were derived on the basis of a 
single pivot material; piezoelectric materials and effects 
were neglected, and the inertia of the arm was assumed to 
be due to a concentrated mass at the end. For the proposed 
AAFE RGG, only the last assumption is completely valid. The 
supporting structure of this sensor is composed of metallic 
end pivots, which constitute approximately 1/64 of the total 
stiffness, and two piezoelectric bimorph benders, which pro- 
vide approximately 63/64 of the total. Since no dimensions 
are involved in eqs . (1), (2), (3), or in their derivation, 
the contribution of each material can be proportioned and 
usea directly. (Additional terms would have to be included 
in a final detailed design since the piezoelectric material 
is supported by still another material.) 


A. Sample Temperature Amplitude Sensitivity Calculation 

A sample calculation demonstrates the problem. Assume 
the arms are of aluminum, 

a » 20 x 10 -6 / # C 
arm 
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the support pivots are of beryllium copper, 

a , = 16.6 x 10* 6 /°C 

Pi 

e. . = -350 x 10" 6 /^ 

and the piezoelectric ’pivots" are 4v toTiO^; 5.7% CaTiO^; 
and 90.3% BaTiO^. (Refer to the discussion of piezoelectric 
material characteristics in the following section of this 
report. ) 


a _ = 1 x 10“ 6 /°C 

P2 

e o p 2 = + 400 x 10 6 /°C (approximate) 

It is also assumed that the sensor Q ^ 100, and the tempera- 
ture change is 1 3 C. 

6 =| [3{ (1/64)16.6 + (63/64)1} - (1/64)350 + (63/64)400 
- 2x20] 10" 6 /°C 

6 % 1.75 x 10 -4 /°C . (4) 

Substituting this value of 6 into eq. (2) , we have 
1 = 

4Q 2 5 2 « 4 x 10 4 x 3 x 10 -8 

26 « 2 x 1.75 x 10" 4 

Total denominator of (2) 


1.0 

= 1.2 x 10 

* 0.3 x 10 
- 1.0015 
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and 


\ 


1.0015 


1/2 


1.0000-0.0007 


Thus , 

A0 

= - 0.0007 = - 0 . 07%/°C . (5) 

o 

This is the amplitude change due to sensor frequency change. 

The characteristic amplitude change due to the sensor 
temperature change is found by substituting the same material 
temperature coefficients into eq. (3) , and we find 

= [2x20-3{ (1/64) 16.6 + (63/64)1) 

-(1/64) (-350)-(63/64) (400) ]10~ 6 . (6) 

* - 3.5 x 10 -4 = - 0 . 0 35%/°C 

due to characteristic amplitude changes. The total amplitude 
change for a 1°C increase in temperature is 

- 0.07 (-0.035) = - 0 .105%/°C. (7) 

*o 

For the earth orbiting RGG it is desired to attain an ac- 
curacy of 0.01 E.U. in the presence of the earth gradient 
of approximately 2,250 E.U. or approximately four parts in 
10®. The above estimate shows a variation in sensor output 
amplitude of 1050 parts in 10® for a 1° temperature change. 


A0. 

4 

~e" 


A0, 

4 
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The above calculations were carried out in detail to 
demons trate the magnitude and difficulty of the problem. 

The dominant terms are those involving the temperature coef- 
ficient of the modulus of elasticity. 

D . Temperature Coefficient of Modulus of Elasticity 

For the metals, this coefficient is approximately 
linear, but it is poorly controlled and sensitive to alloying, 
heat treating, and cold working. However, in this sensor the 
effect of the varying modulus of elasticity of the metals is 
completely overshadowed by the effect of the piezoelectric 
materials, since the piezoelectric transducer provides almost 
all of the mechanical stiffness. 

All of the constants of the piezoelectric ceramic 
materials vary with temperature. The dielectric constant, 
modulus of elasticity, and strain coupling coefficients all 
vary with temperature and age. The shape of the curves for 
some typical compositions is shown on pp. 3-104 and 3-105 of 
Reference 3 reproduced here for convenience as Fig. 1. It 
can be seen from these curves that the percent constant 
change per degree C change can be plus, minus, or zero, de- 
pending on the material and operating temperature. 

In the previous example we picked the 4% PbTiO^, 5.7% 
CaTiO-, and 90.3% BaTiO, at 20°C operating temperature and 
found (AE/E)/°C to be approximately + 400 x 10 /°C. I7hen 

this is compared to the curve of 8% PbTiO^* 5.5% CaTiO^, and 
86.6% BaTiO^, it is apparent that we cannot expect to predict 
the modulus of the piezoelectric ceramic or its temperature 
coefficient to better than ± 10%. 
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In view of these uncertainties it seems unlikely that 
a sensor can be designed with a temperature sensitivity un- 
certainty of less than about 100 x 10~ 6 /°C. Thus simple 
compensation or calibration along with closer temperature 
control appears to be impractical at this time. 

III. AH RGG TEMPERATURE COMPENSATION METHOD 

The equivalent circuit of the RGG sensor with two 
piezoelectric bimorph benders and the two sets of end pivots 
is shown in Fig. 2. 


Nj : I ft N,; I 



Fig. 2. Electrical Equivalent of RGG. 
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M 

o 

C p c ^ 
pi p2 


C ' C ' 

ml' m2 


R x , R 2 


C 01' C 02 

a i' a 2 


R 


moment in sensor due to gravity gradient. 

mechanical compliance for arm end pivots. 

mechanical loss in sensor. 

equivalent mechanical compliance of 
piezoelectric bimorph benders. 

mechanical loss in piezoelectric bimorph 
benders . 

capacitance of piezoelectric bimorph benders. 

electromechanical transformation ratio of 
piezoelectric bimorph benders. 

electrically variable capacitance. 

electrically variable resistance. 


In this circuit the mechanical compliance (C^) of the piezo- 
electric benders will vary with temperature as will C ^ and 

C I , M , and cl. and cl _ . Figure 2 can be further sim- 
pz s o 01 02 

plified to Fig. 3, where all parameters e, R, L, and C are 
temperature sensitive. It is desired that the resonant 
frequency of Fig. 3 be equal to twice the spin speed of the 
satellite and the Q maintained constant. By adjusting the 
variable capacitor and the resistor independently these two 
objectives can be attained. 



e out 


Fig. 3. Reduced Equivalent of Fig. 2 . 
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IV. 


FUTURE WORK AND TESTS 


The equivalent circuit of the binorph benders must be 
analyzed in more detail and tests made to insure that the 
equivalent circuit adequately describes the physical system. 
Tests also need to be made with the variable R and C to 
insure that the Q and amplitude can be adjusted over the 
ranges required. 
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GEODESY WITH ORBITING GRAVITY G RADIOMETERS* 


by 

Robert L. Forward 
Exploratory Studies Department 
Hughes Research Laboratories 
Malibu, California 90265 

ABSTRACT 

A new instrument for sensing of the earth's gravity field, 
the rotating gravity gradiometer, has been demonstrated in the 
laboratory. The gravity gradiometer measures the gradient of 
the gravity force field rather them the field itself. The sensor 
does not respond to acceleration and can operate in free fall 
or in accelerating environments where the usual gravity meters 
cannot work. A gradiometer in a spin stabilized satellite in 
a low polar orbit will make a significant contribution to the 
geodetic mapping program presently being carried out by satel- 
lite tracking, since it preferentially senses the higher har- 
monics (> 35) of the earth's field where the doppler tracking 
signals (revised Williams town report) fall off rapidly. 


Preprint of paper to be published in the AGlf Monograph - 
Symposium on the Use of Artificial Satellites in Geodesy, 
19-20 April 1971, Washington, E.C. 
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I. DISCUSSION 


A. Rotating Gravity Gradiometer 

The rotating gravity gradient sensor that has been de- 
veloped, Bell [1970], and operates daily in the Hughes Research 
Laboratories is a device for measurement of the second oruer 
gradient of the total gravity potential field. The sensor con- 
figuration consists of a resonant cruciform mass-spring system 
with a torsional vibrational mode (see Fig. 1) . In operation, 
the sensor is •'orated about its torsionally resonant axis at 
an angular rat*, to which is exactly one-half the torsional reso- 
nant frequency. When a gravitational field is present, the 
differential forces on the senso" resulting from the gradients 
of the gravitational field excite the sensor structure at twice 
the rotation frequency. Forward [1965b]. Only the differential 
torque AT between the sensor arms at the doubled frequency is 
coupled into the sensor output. 

For the simple example shown in Fig. 1, the differential 
torque induced by the mass M at the distance R is 

AT * l — £ - 2 sin 2wt . (1) 

q R J 


The angular resonant deflection between the two quadrupoles 
of the sensor rotating at one-half its torsional resonant 
frequency with an associated quality factor Q is therefore 


e 


2 ATQ _ , GM Q . 

■ ^ * 3 -r * cos 2u>t 9 

I (2u>) * R J (2ta>r 


( 2 ) 
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Fig. 1. Method of Operation of Torsional Gravity 
Gradiometer. 
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2 

where I = ml /2 is the quadrupole inertia. The angle 9 is 
extremely small. Surface gradients produced by the earth will 
produce angular deflections of 9 'v 3 x 10 - ^ rad in typical or- 
bital torsional sensor designs (see Fig. 2) , while useful 
threshold signals of 10 ^ gal/cm (0.01 Eotvos unit (E.U.)) 
produce angular responses of 'v, 10”** rad. 

Although these deflections are small, they are easily 
measured by utilizing piezoelectric strain transducers attached 

to the torsional flexure. The threshold deflections of 10~^' 1 ' rad 

• 8 

produce voltage outputs of 10 V from typical transducers. These 
voltage levels are easily measured by modern amplifiers. 

B. Present Development Status 

The ultimate objective of our work on rotating gravita- 
tional gradient sensors is the development of a class of rugged 
sensors of high sensitivity and precision which may be used to 
measure accurately and rapidly the details of the gravity field 
during airborne or orbital surveys and as a component in an in- 
ertial guidance system to remove the effects of gravitational 
anomalies on the ultimate system performance. 

The objectives of the initial research programs were to 
investigate the engineering feasibility of the basic concept, 
to develop sensor structures which would operate at a high sen- 
sitivity level both in free fall and in 1 G environment, to 
measure the sensor's sensitivity to gravitational fields, and 
to investigate the sources of noise produced by the rotation of 
the sensor. A torsionally flexible structure utilizing 
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piezoelectric readout was found to be a suitable design and 
offers a significant improvement over other possible gradi- 
ometer designs (see extensive bibliography in Bell [1970]). 

It has demonstrated the capability of being operated in an 
earthbound laboratory environment (see Fig. 3) while still 
maintaining a high sensitivity and low signal-to-noise ratio. 

The present noise level of this sensor is ±1 E.U. (1 a at an 
integration time of 10 sec) and is limited by background noise 
in the laboratory. Using this sensor, we recently carried out 
an experimental simulation where we measured in veal time 
gravity gradient fields that had exactly the same magnitude and 
time variation as the gravity gradient signals that would be 
expected from an orbiting vehicle around the moon. Bell [1970]. 

At the present time, the development effort on the sensor 
is heading in two different directions. One program, sponsored 
by the Air Force, is for the development of an airborne gravity 
gradient measurement system. The major thrust of the develop- 
ment effort is to design a suitable hard mounted bearing and 
drive system that will spin the sensor at the desired speed with- 
out introducing excessive amounts of noise and to design a vi- 
bration and rotation isolation system that will isolate the sensor 
system from the aircraft noise and motion. The goal of the pro- 
gram is to develop a gradiometer system capable of measuring 
gravitational gradients at the 1 E.U. level with a 10 sec inte- 
gration time on a moving base, such as an aircraft or submarine. 

The other program, under NASA sponsorship, is for the de- 
sign of an earth orbiting gravity gradient measurement system. 
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For the orbital case, the optimum method is to fabricate a 
sensor with a relatively low resonant frequency (2 to 8 Hz) 

(see Fig. 2) , attach it directly to the spacecraft and spin 
the spacecraft itself at the desired spin speed (1 to 4 rev/sec), 
(see Fig. 4.) This mode of operation has two significant ad- 
vantages. There are no bearing noise problems that are the 
primary source of difficulty in earthbound operation, and most 
importantly, since the spacecraft is rotating along with the 
sensor, the gravity gradient field of the spacecraft is station- 
ary in the frame of reference of the sensor and the sensor does 
not sense the gravity field of the spacecraft, on?y the gravity 
gradient field of the earth. The objective of the program is 
to develop a sensor system capable of measuring gravitational 
gradients at the 0.01 E.U. level with a 30 sec integration time. 

The effective resolution of a gradiometer at an altitude 
of 250 km is approximately 250 km. If a gradiometer were placed 
in a near polar orbit with suitably chosen orbital parameters, 
it would pass within 250 km of every point on the earth in 80 
orbits, thus completely mapping the earth in 5 days. 

C. Geodesy with Orbiting Gradiometers 

The application of orbiting gravity gradient sensors to 
geodesy is straightforward. An objective of geodesy is to 
determine the variations of the earth's gravitational potential, 
which can be expressed in terms of spherical harmonics: 

OO 00 

1 + L, ‘r’" £. P nm (,ln «" lc n»°°* " x + S nn Bin " X, l 
n**z m-u 
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where a is the mean radius of the earth, P „ is the normal- 

nin 

ized Legendre polynomial, C and S m are the coefficients of 

nm nm 

the harmonic terms, and (r,tf,X) are the coordinate positions 
of the instrument. 

In the present satellite geodesy programs, orbital per- 
♦ urbation methods of obtaining the gravitational potential 
harmonics have led to the determination of the harmonics through 
the fourteenth degree and order, Rapp [1968]. In theory, this 
technique can be extended to obtain all higher orders of the 
gravitational potential; however, it is anticipated that it 
will be difficult to obtain the higher order components, 

Kaula [1969] . 

The advantage of gradiometer techniques in obtaining the 
higher order harmonics of the earth's gravitational field is 
straightforward. Terms with increasing n correspond to small 
scale features on or near the surface. Although the contribu- 
tion of these harmonic components to the gravitational potential 
is quite small, their contribution to the gravitational force 
gradient at a point above them is a substantial fraction of the 
gravitational gradient of the entire earth. 

To illustrate the behavior of the gravitational force 
gradient, let us examine the gradient which is predicted for 
higher orders of n. A typical term in the gravitational po- 
tential 

V nm ' r 5 (f )" +1 p nn (,,in * )c m eot nX <4 > 

gives rise to a radial gravity of 
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g r “ 17 “ (n+1) (?) P nm C nm COS mX (5) 

3 

and a radial gravity gradient of 

2 n+3 

r = = (n+1) (n+2) ^ (f ) P nm C cos mX . (6) 

The present technique for measurement of the gravitational field 
from orbit utilizes doppler velocity tracking of the orbiting 
vehicle, either from the ground or from other spacecraft. The 
portion of the differential doppler velocity due to the higher 
orders of the gravitatior. \1 field is given by the time integral 
of the acceleration field. For the radial doppler velocity we 
take the time integral of the radial acceleration 

Av -fq dt » ~ (^) P C sin mX (7) 

r J n av \r/ nm nm 

where we have used the fact that the maximum spatial periodic 
variation (m»n) has a time variation due to the orbital veloc- 
ity v given by 

co8 max mX " 008 nX * c08 ( n a ) “ 008 *) * 

If we assume that the strength of the components follows 
the statistical law ? nm ^ C nm ^ 10**Vn 2 Kaula [1968] , and that 
(2n+l) terms contribute to each order, we can calculate the 
doppler velocity, gravity and gravity gradient as a f v.uion 
of die harmonic order. These are plotted in Fig. S for 250 km 
altitude . 

The doppler velocity data in Fig. 5 are correct, although 
they differ by two orders of magnitude from what would be 
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calculated from Fig. 5-7, page 5-28 of the Williams town report, 

Kaula [1969] as presently published. In recent correspondence, 
William M. Kaula has brought attention to tne fact that the 

right hand ordinate of Fig. 5-7 in the Williamstown report should 

-2 

read 10 mm/sec rather than mm/sec. 

Figure 5 indicates that if satellite-to-satellite doppler track- 
ing techniques attain their anticipated sensitivity level of 
0.05 mm/sec at 30 sec, doppler tracking will be able to extract 
gravity data up to degree 50, and if a gravity gradiometer with 
an 0.01 E.U. sensitivity at 30 sec can be flown it will contribute 
significant information out to degree 75. The comparative signal- 
to-noise of the two techniques crosses over at degree 35. We thus 
see that tne two techniques are complementary ratner than competitive 
since below degree 35 doppler tracking has a better signal level 
while above degree 35 the gradiometer gives better data. 

The average strength of the higher order gravity variations 
predicted in Fig. 5 and Fig. 5-7 of the Williamstown report use 
a statistical model based on the autocovariance analysis of a 
large variety of samples of gravimetry, Kaula (1963, p. 524) and 
are estimated to be correct within ± 30%. A statistical model 
assumes that the phases of the various harmonics are not correlated, 
whereas we might expect some correlation in phases to occur at the 

position of significant geophysical anomalies, such as mountain 

* 

ranges. To obtain some feeling for this possibility, we have also 
looked at tne gravity fields to be expected at altitude for 
reasonable mass anomalies on the surface. 

D. Periodic Mass Models 

In an attempt to study further the relative sensitivity of 
doppler velocity tracking measurements and gravity gradiometer 
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measurements for the higher order gravity fields, a massive 
disc mass model was used to generate gravity data and the sig- 
nals expected for both a single disc and a periodic array of 
discs were calculated. 

In the single disc model we chose a disc radius of 150 km 
or disc diameter of 300 km. The disc mass was chosen so that 
the gravity at the surface was 10 mgal. This particular choice 
of mass is not important to the question of the relative sen- 
sitivity since the gravity, gravity gradient and doppler velocity 
signals all vary directly as the mass, and if the real anomaly 
is down an order of magnitude all the curves should be changed 
by 10. 

A plot of the results is shown in Fig. 6, which indicates 
that a disc with diameter 300 km, thickness 15 km, mass of 
1.7 x 10 16 kg and density difference of 0.015 gm/cc will create 
at an altitude of 250 km the following signals: 

• Vertical gravity of 1.5 mgal peak 

• Vertical gravity gradient variation of 0.11 E.U. 

• Vertical doppler velocity shift of 1.0 xm/nec. 

The analysis of a single disc is, however, not a close 
analogy to the periodic variation in the gravity field that is 
implied by the usual harmonic representation of the field. The 
disc model analysis was therefore expanded to a calculation of 
the signals expected over a periodic array of positive and neg- 
ative disc anomalies. The mass (positive or negative) was 
assumed the same as in the single disc analysis. The curves 
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in Fig. 7 are extrapolations of the data from the center portion 
of the disc array to eliminate end effects. The periodic sig- 
nals were 

• Vertical gravity ±0.65 mgal 

• Vertical gravity gradient ±0.1 E.U. 

• Vertical doppler velocity shift ±0.08 mm/sec. 

We notice that the gravity gradient magnitude is almost 
the same as for the single disc. This is because the gravity 
gradient signal, being the spatial derivative of the accelera- 
tion, has a sharp cutoff, and the signal from an adjacent disc 
of opposite mass actually contributes slightly to the total 
signal. The magnitude of the periodic vertical gravity signal 
has decreased slightly from the single disc signal. The broad 
signature of the vertical gravity signal causes signals from 
adjacent discs of opposite sign to partially cancel. 

Finally, notice the very large decrea w over an order of 
magnitude, in the vertical doppler velocity signal from the 
single disc case to the periodic disc case. This is because 
the doppler velocity signal is the integral of the acceleration 
signal and the integration tends to smooth out the periodic 
variation that we are looking for. 

E. Gradiometer Noise Limit 

The fundamental sensitivity of any sensor is determined 
by the thermal noise limitation. For the past ten years, we 
have been developing gravitational sensors working near their 
thermal noise limit. Forward [1965a] , For ward f 19 681 . Because 
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Fig. 7. Gravity Signals at 250 km Altitude from 
Periodic Disc Array. 
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this basic limit is dependent upon energy considerations, its 
calculation depends only upon very general parameters of the 
sensor, such as its temperature, mass, effective length, and 
time of integration. The results can then be applied to all 
gravity gradient sensors, regardless of their detailed design. 

In our torsional sensor the thermal signal-to-noise energy 
ratio is obtained by comparing the gravitational gradient sig- 
nal energy stored in the sensor to the kT of thermal noise in 
the resonant mode, Be 11 [1970]. 

The minimum gradient that can be measured for a thermally 
limited sensor, (S/N = 1) is 

r = P = 4(^) V2 ' 

where t = Q/oi is the 1/e time constant of the sensor. 

For what might be the desired sensor for earth geodesy, 
one with a total arm mass of 4 kg, and an arm radius of 30 cm, 
the thermal noise equation gives us a kT limit of 0.005 E.U. 
for a 30 sec integration time. We should be able to reduce the 
measured noise below 2 kT with a properly designed structure 
and electronic matching circuit, Forward [1965a]. 

II. CONCLUSIONS 

A new instrument for remote sensing of the earth's gravity 
field, the rotating gravity gradiometer, has been demonstrated 
in the laboratory. Analysis has shown that the sensor could 
provide significant improvement in the determination of the 
higher harmonics and local anomalies of the earth's gravity field. 
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Development of orbital and airborne versions are under way and 
should result in the near future in gravity measurement systems 
that will have a significant impact on geodesy and earth physics, 
navigation and o T : it determination, and earth resources. 
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INTRODUCTION 


I . 

An ideal second order gravity gradiometer operating 
in an inertially fixed frame of reference will produce out- 
put signals proportional to the second order gradient of the 
gravity field only. The output of a nonideal sensor operat- 
ing in a noninertial frame of reference will contain, among 
other noises, errors due to acceleration and rotation of 
the sensor frame of reference. The purpose of this report 
is to present a dynamic analysis of this latter class of er- 
rors using a simple model of the second order gradiometer, 
and then to apply the results of this analysis to the 
specific case of a satellite MASCON* experiment. 


One such experiment is described in HRL Proposal 70M-0691/ 
C2395 Preliminary Proposal for the Apollo Gravity Gradimeter 
Sensor . 
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SENSOR ANALYSIS 


1 1 . 

A. Gradiometer Mode l 

A simple nechanical model, consisting of three tor- 
sional springs and dampers and two inertia arms, will serve 
as the basis for this analysis. As shown in Fig. 1, the 
inertia arms are individually coupled to the case of the 
instrument through the supporting springs K^,K 2 and are 
mutually coupled by the common torsional spring K Q . Each 
spring is assumed to contain viscous damping defined by the 
coefficients D 0 ' D ]/ D 2 anc * *- s assume ^ to be infinitely rigid 
in all directions other than about the common torsional 
axis JC. 

107 1- I 



Fig. 1. Gradiometer Model. 

The inertial tensors of the arms are defined in terms 
of the principal axes of the arms as in eqs. (1) and (2) . 

I * 1 1 + iclcl (1) 

? 2 - i 2 l j 2 i 2 * • (2) 

The unit vectors 1^, J 2 are nominally orthogonal, but they 
become nonorthogonal because of the differential motion of 
the arms. In this analysis, the inertial tensors of (1) and 
(2) will be approximated by (3) and (4), wherein the tensors 
are described in a senaor case-fixed frame TJlc, and the pro- 
ducts of inertia are neglected. 
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[XT + kk] 


(3) 


[33 + kk] (4) 

The common torsional spring Kg also includes a strain 
transducer to sense the differential angular deflection of 
the inertia arms. The electrical output of this strain trans- 
ducer is amplified, phase-sensitive demodulated, and filtered 
to produce the ultimate output of the gradiometer. 

B . Equations of Motion 

Because the ultimate output of the gradiometer is a 
function of the differential angular deflection of the iner- 
tia arms about the common torsional axis, it is of interest 
to derive the dynamic equations of motion which relate this 
parameter to the sensible inputs to the instrument. This 
derivation will be accomplished on the basis of the classi- 
cal equivalence of torque to the time rate of change of the 
angular momentum of each arm. 

The angular momenta of the arms are defined in terms 
of their inertial tensors and angular velocities as (5) and 
(6) , where w is the inertial angular velocity of the sensor 

case, and the scalars 0,.9 are the velocities of the arms 

1 ' 2 

relative to the sensor case . 

^ = % x * [X + k0 1 ] (5) 

H, M 2 * [w + E0 2 ] (6) 

The time rates of change of these momenta are produced 
by external torques acting on the arms as the result of the 
gravity gradient, arm mass unbalance, and the elastic and 
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viscous coupling between the arms and between the arms and 
the case, as shown in (7) and (8) . 


H = M , + M, , + M, . + M , 
1 gl yl 2,1 cl 


H 2 = M g2 + "u2 + M l,2 * M c2 


For our present objectives, we require only the compo- 
nents of these momental rates about the sensor output axis 
k, and from (5) and (6) these scalars may be expressed as (9) 
and (10) . 


k * H 1 = IjO^ + k * [$ 1 • io + uj x(4> 1 • to)] (9) 

• ** _ . _ 
k • H 2 = *2^2 + ^ * ^2 * 141 + 10 x (^2 * ***) 1 • (10 


The elastic and viscous coupling torques acting on the arms 
may be expressed as (11) and (12) . 


k . 

■ 

lM cl 

+ M l,2 ! * - 

[A x + a 0 ] e 1 + a 0 b 2 

(11) 

1 k • 

[ “c2 

+ »2,1> * - 

lA 2 + V e 2 + Vl 

(12) 


where, by definition: 


A o = D c s + K o 


A 1 “ D l s f K 1 


A 2 - D 2 s + K 2 


The dynamic equations of arm motion may be stated in the 
matrix form (14) by combining equations (7) through (12) . 




i 


(I i s + A i + V 


t-V 


("Ag) 


(I 2 S + A 2 + VJ L e 2 


where, by definition: 


"l “ k ' (M gl + M ul ‘ <*1 ' “ 


+ u) x ) i 


m 2 = k • (M g2 + M y2 - (* 2 * w + w X (5 2 *w))] I 

The solution of (14) for the differential angular deflection 
of the arms is stated as (16) . 

(I 2 S 2 + A 2 ) M l - (IjS 2 + 

1 2 I 1 I 2 S 4 + (I 1 (A 2 ♦ A 0 ) + I 2 (A 1 + A 0 )]S 2 + A x A 2 + A 0 (A 1 + Aj) 


A normalized form of (16) is presented as (17) , wherein the 
denominator is factored into two quadratics representing the 
SUM MODE and the DIFFERENCE MODE characteristic frequencies 


0 1,2 and V 


(e 1 - e 2 ) 


1 2 ^1 1 2 2 t ^2 

(S 2 + o 2 S + 6 2 1 j=-\ - [S 2 + S + 6p 

cs 2 + a i,2 s +*l,2 ] \ s2 + ir s + w o) 


where 
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a l 

- V 1 ! 


- Wl 

a 2 

D 2 /,I 2 

S 2 

“ K 2 /:t 2 


D 1 + D 2 

<*2 

K i + K 2 

12 

+ *2 

*1 + X 2 


"i, + I ' 


D, + D " 

1 2 


1 2 

1 1 I 2 

1 

L 1 ! + X 2j 


I, + 1^" 


’k. + K ' 

1 2 

y 4 

1 2 

*1*2 

» m 

K o 

h + *2 

m m 


The sum mode parameters a 2 2 ^1 2 are a PP rox i matei -Y equal to 
their counterparts in the numerator of (17) and become 
exactly equal to them when the spring rates , K 2 » damping 
D^, and inertias 1^, I 2 are matched. Under these ideal 
conditions, the difference mode frequency may be expressed 
as (18) . 


wjj = (2K Q + K)/I (18) 

where K = = K ? and 1 i ^ ■ 1 2 » 

The differential angular deflection 8. - 6., may be 

2 1 * 
scaled by Wq/Q to yield an expression for an equivalent 

gradient signal, as shown in (19) . 

r e " (u, 0 /Q) l6 l " V (19) 

Combining (17) and (19) yields (20) . 
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(S 2 + a 2 S + 6 2 ) 


ft)- 


(S + S + & 


> ©1 


(s2 + “ 1 , 2 s ♦ <2 ( s2 + r s + “o) 


( 20 ) 


At this point it is instructive to expand and normalize (15) 
to obtain the forcing functions of (20) . The result of this 
expansion is (21) . 


M, 


M, 


- r . . - (i 


+ ui . w . + r , ) 

k i D Ul | 


j. . . — to. - to . w . + r 

l] k y2 


where 


( 21 ) 


r „l * (k • V' 1 ! 


r „2 - (k ’ V 7I 2 


The largest gradient error due to arm mass unbalance occurs 
when the individual unbalances are of opposite sign. There- 
fore, in this analysis, the largest mass unbalance error will 
be assumed to be as defined by (22) . 


r 


yi 


A 



r 

ymax 


( 22 ) 


The resulting equivalent gradient is obtained by combining 
(20) , (21) , and (22) and is presented as (23) . 
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Equation (23) illustrates the following three basic gradient 
error sources . 


1. Rotational field (2u).w.) 

1 D 

2. Arm mass unbalance ( 2 r ) 

ymax 

2 2 

3. Sum Mode mismatch (3^ - * 

This equation also demonstrates a potentially large error 
due to excitation of the sum mode frequency 8. To avoid 

1 / A 

this difficulty, care should be taken in the design of the 
instrument for a particular application such as to separate 
the sum mode frequency from the characteristic motion fre- 
quencies of the application. 

Another property of (23) is that errors due to sum 
mode mismatch are attenuated as the square of frequency 
above the sum mode frequency up to the difference mode fre- 
quency u)q. Above the differexxce mode frequency these errors 
are attenuated as the fourth power of frequency until the 
lead of the numerator becomes effective. For light damping 
of the sum mode, the lead becomes effective at u> * 2Q^ 2 

This latter frequency is usually above the difference mode 
frequency by two orders of magnitude or better; beyond that 
frequency the error is attenuated further as the third power 
of frequency. 
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III. SATELLITE SENSOR INTERACTION 

A . Satellite Application 

In order to apply the general result of the gradiometer 
dynamic error analysis to the specific application of the 
satellite MASCON experiment, it is necessary to estimate 
the acceleration and rotation environment in which the sen- 
sor will be operating. This objective may be accomplished 
logically by estimating the satellite motions in its princi- 
pal inertial axes and then transferring the effects of these 
motions to the sensor by considering the orientation of the 
sensor with respect to the satellite's principal inertial 
axes and center of mass. 

B . Satellite Equations of Motion 

The equations of motion of the satellite developed 
here will be based on the assumption of a lunar circular 
orbit. A set of local orbital coordinates XYZ are defined 
such that Z is normal to the orbital plane and such that Y 
coincides with the local lunar gravity vector. The remain- 
ing unit vector X is defined to be normal to Y and Z and 
tangent to the orbital path as shown in Fig. 2. The princi- 
pal axes of the satellite xyz are defined to be rotated with 
respect to the local orbital frame XYZ by the classical 
Eulerian angles as defined by (24) . 

xyz = m 3 f0J x I*J 3 XYZ (24) 

The numerical subscripts of the rotation matrices indicate 
the axes about which the rotation is taken. The equations 
of motion will be based on the classical equivalence of 
torque to time rate of change of angular momentum as defined 
by (25) . 
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(25) 


H s = M r + M d . 

The expansion of (25) in terms of the principal inertial ten- 
sor of the subsatellite ^ and the gravity gradient torque 
Mj. is defined by (26), where other disturbance torques are 
neglected . 

<t> • a) + <ii x [<fc •<*>]= [Y x ($> • Y) ] (26) 

S s u s 

where 



xx A + yy B + zz C 
GM /rJ 

Universal gravitation constant 
Lunar mass 

Mass center distance between moon and 
satellite 

Inertial angular velocity of the satellite. 


A simplified expansion of (26) is given as (27) where the 
nodal angle <)> is assumed to be zero and the inclination angle 
6 is assumed to be small. 

* ort-t 



Fig. 2. Lunar Orbital Path Coordinates. 
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i) + k, 10 id 

x 1 y z 


- 3ft-9k, cos a) t 
0 1 s 


(jj - k-w id 

y 2xz 


Sft^ek- sin id t 
0 2 s 


“z + k 3“x“y 


1 £l!:k- sin 2uj_t 

2 U 3 S 


(27) 


where 


k x = (C - B)/A 


ip = U) t 
s 


k 2 = (C - A)/B 


k 3 = (B - A)/C 


A linear solution for the transverse rates u , u of (27) may 

x y 

be obtained through the assumptions that the subsatellite 

spin axis velocity is constant u> = w and that the trans- 

z s 

verse inertia ratios k^, k 2 are equal. The Laplace trans- 
formation of these linearized differential equations is 
stated in the matrix form of (78) . 


s p 


u (S) 

X 


-p s 


“ y (s) 



to x (0) - 3ftJ6k 


u)y (0) + 3ft q 0k 


2 2 
S + u >“ 
s 


u> 


s 

i : 2 


S + w 


S ) -> 


(28) 


where 


k ^ k l ^ k 2 


P & 
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The solution of (28) is stated as (29) . 


b) (S) 
x 


uj (0) S - pw (0) 

x ~ _ Jr. 

2T 


[s‘ + p j 


[3flQ0k] [S 2 + pto ] 

— 2 2 2 2 ~ 

[S* + p^) [S^ + to‘] 


to (S) 

y 


to (0) S + pw (0) 

X X 

5 

[s^ + p^l 


[3ftj0k] [to g - p] S 

+ — 2 5 — ^ T~ 

[S* + p Z ] [S* + to g ] 


(29) 


The time solution of (29) is stated as (30) . 


to ~ to (0) cos pt - 

X X r 


CO (0) 

y 


3ft 2 8k 

U) s~ + "P 


sin pt 


•( 


^o 9 * 
*8 + P 


to * to (0) sin pt + 

y * 


w y (0) 


3ft 2 0k 
“s + p 


cos pt - 


3ft 2 0k 


(30) 


In order to simplify expansion of the gradient error 
terms, (30) will be simplified as shown in (31) . 

w x “ fl p 008 pt “ n s 8in u) s t 
to.. * (J sin pt - ft. cos <0 t 

y p • * 


(0 * to) 

2 S 


(31) 


where 


sin to. t 

8 


COS 10 t 
s 
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« g = [3k6iig]/ [u g + p] 

<yo> » - H g . 


Substitution of (31) into (27) yields the approximate angular 

accelerations of the satellite as (32) , where w = id and 
i o z s 


U) 


X 



sin pt - u> ft 
s s 


COS Gd t 
S 


U) 


U3 


z 


- pQ cos pt + u> £1 sin a) t 
p s s s 

^3 2 2 

- x— [ sin 2w t - ft sin 2pt + 2ft ft cos (u> 
2 I, 0 s p r ps s 


p) 1 1 . 

(32) 


C. Gradiometer Reference Motion 


The gradiometer case is assumed to be fixed in the 
satellite but that its reference axes are misaligned with 
respect to the principal axes of the satellite and that its 
center of mass is displaced from that of the satellite. The 
angular misalignment of the gradiometer is defined by the 
small, transverse angles a, 6 shown in matrix form as (33) . 


1 0 -a 


ITJc A 


0 


1 


la -0 



(33) 


The gradiometer center of mass is located with respect to 
that of the satellite by the position vector T. This loca- 
tion of the gradiometer results in the translational accel- 
eration defined by (34) . 
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a = tii x A + u x [u xi] 


(34) 


The angular velocities and accelerations of the 
gradiometer may be described in terms of those of the satel- 
lite through the transformation of (33), as shown in (35) 
and ( 36) . 


W . - U) - QUO 

IX z 


a) . - to + 8w 

3 y z 


w, - w + aw — Bw 
k z x y 


(35] 


u) . - to -aw 

IX z 


u,j = U) y + 0u> z 


w k = “z + au) x " 6 “y 


(36) 


D. Expansion of the Gradient Error Driving Functions 


The gradient error driving functions are defined by 

(23) as 2u).w., w. , and 2r • Each of these driving func- 

13 k umax 

tions will be expanded in the time domain, and the ampli- 
tudes and frequencies so derives will be tabulated. 

T> rotational field driving function 2w^a)j from (35) 
may be expressed as (37) 


2u)jWj ■ 

Expansion of (37) 


2w x u»y + 2u> z (Bw x - aw y ) 
from (31) yields (38). 


- 2aBw 


2 

z 


(37) 
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2 

2u> . u> . = 2uM2_ tB cos pt - a sin pt] - 2aBu> 

1 j s p s 

+ ft 2 sin 2pt - 2ft ft cos (w - p) t 
p c p s s c 

2 

+ 2u>' ft„ (a cos a) t - B sin w t) + ft sin 2w t (38) 
s s s s s s 

The amplitudes and frequencies of the rotational field driv- 
ing function are listed in Table J. 


TABLE 1 

ROTATIONAL FIELD DRIVING FUNCTIONS 


Frequency 

Amplitude 

dc 

2 

2aBw 

s 

p 

2w ft /a 2 + B 2 
s p 

2p 

ft 2 

P 

w s " P 

2 Vs 

u s 

2w ft /a 2 + 

8 8 

2u) ft 



tT2T 


The sum mode mismatch driving function from (32) 
and (36) may be expanded as (39) . 

* - olpftp sin pt ♦ w # Q # cos w g tj 

- Btpftp cos pt ♦ sin w g t) 

lc 

♦ ^ (3Qq sin 2w g t - Q 2 sin 2 pt ♦ 2Q p Q g 


cos (w g - p)tl 
( 39 ) 
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The amplitudes and frequencies of the sum mode mismatch driv- 
ing function are listed in Table II. 

TABLE II 

SUM MODE MISMATCH DRIVING FUNCTIONS 


Frequency 

Amplitude 

P 

pflp /? + 6 2 

2p 

— k n 2 
2 K 3 P 

u s ' P 

k,n n 

3 p s 

w s 

ui A 2 + 

s s 

2w s 

7 *3 0 


T226 


The arm mass unbalance driving function (2T ) is 

pnici x 

produced by the action of the gradiometer acceleration a 
acting on the individual mass unbalance of each arm. The 
torques developed about the gradiometer output axis are de- 
fined by (40) in terms of the individual arm mass unbalances 
mT^ f . 


M ui - E ’ tm^x a) 

M y2 » Jc • (nTjXa) (40) 

An approximation to the arm mass unbalance driving function 
is stated as (41), where ^ - 7 2 ^ T. 

* I 2 ■ 17x5,1 • (41 > 
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* 


The expansion of (41) from (34) in the gradiometer reference 
frame is stated as (42) . 


2r pmax ’ I 2 + “j<*Vk> + + Vj” 

+ ~ [UJ?( -6 i A ^) + ^ j (« j^ ± ) + w k (6 j £ i ' 

+ I 2 “V’j' 5 ! 1 ! - W + “i"k‘- { )*k> + “j"k (, iV ' 


(42) 


We are interested in estimating an upper bound for 
the mass unbalance driving function. In the interest of 
simplicity, an extreme upper bound may be defined from (42) 
by neglecting the algebraic signs of the coefficients. This 
inequality is stated as (43), where (m6JL/I) = e and (1../Z.) 

a K 11 K 


= P. 


2 ^ Umax S 2e [ “i T “j + <“i + “j> “k 1 


+ 2ep t2Wj t + w 2 + + 2u> 2 + 2w^<u. 1 (43) 


Because of its complexity, the time function of (43) will 

not be stated here. Father, an estimate of the most sig- I 

nificant terms is presented as Table II I, whe re the mis- 
alignment angle y is app.oximately /a 2 + 8 2 . The dimen- t , 

sionless mass unbalance parameter e » [mi. .. /I] is estimated 1 

-4 -5 K ! 

t to be 10 £ e 1 10 . TI.«» mass unbalance driving functions V i 

are much smaller than the most significant terms of the ro- * 

tational field and sum mode mismatch functions. m 

I 

- i: 

s i 

t ■ 
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TABLE III 

MASS UNBALANCE DRIVING FUNCTION 


Frequency 

Amplitude 

dc 

2 

4e ( y + p)u> 

p 

4eu> ft 

s P 

2p 

2e (y + p) ftp 

w s 1 p 

4e ( y + p)ftpft s 

w s 

4goj ft 
s s 

2u> 

s 

3ek 3 (Y + 2p)ft 3 
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The most significant gradient errors have been deter- 
mined by substitution of the tabulated driving functions into 
(23) . These gradient errors are listed in Table IV according 
to frequency. The magnitudes of these errors have been cal- 
culated for a set of representative parameters which are 

roughly the same as those used by DeBra in his preliminary 

* 

analysis of this problem. These parameters are defined 
by the following. 


w. 


s 

p = 
ft* = 


(6 1,2 /2 “s> 


1 Hz 

0 

=s 

0 . 3 rad 

0.1 Hz 

Y 

= 

0.2 rad 

10" 3 sec -1 

op 

k 3 

s 

-3 

10 sec 

1.3 x 10 -8 sec~^ 

s 

0.1 

CM 

• 

O 

n 

3S 

0.05 


Q 

S 

30 


-1 


Private communication with Dr. Daniel DeBra, Stanford 
University. 


287 





l 


TABLE IV 

REPRESENTATIVE GRADIENT ERRORS 


Frequency 

Error Function 

Error (E.U.) 

dc 

2 2 ,^ 
Y |*> S /Q 


5.3 x 10 7 

P 

2 > w sV Q 


8.4 x 10 4 

2p 

n p/ Q 


33 

w s 1 p 

2 n a /q 

p s' * 


1 

o 

r-4 


2Yu s a s /Q 


1.1 


„ re, 

2 


2w 

s 

i 1 - 5 ^] [ 2 i 


0.3 


L S J 
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It is of interest to observe that the most significant 
error at the tuned resonant frequency * 2w g is produced 

by the sum mode mismatch, and that this error is only 0.3 
E.U. for a 5% mismatch of the square of the sura mode fre- 

.22 j 

quencies 8^, $ 2 * The remaining errors in Table IV are all 
the result of the rotational field driving function 2u)£tiij, 
and the largest of these occur at dc and at the satellite nu- 
tation frequency p. 

To determine the effect of these error functions on the 
ultimate output of the grad iome ter , we must model and examine 

f 

the signal transmission properties of the gradiometer filter 
process. This process consists of the piezoelectric strain 
transducer and its amplifier matching circuitry, a phase- 
sensitive demodulator operating at the tuned resonant fre- 
quency w Q * 2w g , and a second order filter operating on the 
demodulator output. This process is illustrated in the block 


] 

] 

] 

] 

] 

1 

1 

1 

I 
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diagram of Fig. 3 wherein a simple, untuned model of the 
strain transducer is shown with its gain normalized to 
unity at the tuned resonant frequency u)q = 2u) g . The phase 
demodulation process acts as a frequency converter which 
produces the sums and differences of the input frequencies 
with the tuned resonant frequency u)q = 2u ' s * For example, 
input signals at the nutation frequency p will be converted 
to output signals at frequencies 2w ± p, and input signals 
at the tuned resonant frequency will be converted to dc and 
to the second harmonic of the tuned resonant frequency. 

This converted spectrum is then filtered by the second order 
post-demodulation filter to produce the output gradient sig- 


r. 



107 2 -* 


To 




Fig. 3. Gradi ometer Filter Process Block Diagram. 


As previously stated, the large error magnitudes of 
Table IV occur at dc and at the subsatellite nutation fre- 
quency p. The dc term is blocked by the action of the 
strain transducer and does not appear in the gradient output. 
The term at nutation frequency is attenuated by the strain 
transducer an amount dependent on the parameter cXq . The 

C Q , where C Q is the e.’.t-ctr ical capacity of the 
rand Rq is u its electrical load resistance. 


t?ansS 




****•&* *■*•»«, 


1 


i 


normalized strain transducer gain at nutation frequency (for 
this model) is stated as (44). 


p/l + U n /\0 2 

— . (44) 

/ 2 ~ 2 
/p + a Q 

The demodulator converts the nutation frequency to 2u ± p, 

s 

and the second order output filter provides an additional 
attenuation dependent on the time constant t. The gain of 
the output filter at 2w s t p is approximately the same as 
at ?'i' - This gain is stated as (45). 



V 2 “ s ) 


1 

1 + (2u T| 2 

b 


(45) 


The large nutation frequency error term of Table IV may be 
reduced to an acceptable output level by selection of the 

4 

parameters cx^, t. For example, this error term of 8.4 x 10 

E.U. may be reduced to approximately 1.5 E.U. at the output 

by selecting = 2u» g and t * 5 sec. The combined gain of 

(44) and (45) using these parameters with p * 0.1 w and 
—1 _c s 

u) = 2 it sec is approximately 1.8 x 10 . 

It should be emphasized that the filter estimates used 
here are conservative. An untuned model of the strain trans- 
ducer was employed, and a moderate value of the output filter 
time constant was selected. It seems entirely feasible to 
provide an additional order of magnitude of attenuation on 
the nutation frequency error function without significant 
degradation of the proposed gradient signal measurement. 


* 

s 
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IV. 


SUMMARY AND CONCLUSIONS 


A dynamic analysis of a simple model of the second 
order gravity gradiometer has been performed to determine 
the errors which result from acceleration and rotation of 
the gradiometer frame of reference. The results of this 
analysis were applied to the satellite MASCON experiment by 
estimation of the satellite motions when in lunar circular 
orbit . 

Numerical evaluation of the mot ion -induced gradient 
errors indicates one potentially large error due to the 
rotational field produced by subsatellite nutation and 
gradiometer misalignment. It appears quite feasible to re- 
duce this error to well below the 1 E.U. level by gradiometer 
filter design. 



